Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,156 A→T 7.8% intergenic (+34/+67) 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,1560AT 7.8% 265.9 / 13.4 91intergenic (+34/+67)PP_1661/PP_1662dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base T (2/5);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

CAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCA  >  NC_002947/1859009‑1859305
                                                                                                                                                     |                                                                                                                                                     
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: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:1259843/1‑150 (MQ=255)
                                                                                             gTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGt                                                        >  2:1361079/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                                                                           >  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                   >  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:1187223/1‑105 (MQ=255)
                                                                                                                                         gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                           <  1:1187223/105‑1 (MQ=255)
                                                                                                                                         gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                           >  1:1997225/1‑105 (MQ=255)
                                                                                                                                         gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                           <  2:1997225/105‑1 (MQ=255)
                                                                                                                                               ccGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACa        >  1:1405434/1‑150 (MQ=255)
                                                                                                                                                 ggCATTGCCGGAGCATCGTGGGGTTGCCCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATg      <  1:1886164/150‑1 (MQ=255)
                                                                                                                                                  gCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                         <  2:128910/50‑1 (MQ=255)
                                                                                                                                                  gCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                         >  1:128910/1‑50 (MQ=255)
                                                                                                                                                     aTGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCcca  >  2:1401463/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTT‑‑TCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCA  >  NC_002947/1859009‑1859305

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: