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

GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222085‑1222382
                                                                                                                                                     |                                                                                                                                                      
gtgAAAGTCGAAGGAGATGGGGCGAGTAATGTGGGGGAAATTAAACTAGGCGTAAATGGCAGTCTGAATCCAGCGGCCAGCGTCTGGTGAATGTGGGCGTGCAGCGGGGTGATAATCGCTACAATGACAC‑CGCAGTGGTGGTGGGCCTGaa                                                                                                                                                      <  2:726298/151‑1 (MQ=255)
gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGaa                                                                                                                                                      <  2:1246383/151‑1 (MQ=255)
   aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    >  1:465534/1‑150 (MQ=255)
        cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                              <  2:1236495/151‑1 (MQ=255)
         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGtt                                                                                                                                             >  1:728319/1‑151 (MQ=255)
            cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                                          >  1:2200/1‑151 (MQ=255)
            cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGAAAT‑TGGGCGGGCAGATTGGGGATTATGGGTAGAATAAACC‑CGCAGGGATGGGGGGCCTGGAattaaaggt                                                                                                                                             >  1:643957/1‑140 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                         <  1:1165169/151‑1 (MQ=255)
             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                         <  2:91079/151‑1 (MQ=255)
              aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                                        >  2:315023/1‑151 (MQ=255)
                      cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGt                                                                                                                               >  2:852569/1‑151 (MQ=255)
                          tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGc                                                                                                                           >  2:920775/1‑151 (MQ=255)
                          tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGc                                                                                                                           >  2:487281/1‑151 (MQ=255)
                              ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTaa                                                                                                                        >  1:720605/1‑150 (MQ=255)
                                  gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAAc                                                                                                                   >  2:272961/1‑151 (MQ=255)
                                  gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGAGGGCCTGAAATATAAGTTCTGAT‑CACGCCGTTAGCTaaaaa                                                                                                                     >  2:422715/1‑149 (MQ=255)
                                    gAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACcg                                                                                                                 >  1:724829/1‑151 (MQ=255)
                                     aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACcgc                                                                                                                >  2:177445/1‑151 (MQ=255)
                                         tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  2:56268/1‑151 (MQ=255)
                                                ggagTAAAGGCCACCCAGAATCCGGCGGCCGGCGTGGGGGGAAGGGGGGGGTGCAGCGGGGTGATAAGGGCTAAAATGACCC‑GGCGGTGATGGTGGGCCAAAAATATAATTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                     <  1:384824/148‑1 (MQ=255)
                                                  cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  1:609525/1‑151 (MQ=255)
                                                   gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  2:282138/1‑150 (MQ=255)
                                                    tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 >  2:1227774/1‑151 (MQ=255)
                                                         ggCACCCTAAACCCAGCGGCCAGCGTCGGGGGATTGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACCC‑CGCAGTGATGGTGGCCCTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg                                                                                              <  1:56268/149‑1 (MQ=255)
                                                          gCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                            >  1:735465/1‑150 (MQ=255)
                                                            aaCCTGACTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  2:1236093/1‑150 (MQ=255)
                                                              ccTGAATCCAGCGGCCAGCTTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  1:1117561/151‑1 (MQ=255)
                                                               ctgactcaggCGGCCGGCTTGGGGGGAAGGTGGGGGGGCAGGTGGGTAAAAAGGTCAAAAAGGACAC‑CGCGGGGATGGTGTGCTAAAAATATAAGTTTGGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                       <  2:1111554/141‑1 (MQ=255)
                                                                  aacccAGGGCCCGGCGTGGGGGGAATGGGGGGGTGCAGCGGGGTAATAATGGCTAAAATGACCC‑CGCAGTGATGGGGGGCCTAAAATATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                     <  1:504852/146‑1 (MQ=255)
                                                                   acccggCGGCCTGGGTCTGGGAATTGGGGGGGTGGGGGTGGGTGATAAGGGCAACAGT‑ACACGGGGAGTGAGGGGGGGCGTAAAATATATGTTTTGCT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAAat                                                                                        <  2:535830/140‑1 (MQ=255)
                                                                     ccggCGGCCTGCTTTTGGGGAATGTGGGCGTGCGGGTGGTTAATAAGTGCAAAAATGACCC‑CGCAGTGGGGGTGGGCTTAAAATATAAGTTCTGTT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  2:1165625/148‑1 (MQ=255)
                                                                        ggggCCGGCGTGGGGGGAATGGGGGGGTGCGGGTGGGTGATAAGGGCTAAAATAACCC‑CGAGGTGATGGGGGGCTTAAAATATAGTTTCTGTT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                               <  2:619559/147‑1 (MQ=255)
                                                                        gcgcCCTGCGTTGGGGGAAGGTGGGCGTGCAGTTGGTGAATAATGCCTACAATGCCCC‑CGCAGTGAGGGGGGCCCTAAAAAATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                <  1:1202297/148‑1 (MQ=255)
                                                                        gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  1:162905/1‑150 (MQ=255)
                                                                          ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           >  2:525100/1‑151 (MQ=255)
                                                                                    tGGGGATTGGGGGCTTGCAGTGGGGTGATAATGGCTACAATGCCCC‑CGCGGTGATGGTGGCCCTAAAAAATATTTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                                   <  2:766262/149‑1 (MQ=255)
                                                                                    ggggaAATGGGGGCGTGCAGCTGGGTAAAAATGGCTAAAATACCCC‑CGCAGTGATGGTGGGCCTAAAATAAAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa                                                                  <  2:580305/149‑1 (MQ=255)
                                                                                            gggggggtgaggttGGGTAATATGGCCTAAAAGGACCC‑CGCAGTGAGGGTGGGCTAAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                             <  2:1014640/134‑1 (MQ=255)
                                                                                            gggggggggcggtggggTGATAATGGATACAATGACCC‑AGCAGTGGTGGTGGGCATAAAATATAATTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                          <  2:1207248/136‑1 (MQ=255)
                                                                                            ggggggGTGCGGCGGGGTAATAATGGCTAAAATGACCC‑CGCGGTGGTGGTGGGCCTAAAATAAAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                          <  2:587664/148‑1 (MQ=255)
                                                                                               ggCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                      >  1:67136/1‑151 (MQ=255)
                                                                                                 cGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                    >  1:915005/1‑151 (MQ=255)
                                                                                                  gTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                                   >  2:1099391/1‑151 (MQ=255)
                                                                                                       ggggggTAAAAATGGCTAAAATGACCC‑CGCAGTGAGGGGGGGCTTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                               <  1:156464/147‑1 (MQ=255)
                                                                                                       gCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                               >  1:547073/1‑150 (MQ=255)
                                                                                                           ggTAATAAGGGCTACAATGCCCC‑CGAAGTGGGGGGGGGACTAAAATATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:870220/151‑1 (MQ=255)
                                                                                                                 aggggCTAAAATGACAC‑GGAGGTGATGGTGGGCCAGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:818934/142‑1 (MQ=255)
                                                                                                                     gCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                 >  2:790531/1‑150 (MQ=255)
                                                                                                                        aCAATGACCC‑CGGTGAGAAGTTGGGCCTAAAATAAATGTTCTTAT‑TCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                              <  1:789469/150‑1 (MQ=255)
                                                                                                                        aCAATGACAC‑CGAAGTGAGGGGGGGCCTAAAATATAAGTTCGGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCaa                                <  1:866326/148‑1 (MQ=255)
                                                                                                                         cAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:175285/1‑151 (MQ=255)
                                                                                                                          attGACAC‑CGCATTGATGGGGGGCCTAAAATATATGTTCGGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                             <  1:122458/147‑1 (MQ=255)
                                                                                                                          aaTGACCC‑CGCAGTGAGGGTGGGCCTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            <  1:582045/150‑1 (MQ=255)
                                                                                                                          aaTGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                            >  2:610996/1‑150 (MQ=255)
                                                                                                                           acgccac‑CGCAGTTATGGCGCCCTTAAAATATAATTTCTAAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCa                           <  1:422728/146‑1 (MQ=255)
                                                                                                                            tGACAC‑CGCAGTGATGGTGGGCCTGACATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGCGATCAACAGCAtt                         >  2:399880/1‑151 (MQ=255)
                                                                                                                               cac‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  1:584980/1‑151 (MQ=255)
                                                                                                                                 caggatgTGATGGGCGGCCTAAAATATAA‑ATATGATGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGGCATCAATTGCCGTGATCAACAGc                            <  2:1217818/138‑1 (MQ=255)
                                                                                                                                 c‑cGCAGTGATGGTGGGCCTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  1:177433/147‑1 (MQ=255)
                                                                                                                                 c‑cGAAGTGAGGGGGGGCCAAAAATAAAAGTTCTTCT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg                    <  2:465519/151‑1 (MQ=255)
                                                                                                                                   cGCAGTGGGGGTGGGCCTAAAATAAAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   <  1:219179/151‑1 (MQ=255)
                                                                                                                                      aGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                >  1:590093/1‑151 (MQ=255)
                                                                                                                                       gTGATGGTGGGCTTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                <  1:10128/150‑1 (MQ=255)
                                                                                                                                        tGATGGGGGCCCTAAAATATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                <  1:315030/149‑1 (MQ=255)
                                                                                                                                             gtggGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:187130/1‑151 (MQ=255)
                                                                                                                                              tggGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:278187/1‑150 (MQ=255)
                                                                                                                                              tggGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc         >  2:1233179/1‑150 (MQ=255)
                                                                                                                                               ggtccTAAAATATAATTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa        <  1:571750/147‑1 (MQ=255)
                                                                                                                                                ggCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa       >  1:256277/1‑150 (MQ=255)
                                                                                                                                                 gcttGAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc     <  1:758443/148‑1 (MQ=255)
                                                                                                                                                 gCCTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc     <  1:237160/151‑1 (MQ=255)
                                                                                                                                                    tGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:1067852/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACAC‑CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  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: ‑