Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3020TG48.4% 105.5 / 47.0 97H219P (CAC→CCC) yghSputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/45);  new base G (47/0);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 2.94e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133156‑3133450
                                                                                                                                                   |                                                                                                                                                        
gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCg                                                                                                                                                      >  2:1354787/1‑151 (MQ=255)
 tAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgggtgg                                                                                                                                                      >  1:153325/1‑148 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                              <  1:476829/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                              <  2:152116/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                              <  1:1162412/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGgcccacca                                                                                                                                             >  2:254970/1‑146 (MQ=255)
         cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAg                                                                                                                                            <  2:29589/151‑1 (MQ=255)
         cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAg                                                                                                                                            <  1:56861/151‑1 (MQ=255)
            cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTc                                                                                                                                        <  1:1189925/151‑1 (MQ=255)
            cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTcccggt                                                                                                                                          >  1:197048/1‑146 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGCATTTGTAGTGCATCCATTCGTTAAATCATCCCCTATACATGTTGGTGATCAGGAAAGCGCGGCGTGAATGGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  1:1107021/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  2:1344852/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  2:1584608/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  1:1561610/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  2:1094069/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  2:6023/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                                                   <  2:1535488/151‑1 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCCTTCAAGGAt                                                                                                                                  >  2:691821/1‑150 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATTCCGTTCAAAGAt                                                                                                                                  >  1:1402440/1‑150 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGtttcagggt                                                                                                                                  >  1:455758/1‑143 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTTaaagga                                                                                                                                  >  2:919975/1‑146 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGAt                                                                                                                                  >  2:275531/1‑150 (MQ=255)
                      gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCATTTAAGGGTCttt                                                                                                                              >  1:489035/1‑151 (MQ=255)
                      gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATccccctttacggtcttt                                                                                                                               >  2:1721462/1‑136 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTTaaagggctttt                                                                                                                             >  2:546925/1‑141 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCCAGGATCtttt                                                                                                                             >  1:1201343/1‑149 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGGGAATTGCGCGTAGTGACGCTTGCAGAATTTCAACGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgc                                                                                                                            <  2:715505/151‑1 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTGGCCCATCCCGTTCAAGGAGCTTTtcg                                                                                                                           >  1:965945/1‑151 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCTAGGATCtttttg                                                                                                                           >  1:414816/1‑148 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTtcg                                                                                                                           >  1:31569/1‑151 (MQ=255)
                          tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcgc                                                                                                                          <  1:695908/151‑1 (MQ=255)
                          tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcgc                                                                                                                          <  2:726321/151‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGcc                                                                                                                         <  1:696973/151‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCCTCCCATTCaaagcactttttgg                                                                                                                          >  1:187677/1‑138 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCAGGTCCAGGATCtttttgc                                                                                                                          >  1:683656/1‑146 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGCTTAAgtgtattttcgc                                                                                                                          >  1:1126536/1‑139 (MQ=255)
                                cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTTCCGGTCCAGGAGCTTTTTGCCacttt                                                                                                                    >  2:725847/1‑147 (MQ=255)
                                  gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTTaaa                                                                                                                  >  2:279195/1‑150 (MQ=255)
                                  gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAAGATCTTTTCGCCATTGAAt                                                                                                                  >  1:702000/1‑151 (MQ=255)
                                   tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                     <  1:1565225/147‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  1:1493197/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  2:1157374/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  1:262757/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCAGCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  2:317629/151‑1 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTCCAAGATCTTTTCGCCATTTAATTTCa                                                                                                              >  2:606138/1‑151 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGGCa                                                                                                              >  2:645507/1‑151 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCACTTCAAGGATCTTTTCGCCCTTGAAttttag                                                                                                             >  2:891077/1‑146 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGTGCCATTGAATGTCAg                                                                                                             >  1:1529756/1‑151 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGGTCTTTGCGCCCTTGAATGttcg                                                                                                             >  2:433840/1‑148 (MQ=255)
                                            gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGGTCAAGGGTCTTTTCGCCCtttgattgcagtttcg                                                                                                        >  2:773210/1‑137 (MQ=255)
                                            gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGGTCTTTTCGCCATTGTATGttcgtttgg                                                                                                        >  2:990477/1‑143 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGTCCGTTggggcg                                                                                                     >  1:336990/1‑151 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCACCCAGTTCAAAGATATTTGCGGCCTTGAATTTCAGTTgcggcg                                                                                                     >  1:1495852/1‑151 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGAGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCATCCCGGTCCAGGAACTTTTGGCCCTTGAATGtttgtcgggggg                                                                                                     >  2:1647536/1‑140 (MQ=255)
                                                    tCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACg                                                                                                <  1:298214/151‑1 (MQ=255)
                                                     ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                               <  1:97265/151‑1 (MQ=255)
                                                     ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                               <  1:1349076/151‑1 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAAGGTCTTTTCGCCCTTTAATGTCAGTTGCGGGGGAAAGGCGAtt                                                                                         >  2:516839/1‑151 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCCCCCCTCCAGTTCCAAGATCTTTGCGCCCTTGGATGGCAGTTGCGGCGGTACGGCGGtt                                                                                         >  1:17318/1‑151 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCCTGCAGAATTTCATCGGCGGGGGGGGGTCCCAGCAGTTTAATGTTCTTTTCGGCCATGTAAGTCAGTTGTGGGGTTAAGGCgggt                                                                                         >  1:504073/1‑148 (MQ=255)
                                                            ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTTAAGGATCTTTTCGCCATTGGATGTTAGGTTCGGGGTGACGGCGGtttt                                                                                       >  1:289559/1‑151 (MQ=255)
                                                               aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTCAAGGATCTTTTTGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttttttt                                                                                    >  1:147261/1‑150 (MQ=255)
                                                                atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTcc                                                                                   <  2:356210/151‑1 (MQ=255)
                                                                atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCCAGGGTCTTTTCGCCCTTGAATGTCCGGTGCGGCGTAACGGCGGttttttt                                                                                    >  2:1391396/1‑149 (MQ=255)
                                                                 tcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCAACCAGGTCAAGGAGCTTTGCGCCATTGAATGTCAGGTGCGGCGTAAACGGGAttttttt                                                                                    >  1:1392026/1‑148 (MQ=255)
                                                                     ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGATCTTTTCGCCATTGAATGGCAGGTGCGGCGGAACGGCGGTTTTTTTCTTTaa                                                                              >  2:1326230/1‑151 (MQ=255)
                                                                       ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGGGTAACGGGGGTTTTTTTCTTTaaag                                                                            >  2:461386/1‑149 (MQ=255)
                                                                          aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTTGCCCTTTGATGTTAGTTGCGGCGGAACGGCGGTTTTTTTCTGttaag                                                                            >  1:1421224/1‑144 (MQ=255)
                                                                          aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTTAAGGATCTTTTTGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTTTTACg                                                                            >  1:1080261/1‑148 (MQ=255)
                                                                            ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCAGTTCCAGGGTCTTTTTGCCCTTTGATGTCCGTTTCGGCGTAAAGGCGGttttttttttttaccctgc                                                                        >  1:404657/1‑137 (MQ=255)
                                                                                      gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCCCCCCGGTCAAAGGGCTTTGCGCCCTTTTATGTCAGGTGCGGCGGGACAGCGGGTTTTTTCTTTTAAGGCTCCCAttttgggg                                                             >  1:1289812/1‑146 (MQ=255)
                                                                                      gATCAGGAAAGCGCGGCGTGAATTGCGCGCAGTGACGCTTGCAGAATTTCATCGGCGGCGGGGTCCCCCTCCCATTCAAAGGGCCTTTCGCCCTTGTATGTCTGGTTCGGGGGAAAGGGGGATTTTTTCTCTCACGCTGCCCGTtggtggg                                                            >  1:565691/1‑150 (MQ=255)
                                                                                          aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCCAGGATCTTTTTGCCCTTTGATTTTAGGTTCCGCGTTACGGCGGTTTTTTTCTTTTAAACTGCCACTttgtggggagg                                                        >  1:110979/1‑146 (MQ=255)
                                                                                           ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCACCCGTTTCAAGGTCTTTGCGCCCTTGAATGTCAGGTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTGCCAGTTGGTGGTCAggg                                                       >  1:1058478/1‑150 (MQ=255)
                                                                                           ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTGCGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTGCCCGTTGGGGGGCAggg                                                       >  2:1335762/1‑150 (MQ=255)
                                                                                            gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAgg                                                        <  1:279104/149‑1 (MQ=255)
                                                                                                 cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAcg                                                  <  1:1646703/150‑1 (MQ=255)
                                                                                                   cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAcgcg                                                <  1:989934/150‑1 (MQ=255)
                                                                                                   cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACgcgc                                               <  2:1653984/151‑1 (MQ=255)
                                                                                                     gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGGTCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTTACGCTGCCCGGTTGGGGGCAGGGTTTCGCGCaa                                             >  2:1078172/1‑151 (MQ=255)
                                                                                                             tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGtt                                     <  2:31586/151‑1 (MQ=255)
                                                                                                                   tAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGGCCATCCAGTTTAAGGATCTTTTCGCCCTTTAATGGCAGGTGCGGCGTAACGGGGATTTTTTTCTGTGACGCTGCCAGTTGGTGGGCAGGGTTACGCGCAAACGCGGgttgtt                                 >  1:1105354/1‑149 (MQ=255)
                                                                                                                   tAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtc                                <  1:981234/150‑1 (MQ=255)
                                                                                                                        aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  1:792337/150‑1 (MQ=255)
                                                                                                                         cGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTTCAGGGACTTTGTGCCATTGTAAGTCAGTTTTGGCGTAAAGGCGAGTTTTTTCTGTGACGCTGCCCGTTGGTGGGGCGGGTTACGCGCAAAACCGGGTTGTtcatttatt                         >  1:87424/1‑146 (MQ=255)
                                                                                                                           cTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGcc                       <  1:1161170/151‑1 (MQ=255)
                                                                                                                             tGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTCAAGGTCTTTGTGCCCTTTTATGGCAGTTGTGGCGCAAACGCGGTTTTTTTCTGTTATGCCGTCCCGTTGGGGTGCGGGTTTAGGGCCAAAACCGGTTTGTtctccaagccaa                      >  1:968661/1‑141 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGa                   <  2:979391/148‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:1287823/149‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:683988/149‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:1403101/149‑1 (MQ=255)
                                                                                                                                      tCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATca            <  1:589500/151‑1 (MQ=255)
                                                                                                                                        aTCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaac          <  1:275441/151‑1 (MQ=255)
                                                                                                                                          cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaacaa        <  2:1428451/151‑1 (MQ=255)
                                                                                                                                          cggcggGGTGTCGCCCATCCAGTTCAAGGATCTGTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaacaa        <  1:1720574/151‑1 (MQ=255)
                                                                                                                                            gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt       <  2:1072193/150‑1 (MQ=255)
                                                                                                                                             cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:1249175/151‑1 (MQ=255)
                                                                                                                                             cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:1330186/151‑1 (MQ=255)
                                                                                                                                                  gtgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  <  2:1290638/150‑1 (MQ=255)
                                                                                                                                                  gtgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  1:1285037/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133156‑3133450

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