Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,966 C→G 8.3% intergenic (+314/+165) PP_4318 → / ← PP_4319 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9660CG 8.3% 243.6 / 9.7 84intergenic (+314/+165)PP_4318/PP_4319transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base G (4/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCA  >  NC_002947/4909818‑4910115
                                                                                                                                                    |                                                                                                                                                      
ccTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCTGGGTTCCCTCGATAGGGTCGTTCTCCCGGGAGTACGCGCTGCAGGCCccatcc                                                                                                                                                       >  2:1353659/1‑150 (MQ=11)
ccTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcc                                                                                                                                                       >  1:647194/1‑150 (MQ=38)
ccTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcc                                                                                                                                                       <  1:716584/150‑1 (MQ=38)
 cTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGGGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGCTTGCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcca                                                                                                                                                      >  1:1895846/1‑150 (MQ=255)
      cGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTGCGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgggg                                                                                                                                                 <  2:1618070/150‑1 (MQ=18)
       gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGGGTACGCGCTGCAGGCCCCATCCATGGGGc                                                                                                                                                >  1:488959/1‑150 (MQ=17)
       gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGc                                                                                                                                                <  2:75685/150‑1 (MQ=32)
       gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGc                                                                                                                                                >  1:996137/1‑150 (MQ=32)
                 cggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                                     <  2:1941960/150‑1 (MQ=32)
                   ggTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                   <  2:1911360/150‑1 (MQ=32)
                           tttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATc                                                                                                                           <  2:647194/150‑1 (MQ=32)
                               tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATc                                                                                                                           >  1:1190097/1‑146 (MQ=32)
                               tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATc                                                                                                                           <  2:1190097/146‑1 (MQ=32)
                                 aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGcc                                                                                                                     >  2:87191/1‑150 (MQ=32)
                                  aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTTCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGGACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGAATCCATGccc                                                                                                                    >  1:1120992/1‑150 (MQ=12)
                                  aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                    >  1:1467003/1‑150 (MQ=32)
                                     gcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                 <  1:2258347/150‑1 (MQ=32)
                                          aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATGGGGCCtgcagc                                                                                                                                         >  1:58400/1‑116 (MQ=11)
                                          aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATGGGGCCtgcagc                                                                                                                                         <  2:58400/121‑6 (MQ=11)
                                                tAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGGGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCg                                                                                                      <  2:488959/150‑1 (MQ=25)
                                                tAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCg                                                                                                      <  2:1406795/150‑1 (MQ=34)
                                                 aaCCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGc                                                                                                     >  1:935169/1‑150 (MQ=35)
                                                   ccAAAAACCGCTGACCCCTAGCATCCGGCCCCTGCGCGACGCTCCGGGGTGCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                   <  2:450865/150‑1 (MQ=11)
                                                    cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatccat                                                                                                                                                     <  2:1521993/100‑1 (MQ=34)
                                                    cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatccat                                                                                                                                                     >  1:1521993/1‑100 (MQ=34)
                                                         aaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                                     >  1:123631/1‑110 (MQ=33)
                                                         aaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                                     <  2:123631/110‑1 (MQ=33)
                                                          aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTgc                                                                                            <  1:2018583/150‑1 (MQ=38)
                                                              tGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                        >  1:1362389/1‑150 (MQ=38)
                                                                  tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCtg                                                                                                                                             >  1:1264553/1‑93 (MQ=18)
                                                                  tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCtg                                                                                                                                             <  2:1264553/94‑2 (MQ=21)
                                                                     tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                >  2:329311/1‑103 (MQ=11)
                                                                     tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                <  1:329311/103‑1 (MQ=11)
                                                                            cGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTc                                                                          <  2:996137/150‑1 (MQ=39)
                                                                                 ccTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAAt                                                                     >  1:74065/1‑150 (MQ=38)
                                                                                     cgcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                 <  2:1467003/150‑1 (MQ=39)
                                                                                      gcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCg                                                                >  1:2131854/1‑150 (MQ=39)
                                                                                        tACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTc                                                              >  1:88416/1‑150 (MQ=38)
                                                                                        tACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTc                                                              >  2:797991/1‑150 (MQ=39)
                                                                                          cGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                            <  1:1049398/150‑1 (MQ=39)
                                                                                           gCTCCTGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCACTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGcc                                                           <  1:1353659/150‑1 (MQ=255)
                                                                                             tCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTg                                                         >  2:834656/1‑150 (MQ=39)
                                                                                               cGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                       <  1:2322025/150‑1 (MQ=39)
                                                                                                  ggTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTa                                                    <  2:2131854/150‑1 (MQ=39)
                                                                                                   gTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTAc                                                   >  2:377219/1‑150 (MQ=38)
                                                                                                    ttCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTAcc                                                  >  1:197930/1‑150 (MQ=40)
                                                                                                        cTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                     >  1:565175/1‑111 (MQ=38)
                                                                                                        cTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                     <  2:565175/111‑1 (MQ=38)
                                                                                                         tCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcgcg                                             >  2:1719541/1‑150 (MQ=40)
                                                                                                             tACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCaa                                                                      <  2:2126121/121‑1 (MQ=40)
                                                                                                             tACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCaa                                                                      >  1:2126121/1‑121 (MQ=40)
                                                                                                             tACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTa                                         >  1:1670244/1‑150 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      >  1:1864200/1‑150 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      <  1:1719541/150‑1 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      <  2:1864200/150‑1 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCATTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                      <  1:991831/150‑1 (MQ=255)
                                                                                                                  tCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTGCCGCGCGCTTACAAGc                                    >  1:1398223/1‑150 (MQ=255)
                                                                                                                  tCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGc                                    <  1:957414/150‑1 (MQ=255)
                                                                                                                  tCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGc                                    >  1:1217953/1‑150 (MQ=255)
                                                                                                                    ttGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaa                                  >  1:18731/1‑150 (MQ=255)
                                                                                                                     tGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGTCTGAACTACCGCGCGCTTACAAGCaaa                                 >  1:784905/1‑150 (MQ=255)
                                                                                                                     tGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACCGACAGGCAGGCCCGTCAAGTCCCTTCTCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaca                                 >  1:2037240/1‑148 (MQ=14)
                                                                                                                       cTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCATGAACTACCGCGCGCTTACAAGCAAaag                               <  1:2176827/150‑1 (MQ=255)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTAcc                                                  <  2:1891156/128‑1 (MQ=39)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTAcc                                                  >  1:1891156/1‑128 (MQ=38)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcag                            >  2:2037886/1‑150 (MQ=255)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcag                            <  1:797991/150‑1 (MQ=255)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATAAATCCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcag                            <  1:614461/150‑1 (MQ=255)
                                                                                                                            gggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                          >  2:1230767/1‑150 (MQ=255)
                                                                                                                            gggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                          <  2:1217953/150‑1 (MQ=255)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaag                         <  1:377219/150‑1 (MQ=255)
                                                                                                                             ccgtcaagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCGTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaag                         >  1:2010236/8‑150 (MQ=17)
                                                                                                                               aGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATCCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagca                       >  2:543672/1‑150 (MQ=255)
                                                                                                                                  aCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaa                    <  1:2037886/150‑1 (MQ=255)
                                                                                                                                    gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc                  >  1:496267/1‑150 (MQ=255)
                                                                                                                                        tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                        <  2:2020841/76‑1 (MQ=35)
                                                                                                                                        tGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                        >  1:2020841/1‑76 (MQ=34)
                                                                                                                                         gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagc                        <  1:533843/139‑1 (MQ=255)
                                                                                                                                         gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagc                        >  2:533843/1‑139 (MQ=255)
                                                                                                                                          cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagc            <  2:2037240/150‑1 (MQ=255)
                                                                                                                                            ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag          <  1:1878209/150‑1 (MQ=255)
                                                                                                                                             gCCCCATGGATGGGGCCTGAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgc                                                <  2:1830211/111‑1 (MQ=255)
                                                                                                                                             gCCCCATGGATGGGGCCTGAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgc                                                >  1:1830211/1‑111 (MQ=255)
                                                                                                                                             gCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAg                                     <  1:119249/122‑1 (MQ=255)
                                                                                                                                             gCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAg                                     >  2:119249/1‑122 (MQ=255)
                                                                                                                                             gCCCCATCCATGGGTCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagc         >  2:198394/1‑150 (MQ=255)
                                                                                                                                             gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcaga                           >  2:1418848/1‑132 (MQ=255)
                                                                                                                                             gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcaga                           <  1:1418848/132‑1 (MQ=255)
                                                                                                                                             gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagc         >  2:1232768/1‑150 (MQ=255)
                                                                                                                                               cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcag       <  2:74065/150‑1 (MQ=255)
                                                                                                                                                ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcaga      >  1:867140/1‑150 (MQ=255)
                                                                                                                                                ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcaga      <  2:935169/150‑1 (MQ=255)
                                                                                                                                                  atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag          <  2:750091/144‑1 (MQ=255)
                                                                                                                                                  atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag          >  1:750091/1‑144 (MQ=255)
                                                                                                                                                  atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaag    <  2:898964/150‑1 (MQ=255)
                                                                                                                                                    ccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGACTCCGAAAGGCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAgcagaagcagcagaagca  <  2:1355462/150‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                      
CCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCA  >  NC_002947/4909818‑4910115

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