Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,367,726 G→T L63I (CTT→ATT)  pspF ← DNA‑binding transcriptional dual regulator PspF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,7260GT98.0% 159.3 / ‑5.7 49L63I (CTT→ATT) pspFDNA‑binding transcriptional dual regulator PspF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (24/24);  minor base A (0/1);  total (24/25)
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.86e-01

CCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGT  >  NC_000913/1367480‑1367968
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ccAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGa                                                                                                                                                                                                                                                 >  1:244715/1‑250 (MQ=255)
         caTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGc                                                                                                                                                                                                                                       >  2:50529/1‑251 (MQ=255)
                  gCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGcc                                                                                                                                                                                                                                <  1:150541/249‑1 (MQ=255)
                       ggTTGGCGGCCGCCAACGCGCTCCAGTTCACCGTACCCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAAGGAAATAAACGGCCCTTGCCAACGgga                                                                                                                                                                                                                         <  2:244715/251‑1 (MQ=255)
                            gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGgaga                                                                                                                                                                                                                    <  1:109388/251‑1 (MQ=255)
                                 cgccAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATaa                                                                                                                                                                                                               <  1:178693/251‑1 (MQ=255)
                                      aCGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGc                                                                                                                                                                                                            >  1:457636/1‑249 (MQ=255)
                                            tCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATAATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGACGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGATAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAggcggcg                                                                                                                                                                                                    >  2:86164/1‑250 (MQ=255)
                                                              aacccgcgCAATAATGTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCAGCAAGAAATACCGTACCTCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTTCGCACCGGTAAACGCCCCCGCTGCGTGACCAAACAGTTCGGAAGCCATCATATTTTCATTGAACGCCGCGCAGTGAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCTGGCGTCTGGCAATCAGCTCCTTAcc                                                                                                                                                                                  <  2:566843/247‑1 (MQ=255)
                                                                cACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCgg                                                                                                                                                                                >  2:456490/1‑251 (MQ=255)
                                                                             ttCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGcc                                                                                                                                                                   >  1:544313/1‑251 (MQ=255)
                                                                             ttCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGcc                                                                                                                                                                   >  1:396031/1‑251 (MQ=255)
                                                                              tctcCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACAGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGCGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCg                                                                                                                                                                  >  2:26351/1‑251 (MQ=255)
                                                                                tcCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCgat                                                                                                                                                                >  1:2244/1‑251 (MQ=255)
                                                                                    gctccATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAGCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGGGTCCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACTCCGCGCAGTTAATGGAAAGAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGTCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGatgatg                                                                                                                                                            <  1:72214/248‑1 (MQ=255)
                                                                                    gCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGatgatg                                                                                                                                                            >  1:446450/1‑251 (MQ=255)
                                                                                           catcGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGTCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGAGGATGAGCACCg                                                                                                                                                     <  2:341332/251‑1 (MQ=255)
                                                                                              cGGTGCCGTAGCGAGTTCAGCAAGAAATAGCGTACCGCCGTGGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTTCCGGTGCCGCGTTCGCCGATGATGAGCACCGGtt                                                                                                                                                  <  2:419344/251‑1 (MQ=255)
                                                                                                   ccGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGGTTGTcc                                                                                                                                             <  2:458092/251‑1 (MQ=255)
                                                                                                     gTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTcc                                                                                                                                             >  1:8984/1‑249 (MQ=255)
                                                                                                       aGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCGTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTGACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCg                                                                                                                                         <  1:38043/251‑1 (MQ=255)
                                                                                                        gCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCgg                                                                                                                                        <  1:367483/251‑1 (MQ=255)
                                                                                                          gAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTg                                                                                                                                      >  1:293476/1‑251 (MQ=255)
                                                                                                                atcaAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAAGCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGAt                                                                                                                                <  2:354868/251‑1 (MQ=255)
                                                                                                                  caAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGc                                                                                                                              >  1:504181/1‑251 (MQ=255)
                                                                                                                   aaGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCg                                                                                                                             >  1:453982/1‑251 (MQ=255)
                                                                                                                       aagaGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAAc                                                                                                                         <  2:396031/248‑1 (MQ=255)
                                                                                                                          atcttACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCTCCGGTAACCGCCCCCGCGTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCAGTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTTGTTACCGGTGCCTCGGTCGCCGATGATGAGCACCGGTTTGACCAGCGGTGCTAGATGCGAAACCt                                                                                                                       <  2:171947/246‑1 (MQ=255)
                                                                                                                                              ttCAAATCTCCATGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCGGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGACGGCTGGCAATCAGCTCTTTAGCGGAGCGGCGTGCGCCGATGCTGAGCAACGGTTTGGCCAGGGGTGCGAGATGAGAAACCTATTcaggcccttcgcgaacgc                                                                                                  >  1:312057/1‑234 (MQ=255)
                                                                                                                                               tCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCt                                                                                                 >  2:168870/1‑251 (MQ=255)
                                                                                                                                                cAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTg                                                                                                <  1:252144/251‑1 (MQ=255)
                                                                                                                                                 aaaTCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGt                                                                                               <  2:423910/251‑1 (MQ=255)
                                                                                                                                                               gACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGt                                                                                 >  1:101076/1‑251 (MQ=255)
                                                                                                                                                               gACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCGGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGt                                                                                 <  1:359089/251‑1 (MQ=255)
                                                                                                                                                                  gTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTaa                                                                               >  1:64671/1‑250 (MQ=255)
                                                                                                                                                                         cgcACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCtt                                                                       <  1:2733/251‑1 (MQ=255)
                                                                                                                                                                          gcACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCttt                                                                      >  1:154239/1‑251 (MQ=255)
                                                                                                                                                                                  aaaCGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCg                                                                                            <  1:80273/221‑1 (MQ=255)
                                                                                                                                                                                  aaaCGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCg                                                                                            >  2:80273/1‑221 (MQ=255)
                                                                                                                                                                                       cccccGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATGTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATTATGAGCACCGGTTTGTCCAGCTGTGCGAGATGCGAAGCCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCatg                                                         <  1:28095/251‑1 (MQ=255)
                                                                                                                                                                                            gCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGaa                                                    <  1:168870/251‑1 (MQ=255)
                                                                                                                                                                                                     ccAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGCGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCAc                                           >  2:458898/1‑251 (MQ=255)
                                                                                                                                                                                                                        gcagATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCt                         >  1:119365/1‑250 (MQ=255)
                                                                                                                                                                                                                           gATTTTCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACt                     <  1:50529/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAACTTCGCCACTTGTTAGTGTAATTCGCTAACTCATc                 >  1:147507/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttCATTTAACGCCGCGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCAt                  <  2:154239/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cgccgcGCAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGATCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGCTGCGAAACCTGTTCCCGCACTTCGAGCAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGt        >  1:45176/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cgcAGTTAATGGAAATAAACGGCCCTTGCCAGCGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTTGTGTAATTCGCTAACTCATCCTGGCATGTTGCTg   <  2:272764/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gcAGTTAATGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCCGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTGCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGTCATGGTGCTgg  <  2:110645/251‑2 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
CCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGT  >  NC_000913/1367480‑1367968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑