Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,150 C→G 6.8% intergenic (+28/+73) 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,1500CG 6.8% 276.0 / 10.6 89intergenic (+28/+73)PP_1661/PP_1662dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/45);  new base G (2/4);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACA  >  NC_002947/1859001‑1859299
                                                                                                                                                       |                                                                                                                                                     
cgcgTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGGCCTGTCATGGCAACCCCACGATGCTc                                                                                                                                                       >  1:1157917/1‑150 (MQ=255)
 gcgTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTcc                                                                                                                                                      >  2:2280894/1‑150 (MQ=255)
        cAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATg                                                                                                                                               >  2:2027398/1‑150 (MQ=255)
             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:747290/1‑150 (MQ=255)
                    ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCAtttt                                                                                                                                   >  2:1230387/1‑150 (MQ=255)
                     tGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTGCGGCATTGCCGGAGCAtcgtg                                                                                                                                  <  2:993589/150‑5 (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:169978/1‑150 (MQ=255)
                                                  gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                     <  2:1876963/150‑1 (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)
                                                                      gAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCGACCCCACGATGCTCCgg                                                                                                                                                    >  1:712808/1‑83 (MQ=255)
                                                                      gAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCgg                                                                                                                                                    <  2:712808/83‑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                                                                     <  1:1757872/150‑1 (MQ=255)
                                                                                  gggCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                     >  2:1757872/1‑150 (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                                                                                                                        >  2:1082142/1‑89 (MQ=255)
                                                                                            tCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                        <  1:1082142/89‑1 (MQ=255)
                                                                                              ggTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                         >  2:1325364/1‑150 (MQ=255)
                                                                                              ggTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                         >  1:2201215/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                                          >  2:682860/1‑143 (MQ=255)
                                                                                                                      gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGc                                          <  1:682860/143‑1 (MQ=255)
                                                                                                                         tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                                 >  1:1013281/1‑149 (MQ=255)
                                                                                                                         tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                                 <  2:1013281/149‑1 (MQ=255)
                                                                                                                          aGGCCCTGTCATGGCAACCCGACTATGCTCCGGCAATGCCGCAGCATTTTGGGGTGGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTa                               <  2:1060113/150‑1 (MQ=255)
                                                                                                                            gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                     >  2:150342/1‑110 (MQ=255)
                                                                                                                            gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                     <  1:150342/110‑1 (MQ=255)
                                                                                                                             cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                            <  1:2233938/150‑1 (MQ=255)
                                                                                                                             cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                            <  2:2201215/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                           <  2:2003473/136‑1 (MQ=255)
                                                                                                                                            cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                           >  1:2003473/1‑136 (MQ=255)
                                                                                                                                            cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc             <  1:400459/150‑1 (MQ=255)
                                                                                                                                            cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc             >  2:1250553/1‑150 (MQ=255)
                                                                                                                                             ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg            <  1:570396/150‑1 (MQ=255)
                                                                                                                                              caaaaTGCTGCGGCATTGCCGGAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                >  1:2104617/5‑81 (MQ=25)
                                                                                                                                              caaaaTGCTGCGGCATTGCCGGAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                <  2:2104617/77‑1 (MQ=37)
                                                                                                                                              cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGAGGTATTGGTGTGCTTCGTGTTCGt           <  2:1677191/150‑1 (MQ=255)
                                                                                                                                                 gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                     >  2:1187223/1‑105 (MQ=255)
                                                                                                                                                 gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                     <  2:1997225/105‑1 (MQ=255)
                                                                                                                                                 gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                     <  1:1187223/105‑1 (MQ=255)
                                                                                                                                                 gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                     >  1:1997225/1‑105 (MQ=255)
                                                                                                                                                       ccGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACa  >  1:1405434/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACA  >  NC_002947/1859001‑1859299

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: