Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,388 G→A 5.2% intergenic (‑89/‑97) ycbZ ← / → matP putative peptidase/Ter macrodomain organizer matS‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,3880GA 5.2% 248.2 / 2.4 77intergenic (‑89/‑97)ycbZ/matPputative peptidase/Ter macrodomain organizer matS‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base A (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTGGAAATGGA  >  NC_000913/1018242‑1018530
                                                                                                                                                  |                                                                                                                                              
ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTgg                                                                                                                                                <  2:82948/147‑1 (MQ=255)
      aGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATggg                                                                                                                                          <  4:351382/147‑1 (MQ=255)
          aTCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGtt                                                                                                                                      >  1:268310/1‑147 (MQ=255)
          aTCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGACATGGGCGtt                                                                                                                                      >  7:156976/1‑147 (MQ=255)
          aTCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGACATGGGCGtt                                                                                                                                      >  3:182669/1‑147 (MQ=255)
               tATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCt                                                                                                                                >  6:274745/1‑148 (MQ=255)
                   aGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTcaca                                                                                                                            <  5:274745/148‑1 (MQ=255)
                    ggAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAg                                                                                                                           >  3:336631/1‑148 (MQ=255)
                     gAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTc                                                                                                                               <  8:259583/143‑1 (MQ=255)
                     gAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTc                                                                                                                               >  7:259583/1‑143 (MQ=255)
                     gAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGACATGGGCGTTAAAGCTCACAgg                                                                                                                          <  8:156976/148‑1 (MQ=255)
                     gAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGACATGGGCGTTAAAGCTCACAgg                                                                                                                          <  4:182669/148‑1 (MQ=255)
                       aCCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGAt                                                                                                                        <  2:154174/148‑1 (MQ=255)
                           ggTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCTCTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGgag                                                                                                                    <  5:390840/148‑1 (MQ=255)
                             tCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTcac                                                                                                                             <  6:22534/137‑1 (MQ=255)
                             tCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTcac                                                                                                                             >  5:22534/1‑137 (MQ=255)
                                cGCCATGCAAGTTTCGTAATGCTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttct                                                                                                               <  8:321702/148‑1 (MQ=255)
                                  ccATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACAGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttcttt                                                                                                             >  2:442851/1‑148 (MQ=255)
                                      gCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttcttttctt                                                                                                         >  5:80198/1‑148 (MQ=255)
                                      gCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttcttttct                                                                                                          >  7:402181/1‑147 (MQ=255)
                                       cAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTttcttc                                                                                                        <  8:402181/148‑1 (MQ=255)
                                        aaGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCTGTACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTATCTTCa                                                                                                       <  5:85180/148‑1 (MQ=255)
                                            ttCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTttcttc                                                                                                        <  2:385359/143‑1 (MQ=255)
                                            ttCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTttcttc                                                                                                        >  1:385359/1‑143 (MQ=255)
                                            ttCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTgg                                                                                                   >  3:177873/1‑148 (MQ=255)
                                                 aaTGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCTCGATTATAACGTAACCGGAACGACAGCACAAGTAGAAATCACCGAATATACTACTTGCATGGTCTTTAAATATAACATGATGGATATTCTTTTCTTCAATGGACTa                                                                                               >  5:282955/1‑147 (MQ=255)
                                                    ggTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                                                                                      >  4:152472/1‑137 (MQ=255)
                                                    ggTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                                                                                      <  3:152472/137‑1 (MQ=255)
                                                    ggTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                                                                                      >  8:425991/1‑137 (MQ=255)
                                                    ggTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                                                                                      <  7:425991/137‑1 (MQ=255)
                                                      tagaaGTTGTTGTGTTTTAGTCTGTTGTCAAAGCGTCGATTATACCGTAACCGGCACTACAGCACACTTAGAATGCACCGACAATAATCCTGGCATTGGAGTTAAAGCTAACAGGATGGATATTATTTTATTCACTGGCCTAAAAAga                                                                                         <  6:282955/145‑2 (MQ=255)
                                                      tCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGc                                                                                         <  1:430511/148‑1 (MQ=255)
                                                       cAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTTCAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCt                                                                                        <  2:254088/148‑1 (MQ=255)
                                                               aTGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAATCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTc                                                                                 <  1:390699/147‑1 (MQ=255)
                                                                tGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTACCCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTg                                                                               <  5:432136/148‑1 (MQ=255)
                                                                             gttgtCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaaa                                                                                             >  7:223225/1‑121 (MQ=255)
                                                                             gttgtCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaaa                                                                                             <  8:223225/121‑1 (MQ=255)
                                                                             gttgtCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTacac                                                                  >  3:186672/1‑148 (MQ=255)
                                                                                gtCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttctttt                                                                                                            <  7:424859/103‑1 (MQ=255)
                                                                                gtCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttctttt                                                                                                            >  8:424859/1‑103 (MQ=255)
                                                                                 tCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGt                                                               >  4:28782/1‑147 (MQ=255)
                                                                                 tCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGt                                                               >  5:308783/1‑147 (MQ=255)
                                                                                 tCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGt                                                               >  1:44952/1‑147 (MQ=255)
                                                                                   aaaGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTttcttc                                                                                                        >  4:375535/1‑104 (MQ=255)
                                                                                   aaaGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTttcttc                                                                                                        <  3:375535/104‑1 (MQ=255)
                                                                                               ataCCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaaaaa                                                                                           >  4:406827/1‑105 (MQ=255)
                                                                                               ataCCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaaaaa                                                                                           <  3:406827/105‑1 (MQ=255)
                                                                                                taCCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCt                                                <  6:308783/147‑1 (MQ=255)
                                                                                                taCCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCt                                                <  2:268310/147‑1 (MQ=255)
                                                                                                  ccGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTAt                                              >  7:234216/1‑147 (MQ=255)
                                                                                                    gTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGa                                            >  7:246097/1‑147 (MQ=255)
                                                                                                       aCCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAata                                        <  3:28782/148‑1 (MQ=255)
                                                                                                        ccGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGt                                                                       <  6:243011/116‑1 (MQ=255)
                                                                                                        ccGGCACTACAGCACACGTAGAAAGCACCGACAATACTACTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGt                                                                       >  5:243011/1‑116 (MQ=255)
                                                                                                           gCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATcaa                                    <  4:177873/148‑1 (MQ=255)
                                                                                                                aCAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACAtt                                                         >  8:442215/1‑122 (MQ=255)
                                                                                                                aCAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACAtt                                                         <  7:442215/122‑1 (MQ=255)
                                                                                                                   gCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTg                             >  4:325748/1‑147 (MQ=255)
                                                                                                                         gTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCt                      >  6:297831/1‑148 (MQ=255)
                                                                                                                         gTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCt                      <  8:234216/148‑1 (MQ=255)
                                                                                                                          tAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCtt                     <  1:442851/148‑1 (MQ=255)
                                                                                                                          tAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCtt                     <  2:170077/148‑1 (MQ=255)
                                                                                                                           aGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCtt                     >  8:292442/1‑147 (MQ=255)
                                                                                                                            gAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGa                   <  6:80198/148‑1 (MQ=255)
                                                                                                                              aaGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGa                   >  6:21122/1‑146 (MQ=255)
                                                                                                                                gCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGc               <  3:325748/148‑1 (MQ=255)
                                                                                                                                   ccGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGt            <  4:336631/148‑1 (MQ=255)
                                                                                                                                   ccGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGt            <  4:186672/148‑1 (MQ=255)
                                                                                                                                   ccGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGt            <  8:246097/148‑1 (MQ=255)
                                                                                                                                   acGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGt            >  2:27515/2‑148 (MQ=255)
                                                                                                                                     gACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTg          >  4:330955/1‑148 (MQ=255)
                                                                                                                                      aCAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACt                                                                                                     >  1:361209/1‑56 (MQ=255)
                                                                                                                                      aCAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACt                                                                                                     <  2:361209/56‑1 (MQ=255)
                                                                                                                                        aaTACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTGGaa       <  7:292442/148‑1 (MQ=255)
                                                                                                                                        aaTACTCCAGGCCTTGGCGGTAAAGCTCACCGGAGGGATATTCTTTTCTTCCATGGCCTAACAAGCTGCTATTCCGGAAAGACGTACCCGGGAACATTGCGATCGCGTTGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTGGaa       <  6:405930/148‑1 (MQ=255)
                                                                                                                                              ccTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGAt                                                    <  5:215389/97‑1 (MQ=255)
                                                                                                                                              ccTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGAt                                                    >  6:215389/1‑97 (MQ=255)
                                                                                                                                               cTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTGGAAATGGa  <  5:297831/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTACACGGTAACATTGAGATCGCTATGAAATATCAACAACTTGAAAATCTTGAAAGCGGTTGGAAATGGA  >  NC_000913/1018242‑1018530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: