Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,688 A→T 9.1% intergenic (+31/+21) fadR → / ← ycgB fatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,6880AT 9.1% 137.9 / 2.9 53intergenic (+31/+21)fadR/ycgBfatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base T (2/4);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGCGTTTGCTTCCGGCAACCCGATTTACGGTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCA  >  NC_000913/1235389‑1235988
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
gCCTGGCGGTTGGTTCCGTCAACCCGATTTATGGTTTTATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAAc                                                                                                                                                                                                                                                                                                             <  1:626033/301‑1 (MQ=255)
 ccTGGCGTTTGCTTCCGGCAACCCGATTTACGGTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAAAGCAAACCCCTCAAAca                                                                                                                                                                                                                                                                                                            >  2:228775/1‑300 (MQ=255)
    ggggTTTGCTTCCGGCAACCCGATTTACGGTCTTATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAgg                                                                                                                                                                                                                                                                                                         <  1:31566/298‑1 (MQ=255)
            tcagacgtttgttcttccgatctCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgct                                                                                                                                                                                                                                                                                                 <  2:169557/279‑4 (MQ=255)
              tagacgtgtgctcttccgatcttGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgctct                                                                                                                                                                                                                                                                                               <  2:461234/279‑6 (MQ=255)
               ccGGCAACCCGTTTTTCGGTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAAGCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgtt                                                                                                                                                                                                                                                                                              <  1:611127/301‑1 (MQ=255)
                cGGCAACCCTATTTTCGGTCTGATTCTTAACGGGATGAAAGGGCTGTTTACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgttg                                                                                                                                                                                                                                                                                             <  2:353252/301‑1 (MQ=255)
                cGGCAACCCGATTTACGGTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgttg                                                                                                                                                                                                                                                                                             <  2:358948/301‑1 (MQ=255)
                                ggTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTccc                                                                                                                                                                                                                                                                             >  2:599456/1‑301 (MQ=255)
                                ggTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAAAGAGGGGTTTTTTGTTGTTTTTACAGATTTccc                                                                                                                                                                                                                                                                             >  1:496647/1‑301 (MQ=255)
                                ggTCGGATTTTTGACGTGATGATAGGTCTGTGTACGCGTGTTGGTCGTGACTTTTTCTCCAATCCGTAATCTCGCAGTCTGGCGCTTGGCTTCTACCACTAGCCGTCGGCGTTGTGCGTTGATGGCGCGCACGATCATGTGTACAGAACAGTGCGTCGCGAGTGGCATGAGAGTGGCGAGATGTGTCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATGATCCTTTTCGTTTGAAGAGCAAACCCCTCAAACGAGGTGTTTTTTGTTGTTTTTACAGATTGccc                                                                                                                                                                                                                                                                             <  1:178541/301‑1 (MQ=255)
                                   cTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGAAAACCCCTCGTTTGAGGGGtttgctagatcggaagagcggcgggtagg                                                                                                                                                                                                                                                                          >  1:169360/1‑275 (MQ=255)
                                    tGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCACGTTTGAGGGGtttgctatagatcggaagagcgtcgtgcag                                                                                                                                                                                                                                                                         >  1:460670/1‑274 (MQ=255)
                                                       atggCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGTGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAg                                                                                                                                                                                                                                                      <  2:323580/299‑1 (MQ=255)
                                                           cTGTATACGTGTATTGGGCGTCACTATTTCGCCAAGCCGGAAGCGCGCGGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGGTCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAAc                                                                                                                                                                                                                                                  <  1:456548/301‑1 (MQ=255)
                                                               ataCGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCAATTCTTGGCGGGCAACGTTCCAGCAACTCGa                                                                                                                                                                                                                                              >  2:329675/1‑301 (MQ=255)
                                                                          ggTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCt                                                                                                                                                                                                                                   >  2:516351/1‑301 (MQ=255)
                                                                                                  gAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGCTTTCCCATTCTTGGCGGTCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacc                                                                                                                                                                                                           >  1:235459/1‑300 (MQ=255)
                                                                                                    ggcgcgcAGTCTGGCTCGGGGGTTCTACCCCAAACGTGCGGGGTTGTGCAGTGAAGGCGCGCACGAGCAGGTGTACGAAACAGGGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCTATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGTGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatc                                                                                                                                                                                                         <  1:458543/300‑1 (MQ=255)
                                                                                                       gcgcAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCaaa                                                                                                                                                                                                      >  2:275724/1‑301 (MQ=255)
                                                                                                              ggcgggtcttccgatttccACAAACGGTCGGCGTTGTGCAGTGAAGGCGCGCACGGTCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCcac                                                                                                                                                                                               <  1:20525/284‑1 (MQ=255)
                                                                                                               ttgcgcTGGGCTTCGACCACAAACGGTCGGCGTTTTGCAATGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTTATTTAGCCATTCAGGGGCGATAATCACTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCcaca                                                                                                                                                                                              <  2:614628/299‑1 (MQ=255)
                                                                                                                 gcgcGGGGGGTCTACCCCACATTGTCGGGGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAgg                                                                                                                                                                                            <  1:598695/301‑1 (MQ=255)
                                                                                                                  ggcggggCTTCTACCGCAGTGTGGTGGCGGTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGGGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGc                                                                                                                                                                                           <  1:97325/297‑1 (MQ=255)
                                                                                                                  cgcTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCAAGCATCACATCAAATCCCCACCGGc                                                                                                                                                                                           >  1:446092/1‑301 (MQ=255)
                                                                                                                   gcTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCg                                                                                                                                                                                          <  2:489451/301‑1 (MQ=255)
                                                                                                                             tACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGc                                                                                                                                                                                <  2:618545/301‑1 (MQ=255)
                                                                                                                               ccACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAgatcggaagagcacaa                                                                                                                                                                              >  2:20576/1‑286 (MQ=255)
                                                                                                                                cacaAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATCATCCCTTCCGTTCAAAGAGCACACCCCTCCACCGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGTCGGGCAACGTTCCAGCCACTCGATGCTGCCGCCTTACTTTTGCTGTTTGCGACTCACATCACATCACCACAAGCCATGCAGATGCttt                                                                                                                                                                             >  1:337874/1‑300 (MQ=255)
                                                                                                                                cacaAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTc                                                                                                                                                                             <  2:235725/301‑1 (MQ=255)
                                                                                                                                              gTTTTGCATTGAAGGGGCGTAGTATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcg                                                                                                                                                               <  1:275421/301‑1 (MQ=255)
                                                                                                                                                      gTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGAt                                                                                                                                                       <  2:446656/301‑1 (MQ=255)
                                                                                                                                                                          ttgtACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGTTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTAtgt                                                                                                                                   <  2:497294/300‑1 (MQ=255)
                                                                                                                                                                                    cAGTGCGTCGCTACGGGCATTAGATTTGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAAc                                                                                                                         <  2:435601/301‑1 (MQ=255)
                                                                                                                                                                                      ttccctacacgacgctcgtctgatctcGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACGCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCGTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGc                                                                                                                       <  1:18886/275‑1 (MQ=255)
                                                                                                                                                                                        gCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAg                                                                                                                     >  2:599461/1‑301 (MQ=255)
                                                                                                                                                                                                   ggCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAAc                                                                                                          >  2:30786/1‑301 (MQ=255)
                                                                                                                                                                                                                cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTCAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAgatcggaagagcacacgtctgaact                                                                                             >  2:18937/1‑277 (MQ=255)
                                                                                                                                                                                                                cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTAGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCCCCGCGCac                                                                                             >  1:635801/1‑300 (MQ=255)
                                                                                                                                                                                                                              ggATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGa                                                                               >  2:259409/1‑301 (MQ=255)
                                                                                                                                                                                                                                            ttGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGc                                                                 <  2:636678/300‑1 (MQ=255)
                                                                                                                                                                                                                                                 gggTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGACTATAACGCAGCGTCAGCGAACGGTCGCCGCGCCAATCCACGTTCCAGATCTGAATATTCGTCTCCAg                                                            >  1:520673/1‑301 (MQ=255)
                                                                                                                                                                                                                                                       tttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACt                                                      >  2:506964/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         tAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTa                                                    >  2:147067/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            ccATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGt                                                 <  2:802/301‑1 (MQ=255)
                                                                                                                                                                                                                                                               gtcgggggCGATATTCCCTTCCGTTTAAAGTGCAAGTCCCTCAAACGAGTGGTTTTTTGTTTTTTTTACAGCTTTCCCATTCTTGGCGGGCAACGTGCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTtat                                              <  1:30735/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAAACCACGTTCCAGATCTGAATATTCGGCCCCAGAATAATTAAGTtaacatt                                          >  1:200149/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                         ttatttcttttcGTTTAGTGAGCAAACCCCTCAACGTAGGGGTTTTTTGTTTTTTTTAAAGATTGCGCAGTCTTGGCGGGCAACGTTACAGCAACTCGAGGCGGCCGGCTGCGTTTTGCTGTTCGGGCATAACATCAAATCCCCACAGGAGATGCACATGCTGCAGGACTTCTTGGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGcgacga                                    <  1:228532/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ttCCGTTTTAGCGCAAACCCCTCAATCTGTGGGTTTTTTGTTGTTTTTACAGTTTTCCCATTCTTGGCGGGCAACGGTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGCGACGATAACCGGtt                           <  1:284264/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     ggtgagaggggagatcccTGAAAGGGGGGGTTTGTTGTTGTTTTTACAGATTTCCCAGGTGTGGCGGGAAACGTGCCAGCAACGCGAGGCTGGCGTCTTCGTGTTGCGGTTCGAGCATCACAGCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCGGTGCGCCCCCGATCCAGCGGTGCGCGATTATGGGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCGCCAGATTACTTAAGTTATATTGCGACGATAACCGGTTACg                        <  1:515679/286‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTAAATTGCGACGATAACCGGTTagatcggaa                  >  2:284567/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                                                  cccccTCAATCGATGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATaa           <  2:200392/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      cTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGACCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTAAGCGAACGGTCGCCCCGCAAATCCACGTTCCAGATCTGAATCTTCGGCTCCAGATTAATTAAGTTATATTGAGACGATAACCGGTTAAGGAACTCCCGATaaacat       >  1:498312/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                                                           aCGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCa  <  2:498959/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
GCCTGGCGTTTGCTTCCGGCAACCCGATTTACGGTCTGATTCTTAACGGGATGAAAGGGCTGTATACGCGTATTGGTCGTCACTATTTCGCCAATCCGGAAGCGCGCAGTCTGGCGCTGGGCTTCTACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGATTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCA  >  NC_000913/1235389‑1235988

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: