Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,162 C→G 7.3% intergenic (+40/+61) PP_1661 → / ← PP_1662 dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,1620CG 7.3% 296.5 / 13.5 97intergenic (+40/+61)PP_1661/PP_1662dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (43/46);  new base G (2/5);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 101 total aligned reads displayed.

AGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGT  >  NC_002947/1859014‑1859310
                                                                                                                                                      |                                                                                                                                                    
aGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGAAGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCgca                                                                                                                                                     <  1:366844/150‑1 (MQ=255)
      gTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGACCTGTCATGGCAACCCCCCGATGCTCCGGCAATGCAGCAGCAttt                                                                                                                                               >  1:559794/1‑150 (MQ=255)
       ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGTTCCGGCAATGCCGCAGCAtttt                                                                                                                                              <  2:299995/150‑1 (MQ=255)
       ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCAtttt                                                                                                                                              >  2:1230387/1‑150 (MQ=255)
       ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCAtttt                                                                                                                                              >  2:747290/1‑150 (MQ=255)
        tGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTg                                                                                                                                             <  1:1536742/150‑1 (MQ=255)
         gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTGCGGCATTGCCGGAGCAtcgtgg                                                                                                                                            <  2:2174338/150‑6 (MQ=255)
         gCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgg                                                                                                                                            <  2:880726/150‑1 (MQ=255)
            gCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCATGGCATGGCAACCCCACGATGCTCCGGCAAGGCGGCGGCACTTTGGGGt                                                                                                                                         >  1:1060113/1‑150 (MQ=255)
             cTGGGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTTGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGtt                                                                                                                                        <  2:8667/150‑1 (MQ=255)
              tgtgCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTg                                                                                                                                       <  1:1490019/150‑1 (MQ=255)
                tgCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGcc                                                                                                                                     <  2:505721/150‑2 (MQ=37)
                 gCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTc                                                                                                                                    >  2:385821/1‑150 (MQ=255)
                   aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGGT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCAc                                                                                                                                  >  2:857139/1‑150 (MQ=255)
                    tCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATTGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACt                                                                                                                                 <  1:811453/150‑1 (MQ=255)
                      aaTGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGa                                                                                                                               <  2:1157917/150‑1 (MQ=255)
                        tGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGGAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGc                                                                                                                             <  2:119267/150‑1 (MQ=255)
                          acgacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGccc                                                                                                                           <  1:385821/150‑1 (MQ=255)
                           cgacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCa                                                                                                                          >  2:570396/1‑150 (MQ=255)
                             acCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAgg                                                                                                                        <  2:170597/150‑1 (MQ=255)
                              ccTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                                       <  1:897306/150‑1 (MQ=255)
                                tGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGa                                                                                                                     >  1:1876963/1‑150 (MQ=255)
                                     gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                                <  2:1876963/150‑1 (MQ=255)
                                     gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                                >  2:169978/1‑150 (MQ=255)
                                          ttACGGGGCTGGCCTGACCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTCACCGAGTCTACCTTCGGCCTG‑‑ATGACAACCCCACGATGATCCGGCACTGCCGAAAACCTCTGGTCTGTCGCACTGATCCCAGGCGAGTCAGGGGtt                                                                                                           <  1:2280894/150‑1 (MQ=255)
                                              aacacTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgcc                                                                                                       >  1:81334/5‑150 (MQ=255)
                                                  ccggCCTGATCCATATCCTGCGCCGCGTGCGGGTAGAGGTGCGTT‑‑TCGGCGCGCTGTAGGCCCTGTCATGGCAACCCCACGATTCTCTGGCAATGCCGCAGCCTTTTGCGATTGCTCACTGTCCCCAGGCGCGTCAGGGGTTCGCCGCaa                                                                                                   <  2:923511/148‑1 (MQ=255)
                                                        tGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                                             >  2:1618641/1‑150 (MQ=255)
                                                         gagcctaaaCCTCGGCATCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAg                                                                                            <  2:1333668/143‑1 (MQ=255)
                                                            cGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGAAACGCAGGGGTTCGCCGCAATGGGCAGGGt                                                                                         <  1:1230387/150‑1 (MQ=255)
                                                             gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTc                                                                                        >  2:2101375/1‑150 (MQ=255)
                                                                  ccTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACTATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgt                                                                                   >  2:616753/1‑150 (MQ=255)
                                                                  ccTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTCGCCCACTGAGCCCAGGCGAGTCAGGGGTTCTCCGCAATGGGCAGGGTCATtgt                                                                                   <  1:1696484/150‑1 (MQ=255)
                                                                  ccTGGCCCGCCTGTCGGTAGAGGTGCGTT‑‑TCGGTGCGGTGTAGGCCCTGTCAGGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgt                                                                                   <  1:1208333/150‑1 (MQ=255)
                                                                   cTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtg                                                                                  >  2:1917068/1‑150 (MQ=255)
                                                                    tGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTGATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgt                                                                                 <  2:2127472/150‑1 (MQ=255)
                                                                     gggCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                                >  2:1757872/1‑150 (MQ=255)
                                                                     gggCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                                <  1:1757872/150‑1 (MQ=255)
                                                                       gCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGgc                                                                              <  1:169978/150‑1 (MQ=255)
                                                                          gCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAtt                                                                                     >  2:1793640/1‑140 (MQ=255)
                                                                          gCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAtt                                                                                     <  1:1793640/140‑1 (MQ=255)
                                                                          gCCTGTCGGTAGAGGTGCGTG‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGgcgcg                                                                           >  2:2233938/1‑150 (MQ=255)
                                                                            cTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCt                                                                         >  2:1275025/1‑150 (MQ=255)
                                                                             tGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTg                                                                        <  1:1647434/150‑1 (MQ=255)
                                                                             tGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGAGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTg                                                                        <  2:81334/150‑1 (MQ=255)
                                                                               tCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                                   <  1:1082142/89‑1 (MQ=255)
                                                                               tCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                                   >  2:1082142/1‑89 (MQ=255)
                                                                                 ggTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                                    >  1:2201215/1‑150 (MQ=255)
                                                                                 ggTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                                    >  2:1325364/1‑150 (MQ=255)
                                                                                  gTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGtttt                                                                   >  2:2088808/1‑150 (MQ=255)
                                                                                        gTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGt                                                             >  2:1361079/1‑150 (MQ=255)
                                                                                        gTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGt                                                             >  2:1259843/1‑150 (MQ=255)
                                                                                            gTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTTACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTCCCAGGTCTTc                                                         >  2:1718321/1‑150 (MQ=255)
                                                                                              t‑‑tCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTTACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTCCCAGGTCTTCaa                                                       <  1:1718321/150‑1 (MQ=255)
                                                                                                 tCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAg                                                      <  1:316912/150‑1 (MQ=255)
                                                                                                      gCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCaa                                                 >  2:936028/1‑150 (MQ=255)
                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAg                                                                                            >  2:962321/1‑104 (MQ=255)
                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAg                                                                                            <  1:962321/104‑1 (MQ=255)
                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGc                                                     <  1:682860/143‑1 (MQ=255)
                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGc                                                     >  2:682860/1‑143 (MQ=255)
                                                                                                            tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                                            >  1:1013281/1‑149 (MQ=255)
                                                                                                            tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                                            <  2:1013281/149‑1 (MQ=255)
                                                                                                             aGGCCCTGTCATGGCAACCCGACTATGCTCCGGCAATGCCGCAGCATTTTGGGGTGGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTa                                          <  2:1060113/150‑1 (MQ=255)
                                                                                                               gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                                <  1:150342/110‑1 (MQ=255)
                                                                                                               gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                                >  2:150342/1‑110 (MQ=255)
                                                                                                                cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                       <  2:2201215/150‑1 (MQ=255)
                                                                                                                cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                       <  1:2233938/150‑1 (MQ=255)
                                                                                                                 ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                                      <  1:1325364/150‑1 (MQ=255)
                                                                                                                   tcagtgagCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCCCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtg                                    >  2:1886164/8‑150 (MQ=25)
                                                                                                                         ggCAACCCCACTATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                              <  1:616753/150‑1 (MQ=255)
                                                                                                                          gCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCg                             >  2:1245205/1‑150 (MQ=255)
                                                                                                                           cAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgt                            <  1:1618641/150‑1 (MQ=255)
                                                                                                                             aCCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtgt                          >  2:1461126/1‑150 (MQ=255)
                                                                                                                               cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                                      >  1:2003473/1‑136 (MQ=255)
                                                                                                                               cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                                      <  2:2003473/136‑1 (MQ=255)
                                                                                                                               cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                        >  2:1250553/1‑150 (MQ=255)
                                                                                                                               cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                        <  1:400459/150‑1 (MQ=255)
                                                                                                                                ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg                       <  1:570396/150‑1 (MQ=255)
                                                                                                                                 caaaaTGCTGCGGCATTGCCGGAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                           <  2:2104617/77‑1 (MQ=37)
                                                                                                                                 caaaaTGCTGCGGCATTGCCGGAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                           >  1:2104617/5‑81 (MQ=25)
                                                                                                                                 cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGAGGTATTGGTGTGCTTCGTGTTCGt                      <  2:1677191/150‑1 (MQ=255)
                                                                                                                                    gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                                <  2:1997225/105‑1 (MQ=255)
                                                                                                                                    gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                                <  1:1187223/105‑1 (MQ=255)
                                                                                                                                    gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                                >  2:1187223/1‑105 (MQ=255)
                                                                                                                                    gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                                >  1:1997225/1‑105 (MQ=255)
                                                                                                                                          ccGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACa             >  1:1405434/1‑150 (MQ=255)
                                                                                                                                            ggCATTGCCGGAGCATCGTGGGGTTGCCCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATg           <  1:1886164/150‑1 (MQ=255)
                                                                                                                                             gCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                              >  1:128910/1‑50 (MQ=255)
                                                                                                                                             gCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                              <  2:128910/50‑1 (MQ=255)
                                                                                                                                                aTGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCcca       >  2:1401463/1‑150 (MQ=255)
                                                                                                                                                 tGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGAccac      <  2:559794/150‑1 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                       >  2:769317/1‑114 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                       <  1:769317/114‑1 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAg   >  2:384494/1‑150 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAg   <  1:2088808/150‑1 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAg   >  1:1522563/1‑150 (MQ=255)
                                                                                                                                                    cGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAg   >  2:1219706/1‑150 (MQ=255)
                                                                                                                                                     gcagcaTTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGaa               >  1:1656190/1‑137 (MQ=255)
                                                                                                                                                     gcagcaTTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGaa               <  2:1656190/137‑1 (MQ=255)
                                                                                                                                                     gcagcaTTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGt  <  1:857139/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGT  >  NC_002947/1859014‑1859310

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