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

GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133156‑3133450
                                                                                                                                                   |                                                                                                                                                        
gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCg                                                                                                                                                      >  2:1354787/1‑151 (MQ=255)
 tAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑Ggggtgg                                                                                                                                                      >  1:153325/1‑148 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATcc                                                                                                                                              <  1:476829/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATcc                                                                                                                                              <  2:152116/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATcc                                                                                                                                              <  1:1162412/150‑1 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑gcccacca                                                                                                                                             >  2:254970/1‑146 (MQ=255)
         cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAg                                                                                                                                            <  2:29589/151‑1 (MQ=255)
         cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAg                                                                                                                                            <  1:56861/151‑1 (MQ=255)
            cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTc                                                                                                                                        <  1:1189925/151‑1 (MQ=255)
            cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTcccggt                                                                                                                                          >  1:197048/1‑146 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGCATTTGTAGTGCATCCATTCGTTAAATCATCCCCTATACATGTTGGTGATCAGGAAAGCGCGGCGTGAATGGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  1:1107021/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  2:1344852/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  2:1584608/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  1:1561610/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  2:1094069/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  2:6023/151‑1 (MQ=255)
                 cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGa                                                                                                                                   <  2:1535488/151‑1 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCCC‑TTCAAGGAt                                                                                                                                  >  2:691821/1‑150 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATTCCG‑TTCAAAGAt                                                                                                                                  >  1:1402440/1‑150 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCCG‑tttcagggt                                                                                                                                  >  1:455758/1‑143 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCAG‑TTTaaagga                                                                                                                                  >  2:919975/1‑146 (MQ=255)
                   gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCAG‑TTCAAGGAt                                                                                                                                  >  2:275531/1‑150 (MQ=255)
                      gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCCA‑TTTAAGGGTCttt                                                                                                                              >  1:489035/1‑151 (MQ=255)
                      gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATccccctttacggtcttt                                                                                                                               >  2:1721462/1‑136 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTTG‑CCCCTCCCG‑TTTaaagggctttt                                                                                                                             >  2:546925/1‑141 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCAG‑TTCCAGGATCtttt                                                                                                                             >  1:1201343/1‑149 (MQ=255)
                        aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGGGAATTGCGCGTAGTGACGCTTGCAGAATTTCAACGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTgc                                                                                                                            <  2:715505/151‑1 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTGG‑CCCATCCCG‑TTCAAGGAGCTTTtcg                                                                                                                           >  1:965945/1‑151 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCCG‑TTCTAGGATCtttttg                                                                                                                           >  1:414816/1‑148 (MQ=255)
                         cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCTCCAG‑TTCAAGGATCTTTtcg                                                                                                                           >  1:31569/1‑151 (MQ=255)
                          tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTgcgc                                                                                                                          <  1:695908/151‑1 (MQ=255)
                          tATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTgcgc                                                                                                                          <  2:726321/151‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGcc                                                                                                                         <  1:696973/151‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGGCG‑CCCCTCCCA‑TTCaaagcactttttgg                                                                                                                          >  1:187677/1‑138 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTTG‑CCCCTCCAG‑GTCCAGGATCtttttgc                                                                                                                          >  1:683656/1‑146 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCATCCAG‑CTTAAgtgtattttcgc                                                                                                                          >  1:1126536/1‑139 (MQ=255)
                                cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTTCCG‑GTCCAGGAGCTTTTTGCCacttt                                                                                                                    >  2:725847/1‑147 (MQ=255)
                                  gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCCG‑TTCAAGGATCTTTTCGCCCTTTaaa                                                                                                                  >  2:279195/1‑150 (MQ=255)
                                  gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCCG‑TTCAAAGATCTTTTCGCCATTGAAt                                                                                                                  >  1:702000/1‑151 (MQ=255)
                                   tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTg                                                                                                                     <  1:1565225/147‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  1:1493197/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  2:1157374/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  1:262757/151‑1 (MQ=255)
                                    ttgagtttgTAGTGCATCCATTCGTTAAATCAGCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGt                                                                                                                <  2:317629/151‑1 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCTCCCG‑GTCCAAGATCTTTTCGCCATTTAATTTCa                                                                                                              >  2:606138/1‑151 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCTCCAG‑TTCAAGGATCTTTTCGCCATTGAATGGCa                                                                                                              >  2:645507/1‑151 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCCTCCAC‑TTCAAGGATCTTTTCGCCCTTGAAttttag                                                                                                             >  2:891077/1‑146 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCAAGGATCTTTGTGCCATTGAATGTCAg                                                                                                             >  1:1529756/1‑151 (MQ=255)
                                       agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCCG‑TTCAAGGGTCTTTGCGCCCTTGAATGttcg                                                                                                             >  2:433840/1‑148 (MQ=255)
                                            gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTTG‑CCCCTCCCG‑GTCAAGGGTCTTTTCGCCCtttgattgcagtttcg                                                                                                        >  2:773210/1‑137 (MQ=255)
                                            gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCCG‑TTCAAGGGTCTTTTCGCCATTGTATGttcgtttgg                                                                                                        >  2:990477/1‑143 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCCTCCAG‑TTCAAGGATCTTTTCGCCATTGAATGTCCGTTggggcg                                                                                                     >  1:336990/1‑151 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCACCCAG‑TTCAAAGATATTTGCGGCCTTGAATTTCAGTTgcggcg                                                                                                     >  1:1495852/1‑151 (MQ=255)
                                               gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGAGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGTG‑CCCATCCCG‑GTCCAGGAACTTTTGGCCCTTGAATGtttgtcgggggg                                                                                                     >  2:1647536/1‑140 (MQ=255)
                                                    tCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACg                                                                                                <  1:298214/151‑1 (MQ=255)
                                                     ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                               <  1:97265/151‑1 (MQ=255)
                                                     ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                               <  1:1349076/151‑1 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCAAAGGTCTTTTCGCCCTTTAATGTCAGTTGCGGGGGAAAGGC‑GAtt                                                                                         >  2:516839/1‑151 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCC‑CCCCTCCAG‑TTCCAAGATCTTTGCGCCCTTGGATGGCAGTTGCGGCGGTACGGC‑GGtt                                                                                         >  1:17318/1‑151 (MQ=255)
                                                          cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCCTGCAGAATTTCATCGGCGG‑GGGGGGG‑TCCCAGCAG‑TTTAATGTTCTTTTCGGCCATGTAAGTCAGTTGTGGGGTTAAGGC‑gggt                                                                                         >  1:504073/1‑148 (MQ=255)
                                                            ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCTCCAG‑TTTAAGGATCTTTTCGCCATTGGATGTTAGGTTCGGGGTGACGGC‑GGtttt                                                                                       >  1:289559/1‑151 (MQ=255)
                                                               aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCATCCAG‑TTCAAGGATCTTTTTGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GAttttttt                                                                                    >  1:147261/1‑150 (MQ=255)
                                                                atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTcc                                                                                   <  2:356210/151‑1 (MQ=255)
                                                                atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCCAGGGTCTTTTCGCCCTTGAATGTCCGGTGCGGCGTAACGGC‑GGttttttt                                                                                    >  2:1391396/1‑149 (MQ=255)
                                                                 tcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGTG‑CCCAACCAG‑GTCAAGGAGCTTTGCGCCATTGAATGTCAGGTGCGGCGTAAACGG‑GAttttttt                                                                                    >  1:1392026/1‑148 (MQ=255)
                                                                     ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTTAAGGATCTTTTCGCCATTGAATGGCAGGTGCGGCGGAACGGC‑GGTTTTTTTCTTTaa                                                                              >  2:1326230/1‑151 (MQ=255)
                                                                       ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGGGTAACGGG‑GGTTTTTTTCTTTaaag                                                                            >  2:461386/1‑149 (MQ=255)
                                                                          aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCAAGGATCTTTTTGCCCTTTGATGTTAGTTGCGGCGGAACGGC‑GGTTTTTTTCTGttaag                                                                            >  1:1421224/1‑144 (MQ=255)
                                                                          aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑GTTAAGGATCTTTTTGCCATTGAATGTCAGTTGCGGCGGAACGGC‑GATTTTTTTCTTTTACg                                                                            >  1:1080261/1‑148 (MQ=255)
                                                                            ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCCCCAG‑TTCCAGGGTCTTTTTGCCCTTTGATGTCCGTTTCGGCGTAAAGGC‑GGttttttttttttaccctgc                                                                        >  1:404657/1‑137 (MQ=255)
                                                                                      gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGTG‑CCCCCCCCG‑GTCAAAGGGCTTTGCGCCCTTTTATGTCAGGTGCGGCGGGACAGCGGGTTTTTTCTTTTAAGGCTC‑CCAttttgggg                                                             >  1:1289812/1‑146 (MQ=255)
                                                                                      gATCAGGAAAGCGCGGCGTGAATTGCGCGCAGTGACGCTTGCAGAATTTCATCGGCGGCGGGGTC‑‑CCCCTCCCA‑TTCAAAGGGCCTTTCGCCCTTGTATGTCTGGTTCGGGGGAAAGGG‑GGATTTTTTCTCTCACGCTG‑CCCGTtggtggg                                                            >  1:565691/1‑150 (MQ=255)
                                                                                          aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCCTCCCG‑TTCCAGGATCTTTTTGCCCTTTGATTTTAGGTTCCGCGTTACGGC‑GGTTTTTTTCTTTTAAACTG‑CCACTttgtggggagg                                                        >  1:110979/1‑146 (MQ=255)
                                                                                           ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCGCCCCACCCGT‑TTCAAGG‑TCTTTGCGCCCTTGAATGTCAGGTGCGGCGGAACGGC‑GGTTTTTTTCTGTAACGCTG‑CCAGTTGGTGGTCAggg                                                       >  1:1058478/1‑150 (MQ=255)
                                                                                           ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCATCCCG‑TTCAAGGATCTTTGCGCCCTTGAATGTCAGTTGCGGCGGAACGGC‑GGTTTTTTTCTGTAACGCTG‑CCCGTTGGGGGGCAggg                                                       >  2:1335762/1‑150 (MQ=255)
                                                                                            gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAgg                                                        <  1:279104/149‑1 (MQ=255)
                                                                                                 cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTAcg                                                  <  1:1646703/150‑1 (MQ=255)
                                                                                                   cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTAcgcg                                                <  1:989934/150‑1 (MQ=255)
                                                                                                   cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACgcgc                                               <  2:1653984/151‑1 (MQ=255)
                                                                                                     gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAG‑TTCAAGGGTCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGC‑GGTTTTTTTCTGTTACGCTG‑CCCGGTTGGGGGCAGGGTTTCGCGCaa                                             >  2:1078172/1‑151 (MQ=255)
                                                                                                             tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGtt                                     <  2:31586/151‑1 (MQ=255)
                                                                                                                   tAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑GCCATCCAG‑TTTAAGGATCTTTTCGCCCTTTAATGGCAGGTGCGGCGTAACGGG‑GATTTTTTTCTGTGACGCTG‑CCAGTTGGTGGGCAGGGTTACGCGCAAACGCGGgttgtt                                 >  1:1105354/1‑149 (MQ=255)
                                                                                                                   tAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtc                                <  1:981234/150‑1 (MQ=255)
                                                                                                                        aCGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  1:792337/150‑1 (MQ=255)
                                                                                                                         cGCTTGCAGAATTTCATCGGCGG‑GGGGGCG‑CCCATCCAG‑TTTCAGGGACTTTGTGCCATTGTAAGTCAGTTTTGGCGTAAAGGC‑GAGTTTTTTCTGTGACGCTG‑CCCGTTGGTGGGGCGGGTTACGCGCAAAACCGGGTTGTtcatttatt                         >  1:87424/1‑146 (MQ=255)
                                                                                                                           cTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGcc                       <  1:1161170/151‑1 (MQ=255)
                                                                                                                             tGCAGAATTTCATCGGCGG‑GGGGTCG‑CCCCTCCAGTTTCAAGG‑TCTTTGTGCCCTTTTATGGCAGTTGTGGCGCAAACGC‑GGTTTTTTTCTGTTATGCCGTCCCGTTGGGGTGCGGGTTTAGGGCCAAAACCGGTTTGTtctccaagccaa                      >  1:968661/1‑141 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGa                   <  2:979391/148‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:1287823/149‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:683988/149‑1 (MQ=255)
                                                                                                                                  aaTTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  2:1403101/149‑1 (MQ=255)
                                                                                                                                      tCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATca            <  1:589500/151‑1 (MQ=255)
                                                                                                                                        aTCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaac          <  1:275441/151‑1 (MQ=255)
                                                                                                                                          cggcgg‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaacaa        <  2:1428451/151‑1 (MQ=255)
                                                                                                                                          cggcgg‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTGTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaacaa        <  1:1720574/151‑1 (MQ=255)
                                                                                                                                            gcgg‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt       <  2:1072193/150‑1 (MQ=255)
                                                                                                                                             cgg‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:1249175/151‑1 (MQ=255)
                                                                                                                                             cgg‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:1330186/151‑1 (MQ=255)
                                                                                                                                                  gtgtCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  <  2:1290638/150‑1 (MQ=255)
                                                                                                                                                  gtgtCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  1:1285037/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGG‑GGTGTCG‑CCCATCCAG‑TTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGC‑GATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  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: ‑