Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,368,956 C→T 20.0% L98L (CTA→TTA)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,368,9560CT20.0% 144.3 / 5.6 70L98L (CTA→TTA) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base C (18/36);  new base T (12/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 2.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTG  >  NC_000913/4368721‑4369194
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
tttttttGAGGTTTCCTTTGTTTTGGTGGTATTGCTGCCCATCGGGCTTGTTAAATTCACGGCAGTTATCGGTATGTCACCGGTACGTAACCAGGGCTTTAAGAATTTCGTGCAGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTAc                                                                                                                                                                                                                                                   <  2:101647/238‑1 (MQ=255)
          gTTGCCCATGTATTGGGGGAATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAACAATGGAGGCGGGTGAAAACCCAGCGTCCGGGGAGAGTTAAAAGGGTTTCCTGTGCATTGTTGGTTTGTGTGTTTTTTTACCGTGCGTTTTTACTCTCACCCTCGGCCTTTTACCTTTTTTTgc                                                                                                                                                                                                                                          >  2:526732/1‑239 (MQ=255)
               ccATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGATTAAAAGGGTTTCGCTGATCATTGCTTGGTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTTCTCTTTTATTtccgcccg                                                                                                                                                                                                                                    >  1:38883/1‑237 (MQ=255)
                    ttttGGGGGTTTTGGTTGCCCTCTTGTTGGTTTTTTTTACGGTATTTATCGGTATGGCACTGGTACGGAACCAGGGGTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  2:274777/237‑1 (MQ=255)
                    ttttGGGGGTTTTGCTGCCCCGTGTGTTGGTTTTTTTCCAGTCATTTATCGGGATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  2:678271/237‑1 (MQ=255)
                     aTTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  1:647178/238‑1 (MQ=255)
                              aTTGCTGACCCTGGTGCTGGTTATAGTCACGTCAGTTACCGGTATGCCACTGATACGTAACCAGGGATTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGTTGAAATCCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGa                                                                                                                                                                                                                    <  1:526732/239‑1 (MQ=255)
                              aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGATTAAAAAGATTTTGCTGATCTTTTCTGCGTTGCTGCCTTTTTTTACGGGCGTTTTTTCCCCCTTCCCCGGGCTTTTTCTTTTTTTTccgccgtggtcaaaacacccg                                                                                                                                                                                                                      >  2:647178/1‑223 (MQ=255)
                              aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGGTCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGGCTTTTACTTTTTTTTCCGCCCGTGCAAAAACATCTg                                                                                                                                                                                                                     >  1:533747/1‑238 (MQ=255)
                                    taCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGa                                                                                                                                                                                                              <  1:4972/238‑1 (MQ=255)
                                     aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGGTGATTTAAAAAGTTTCGCCGATCATTGCTGGGTTGCTGCTTTTTTTTCCGGGCTTTTTTATCGACTTCCTCGGGCTTTTTCTTTTTTTTCCCCCCGGGCAAAAGAACCTGaaatggaa                                                                                                                                                                                                             >  2:90846/1‑232 (MQ=255)
                                          cGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGAGGTGATTTAAAATGGTTTGCTGGTCATTGCTGGTTTGGTGCTTTTATTTACGGGCTTTTTTTACGCCTTTCTCGGGCTTTTCCTTTTTTTTGCCCCGGTCGTAAAACACTTGCCAGTGAAGGTgg                                                                                                                                                                                                       >  2:571291/1‑238 (MQ=255)
                                             ggtGGTTTTTTTCACGTAAGTTTTGGGTATGGTATGTGGACGTAACCAGGGCTTTAAGAATTGCGTGCGGATGCAGCAAAAAATGGCGGCGGGTGAAAGCCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTGTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGc                                                                                                                                                                                                    <  2:20919/237‑1 (MQ=255)
                                                 gTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAAGGTTTTGCTGATCATTGCTGGGTTGCTGCTTTTTTTACCGGGGTTTTTTACCGACTTCCCCCGTCTTTTACTTTTTTTGCCCCCGGTGCCAAAGCATCTGACAGTGGAGTTTATGCCGCa                                                                                                                                                                                                >  1:660842/1‑239 (MQ=255)
                                                       tttccGTTAGTTTTCGGTTTGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGAAAAAAAGGGCGGGGGGTGAAAACCCAGGGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTCCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                                                                                                          <  2:281888/235‑1 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTTCCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAACACATTGTCAGTGAAGTTTATGCCGCATTTTGGCttt                                                                                                                                                                                      >  1:911200/1‑239 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCCGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCCCCGGTGCAAAAACATCTGCCCGTGAAGTTGATTCCGCCTTTtcggctt                                                                                                                                                                                      >  1:281888/1‑235 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCCTTTATTTACGGGCCTTTTTACCGACTTCCTCGGTCTTTTACTTTTATTGCCCCCGGTGCAAAAACATCTTACAGTGAAGTTGATTCCGCATTTTCGCttt                                                                                                                                                                                      >  1:825937/1‑239 (MQ=255)
                                                               aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgcc                                                                                                                                                                                                                                      >  1:642852/1‑189 (MQ=255)
                                                               aGTTATCGGTATGGCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgcc                                                                                                                                                                                                                                      <  2:642852/189‑1 (MQ=255)
                                                                 ttATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAATGGTTCGCTGATCATTGATGGTTTTGTGCTTTTATTACCCGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCGCCCGTGCCAAAAAATCTTCCAATGAAGTTTGTGTCGCAATTTGGCtttttttc                                                                                                                                                                                 >  2:323600/1‑234 (MQ=255)
                                                                      ggTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAAAGTTTCGCTGATCATTGCTGGTTTGCGGCTTTTATTACCGGGCTTTTTTTCCTCCTTCCTCGGGCTTCTACTTTTTTTTCCGCCGGTGCAAAAACACCTGACACTGAAGGTCATGCCCCATTTCTGCTTTTTTCGCAAGcc                                                                                                                                                                           >  2:160767/1‑239 (MQ=255)
                                                                           gTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCCTTTTTACCGACTTCCTCGGTCTTTTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGGGAAGGTGATGCCGCATTTTGGCTTTtcttcccttcctgggg                                                                                                                                                                      >  1:464176/1‑226 (MQ=255)
                                                                              aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGTAAAACCCCGCGGCGGAGGTGGTAAAAAGGTGTTCGCTGATCATTGCTGGTTTGTTGCTTTTATTACTGGGCTTTTTTACCGACTTCCCCGGGCTTTTACTTTTTTTTCCGCCGGGGCAAAAACATTTGCCCGTGAAGTTGATACCCCATTTCTGCTTTTCTTCCATGcgctgggggg                                                                                                                                                                   >  2:802067/1‑230 (MQ=255)
                                                                                    aCGTAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTa                                                                                                                                                            <  1:571291/239‑1 (MQ=255)
                                                                                             gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGGCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACCGTGAAGTTTATGCCGCACTTTCGCTTTTTTCGCATGCCTGGCGGCGGGTTTTGAGCCggg                                                                                                                                                   >  1:162169/1‑238 (MQ=255)
                                                                                             gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGATGATTAAAAGAGGTTCGCTGGTCATTGCTGGTTTTCTGCTTTTATTACCGGGCGTTTTTACCGACTTCCTCGGTCTTATACTTTTATTTCCGCCGTGGTAAAAACATCTGCAAGGGAGATTTAATCCCCACTTTGGGCTTTTTCGCGTGCCCGGCGGCGGGGTTATCGCCcggg                                                                                                                                                  >  2:522794/1‑238 (MQ=255)
                                                                                              ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGGTTGCTGCTTTTTTTATCGGGCTTTTTTACCGGCTTCCTCGGGCTTTTACTTTTTTTTCCGCCGGTGCCAAAGCACCTTGCAGGGAATGTGAGTCCCCATTTTCGCTTTTCTCGCATGCCCGGCGGCGGGTTTAAGCGCCGGGa                                                                                                                                                  >  2:31356/1‑239 (MQ=255)
                                                                                               gCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGa                                                                                                                                                  <  1:249991/237‑1 (MQ=255)
                                                                                                  tttataatttggggtggATGCAGCAAAAATGGGTGGGGGCGGAAACACCAGGTGCGGAGATGATTAAAAGTGTGTAGCTGATAGTGGCCGGTTTGCTGCTTTTATGACCGGGCTTTTATGACGACTTCCTCGGCCTTCTACTTGTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGGGCTTTTCTCGCATGACTGGCGGCGGTTTTAGCGCCGTGACCGGt                                                                                                                                             <  2:635937/223‑1 (MQ=255)
                                                                                                  ttAACAATTGCGTGTTGAGGCACAAAAAACTGGAGGCGGCTGAAAAGCCAGGGGATGTGATGTTGAAAAGAGTGTCGTTGTTCCTTGCTGGTTTGCGTTTTTTATTTCCGGGCTTTTTGACCGCCTTCATCGGTCTTCTTCTTTTATTGCCGCCGGTGCAAAAACATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                                                                                                             <  2:531153/239‑1 (MQ=255)
                                                                                                   tAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTTTTACCGGGCTTTTTTTCCGACTTCCTCGGGCTTTTACTTTTATTTCCGCCCGTGCAAAAACACCTGCCAGGGAAGTTGATGCCCCATTTTCGCCTTTTTCCCATGCCTGGGGGGGGGTTTATAGCCCGGGCCgggg                                                                                                                                            >  2:183625/1‑237 (MQ=255)
                                                                                                        aTTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTg                                                                                                                                                                                                               <  1:127530/170‑1 (MQ=255)
                                                                                                        aTTCCGTGCTGAGGCAGAAAAAAATGGCGGGGGGGGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTTCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGt                                                                                                                                       <  1:737931/239‑1 (MQ=255)
                                                                                                         tttCGTGGTGATGAAGCAAAAAATGGCGGCGGGTGAAACCCCAGCGGCGGAGATGATTAAAAGTGTTTTGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTa                                                                                                                                      <  1:837620/239‑1 (MQ=255)
                                                                                                           tcgtgcgtgtGCAGAAAAAAATGGGGGCGGGTGAAACCCCGGGGGGGGAGATGATTAAAAGTGTTTCCCTGACCATTGCTGGTTTGCTGTTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTTCTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAAt                                                                                                                                    <  2:464176/230‑1 (MQ=255)
                                                                                                            cGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTTGCGCCCGGGACCGTGGGGGGTata                                                                                                                                    >  1:173737/1‑237 (MQ=255)
                                                                                                                  gaggcagcaAAAAATGGCGGCGGGTGAAACCCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                                                                                                             <  1:160767/236‑1 (MQ=255)
                                                                                                                       tcaaaaaatTGGCGGCGGGTTAAAACCCGGGGGCGGAGTTTATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:660842/236‑1 (MQ=255)
                                                                                                                        gaaaaaaaTGGGGGGGGGTGAAACCCCGGCGGGGGATATGATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:247574/236‑1 (MQ=255)
                                                                                                                            aaaaTGGGGGCGGGTGAAAACCCGGGGGTGGAGATGATTAAAAGTTTTTCCGTGATTATTGCGTGTTTTCTGTTTTTATTACCGGGCTTTTTGACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTa                                                                                                                   <  2:162169/239‑1 (MQ=255)
                                                                                                                             aaaTGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                                                                                                                  <  1:653680/239‑1 (MQ=255)
                                                                                                                                  gggggggggtaaacCCCCGCGGTGGAGTTGATTAAAGGTGTTCTGCTGAGCATTGCTGGTTTTCTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  2:38883/229‑1 (MQ=255)
                                                                                                                                  gcggcgGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  1:183625/239‑1 (MQ=255)
                                                                                                                                   cgggggGTGAAACCCCATGGGGGGGTATGTTTAAAAGTTTTTCCCTTATCATTGCTGTGTTTCTGCTTGTATTACCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTTCAAAAGCATCTTACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  2:204441/238‑1 (MQ=255)
                                                                                                                                                   aGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTTACGCCGGGAGCCGGGGCGGGGATACTTTTTATGGTGAGGACCAGCGAAAGGGAGATGAggggg                                                                                            >  2:517579/1‑235 (MQ=255)
                                                                                                                                                    ggggggggttttgtttaaagttttttCGCTGATCTTTGCTGGTTTTCTGTTTTTATTCCCGGGCTTTTTGCCGGACTTCCTCGGTCTTTTATTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGGTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACTAGCGAAAGGATGATGAGcgcg                                                                                            <  2:173737/216‑1 (MQ=255)
                                                                                                                                                    gcggcgGGTATTATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                           <  1:473478/239‑1 (MQ=255)
                                                                                                                                                    gcggcgGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                           <  1:132704/239‑1 (MQ=255)
                                                                                                                                                             gattatTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                                                                                    <  1:323600/238‑1 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        >  1:352647/1‑188 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:352647/188‑1 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCCTAAAGACGa                                                                     >  1:632078/1‑238 (MQ=255)
                                                                                                                                                                              cGCTGATCGTTGCTTGTTTGCTGCTTTTATTACCGGGCTTTTTTCCCGATTTCCTCGGTCTGCTGCTTTTATGGCCTCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCGGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGc                                                                 <  1:66480/239‑1 (MQ=255)
                                                                                                                                                                               gCTGATCTTGGCTGGTTTGCGGCTTTTTTTCCCGGGTTTTTTTACCGACTTCCGCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTCAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                                                <  1:31356/239‑1 (MQ=255)
                                                                                                                                                                                   aTCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGGGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAgag                                                            >  1:867231/1‑237 (MQ=255)
                                                                                                                                                                                     cATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGTGTACCAGCGAAAGGATGATGACGGCGACCGCCTTGATCATAAAGACGATCGCCAGGGtt                                                          >  1:840498/1‑239 (MQ=255)
                                                                                                                                                                                        tGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGTTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGCCCGGTGGCGGGAATACTTTTTTTGGGGAGTACCATCGAAAGGATGAGTAGCGCGCCCGCCCTTATCACAAAAAAGAGCGCCCggggtaag                                                       >  2:751622/1‑233 (MQ=255)
                                                                                                                                                                                             ttttgttttttttttTGCGGGGTTTTTTTCCCGTCTTCCTCGGTCTTCTACTTTTATTCCCGCCGGTGCAAAAGGATCTGACAGGGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCg                                                  <  2:533747/229‑1 (MQ=255)
                                                                                                                                                                                             gTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGAAGATGAGCGGGACCGCCCTGATCATAAAGACGATCGCCAGGGTTAATGTCGa                                                 >  1:43801/1‑239 (MQ=255)
                                                                                                                                                                                                tgctgcTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCCGGATTAATGTCGaaa                                               >  1:715531/1‑238 (MQ=255)
                                                                                                                                                                                                 gttttttttttttCCGGGTTTTTTTCCCGACTTCCTCGGTCTTCTACTTTTTGTGCCCCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                             <  2:715531/234‑1 (MQ=255)
                                                                                                                                                                                                 gctgctTTTATTCCCGGTCTTTTTTCCCGACTTCCGCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACg                                             <  1:522794/239‑1 (MQ=255)
                                                                                                                                                                                                  ctgctTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAACCGATCGCCAGGATTAATGTCGcaaag                                             >  1:599014/1‑236 (MQ=255)
                                                                                                                                                                                                   tttttttATTTCCGGGTTTTTTGACCGCCGTCCCCGTGCTTTTCCTTTGTTGGCCGCCGGTGCAAAAGCATCTGTCGGGGAAGGTGATGCCGCATTTGAGCTTTTCTCGCATGCCGGGGGGCGGTTTTAGCGCTGGGACCGGTGGCGGTAAGCCTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                           <  2:261016/236‑1 (MQ=255)
                                                                                                                                                                                                    tttttttttctccggtttttttACCCCCTTCCCCGGGTTTTTTTTTTTTGTGCCCGCCGGTGCAAAAGCATCTGCCAGTGAATTTGATGCCGCATTTGCGCTTTTCTCGCAGGCCTGGGGGCGTTTTTGGCGCCGGGACCGGTGGCGGTAATCCTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGcc                                           <  2:867231/223‑1 (MQ=255)
                                                                                                                                                                                                               cGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGGGGGTCCAGCGGAAGGAGGATTGGCGCGAACCCCTTGTTCATAAAAAACAATCCCAGGATTAATGTTCAAAACGCCggttatagggg                                >  2:304309/1‑230 (MQ=255)
                                                                                                                                                                                                                  gCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATGGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTa                             <  1:517579/238‑1 (MQ=255)
                                                                                                                                                                                                                      ttttACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGGGGTACCAGGGAAAGGGTGAGGAGGGCGACCGCCTTGAGCATAAAGACGATCGCCAGGATTAATGTCGAAACACCGGGTTATGTGGGTATGCCa                        >  2:203637/1‑239 (MQ=255)
                                                                                                                                                                                                                           ccGATTTCCTCGGTCTTCTTCTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTc                   <  1:90846/239‑1 (MQ=255)
                                                                                                                                                                                                                                             ctACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCCCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGCTTATGTGGTTATGCCATTTTTCCGGGGTTTTCGTTTTg  >  2:842637/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
CTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTG  >  NC_000913/4368721‑4369194

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