Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA40.3% 85.6 / 36.6 72pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (36/7);  new base A (0/29);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.46e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑1222382
                                                                                                                                                     |                                                                                                                                                      
gtgAAAGTCGAAGGAGATGGGGCGAGTAATGTGGGGGAAATTAAACTAGGCGTAAATGGCAGTCTGAATCCAGCGGCCAGCGTCTGGTGAATGTGGGCGTGCAGCGGGGTGATAATCGCTACAATGACACCGCAGTGGTGGTGGGCCTGaa                                                                                                                                                      <  2:726298/151‑1 (MQ=255)
gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                                      <  2:1246383/151‑1 (MQ=255)
   aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    >  1:465534/1‑150 (MQ=255)
        cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                              <  2:1236495/151‑1 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                             >  1:728319/1‑151 (MQ=255)
            cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                                          >  1:2200/1‑151 (MQ=255)
            cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGAAATTGGGCGGGCAGATTGGGGATTATGGGTAGAATAAACCCGCAGGGATGGGGGGCCTGGAattaaaggt                                                                                                                                             >  1:643957/1‑140 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                         <  1:1165169/151‑1 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                         <  2:91079/151‑1 (MQ=255)
              aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                                        >  2:315023/1‑151 (MQ=255)
                      cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGt                                                                                                                               >  2:852569/1‑151 (MQ=255)
                          tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                           >  2:920775/1‑151 (MQ=255)
                          tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                           >  2:487281/1‑151 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                                                                                                                        >  1:720605/1‑150 (MQ=255)
                                  gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAAc                                                                                                                   >  2:272961/1‑151 (MQ=255)
                                  gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGAGGGCCTGAAATATAAGTTCTGATCACGCCGTTAGCTaaaaa                                                                                                                     >  2:422715/1‑149 (MQ=255)
                                    gAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcg                                                                                                                 >  1:724829/1‑151 (MQ=255)
                                     aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgc                                                                                                                >  2:177445/1‑151 (MQ=255)
                                         tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  2:56268/1‑151 (MQ=255)
                                                ggagTAAAGGCCACCCAGAATCCGGCGGCCGGCGTGGGGGGAAGGGGGGGGTGCAGCGGGGTGATAAGGGCTAAAATGACCCGGCGGTGATGGTGGGCCAAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                     <  1:384824/148‑1 (MQ=255)
                                                  cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  1:609525/1‑151 (MQ=255)
                                                   gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  2:282138/1‑150 (MQ=255)
                                                    tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 >  2:1227774/1‑151 (MQ=255)
                                                         ggCACCCTAAACCCAGCGGCCAGCGTCGGGGGATTGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGCCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg                                                                                              <  1:56268/149‑1 (MQ=255)
                                                          gCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                            >  1:735465/1‑150 (MQ=255)
                                                            aaCCTGACTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  2:1236093/1‑150 (MQ=255)
                                                              ccTGAATCCAGCGGCCAGCTTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  1:1117561/151‑1 (MQ=255)
                                                               ctgactcaggCGGCCGGCTTGGGGGGAAGGTGGGGGGGCAGGTGGGTAAAAAGGTCAAAAAGGACACCGCGGGGATGGTGTGCTAAAAATATAAGTTTGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  2:1111554/141‑1 (MQ=255)
                                                                  aacccAGGGCCCGGCGTGGGGGGAATGGGGGGGTGCAGCGGGGTAATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                     <  1:504852/146‑1 (MQ=255)
                                                                   acccggCGGCCTGGGTCTGGGAATTGGGGGGGTGGGGGTGGGTGATAAGGGCAACAGTACACGGGGAGTGAGGGGGGGCGTAAAATATATGTTTTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAAat                                                                                        <  2:535830/140‑1 (MQ=255)
                                                                     ccggCGGCCTGCTTTTGGGGAATGTGGGCGTGCGGGTGGTTAATAAGTGCAAAAATGACCCCGCAGTGGGGGTGGGCTTAAAATATAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:1165625/148‑1 (MQ=255)
                                                                        ggggCCGGCGTGGGGGGAATGGGGGGGTGCGGGTGGGTGATAAGGGCTAAAATAACCCCGAGGTGATGGGGGGCTTAAAATATAGTTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                               <  2:619559/147‑1 (MQ=255)
                                                                        gcgcCCTGCGTTGGGGGAAGGTGGGCGTGCAGTTGGTGAATAATGCCTACAATGCCCCCGCAGTGAGGGGGGCCCTAAAAAATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  1:1202297/148‑1 (MQ=255)
                                                                        gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  1:162905/1‑150 (MQ=255)
                                                                          ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           >  2:525100/1‑151 (MQ=255)
                                                                                    tGGGGATTGGGGGCTTGCAGTGGGGTGATAATGGCTACAATGCCCCCGCGGTGATGGTGGCCCTAAAAAATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                                   <  2:766262/149‑1 (MQ=255)
                                                                                    ggggaAATGGGGGCGTGCAGCTGGGTAAAAATGGCTAAAATACCCCCGCAGTGATGGTGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa                                                                  <  2:580305/149‑1 (MQ=255)
                                                                                            gggggggtgaggttGGGTAATATGGCCTAAAAGGACCCCGCAGTGAGGGTGGGCTAAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                             <  2:1014640/134‑1 (MQ=255)
                                                                                            gggggggggcggtggggTGATAATGGATACAATGACCCAGCAGTGGTGGTGGGCATAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                          <  2:1207248/136‑1 (MQ=255)
                                                                                            ggggggGTGCGGCGGGGTAATAATGGCTAAAATGACCCCGCGGTGGTGGTGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                          <  2:587664/148‑1 (MQ=255)
                                                                                               ggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                      >  1:67136/1‑151 (MQ=255)
                                                                                                 cGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                    >  1:915005/1‑151 (MQ=255)
                                                                                                  gTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                                   >  2:1099391/1‑151 (MQ=255)
                                                                                                       ggggggTAAAAATGGCTAAAATGACCCCGCAGTGAGGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                               <  1:156464/147‑1 (MQ=255)
                                                                                                       gCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                               >  1:547073/1‑150 (MQ=255)
                                                                                                           ggTAATAAGGGCTACAATGCCCCCGAAGTGGGGGGGGGACTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:870220/151‑1 (MQ=255)
                                                                                                                 aggggCTAAAATGACACGGAGGTGATGGTGGGCCAGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:818934/142‑1 (MQ=255)
                                                                                                                     gCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                 >  2:790531/1‑150 (MQ=255)
                                                                                                                        aCAATGACCCCGGTGAGAAGTTGGGCCTAAAATAAATGTTCTTATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                              <  1:789469/150‑1 (MQ=255)
                                                                                                                        aCAATGACACCGAAGTGAGGGGGGGCCTAAAATATAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCaa                                <  1:866326/148‑1 (MQ=255)
                                                                                                                         cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:175285/1‑151 (MQ=255)
                                                                                                                          attGACACCGCATTGATGGGGGGCCTAAAATATATGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                             <  1:122458/147‑1 (MQ=255)
                                                                                                                          aaTGACCCCGCAGTGAGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            <  1:582045/150‑1 (MQ=255)
                                                                                                                          aaTGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:610996/1‑150 (MQ=255)
                                                                                                                           acgccacCGCAGTTATGGCGCCCTTAAAATATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCa                           <  1:422728/146‑1 (MQ=255)
                                                                                                                            tGACACCGCAGTGATGGTGGGCCTGACATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGCGATCAACAGCAtt                         >  2:399880/1‑151 (MQ=255)
                                                                                                                               cacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  1:584980/1‑151 (MQ=255)
                                                                                                                                 caggatgTGATGGGCGGCCTAAAATATAAATATGATGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGGCATCAATTGCCGTGATCAACAGc                            <  2:1217818/138‑1 (MQ=255)
                                                                                                                                 ccGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  1:177433/147‑1 (MQ=255)
                                                                                                                                 ccGAAGTGAGGGGGGGCCAAAAATAAAAGTTCTTCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg                    <  2:465519/151‑1 (MQ=255)
                                                                                                                                   cGCAGTGGGGGTGGGCCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   <  1:219179/151‑1 (MQ=255)
                                                                                                                                      aGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                >  1:590093/1‑151 (MQ=255)
                                                                                                                                       gTGATGGTGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                <  1:10128/150‑1 (MQ=255)
                                                                                                                                        tGATGGGGGCCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                <  1:315030/149‑1 (MQ=255)
                                                                                                                                             gtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:187130/1‑151 (MQ=255)
                                                                                                                                              tggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:278187/1‑150 (MQ=255)
                                                                                                                                              tggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:1233179/1‑150 (MQ=255)
                                                                                                                                               ggtccTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa        <  1:571750/147‑1 (MQ=255)
                                                                                                                                                ggCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa       >  1:256277/1‑150 (MQ=255)
                                                                                                                                                 gcttGAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc     <  1:758443/148‑1 (MQ=255)
                                                                                                                                                 gCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc     <  1:237160/151‑1 (MQ=255)
                                                                                                                                                    tGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:1067852/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑1222382

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