Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,555,047 2 bp→TG coding (1963‑1964/2706 nt) malT → DNA‑binding transcriptional activator MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,555,0470GT92.3% 123.0 / 2.0 52V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (17/31);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01
*NC_0009133,555,0480TG94.2% 131.4 / 2.1 52V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (17/32);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

CGCATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAAT  >  NC_000913/3554804‑3555289
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
cgcATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGt                                                                                                                                                                                                                                              <  1:493335/251‑1 (MQ=255)
    atCGTGCGCAGGTGGTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAAGATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGa                                                                                                                                                                                                                                            <  2:435557/248‑1 (MQ=255)
                 gTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACTGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTTCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATg                                                                                                                                                                                                                              <  1:500433/250‑1 (MQ=255)
                 gTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGGCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGGGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGa                                                                                                                                                                                                                             <  1:250779/251‑1 (MQ=255)
                       ggctcgtggggggcgggagggaGCCGAAGCGGGGGCGCGTGGCGGGATTGAAGGTTTGGAGTGTGATCAGCCAAAGCAACAGCTGCAGTGCCTGGCAATGGTGATTAAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGGAGCCAGATGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGAATGGAGCGCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCg                                                                                                                                                                                                                       <  2:182170/230‑1 (MQ=255)
                         ccTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGAt                                                                                                                                                                                                                     <  1:114810/251‑1 (MQ=255)
                          cTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGAt                                                                                                                                                                                                                     >  1:130464/1‑250 (MQ=255)
                                      ggATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTTAAGTCTTGTCGTCTTACCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTg                                                                                                                                                                                                        <  1:490455/251‑1 (MQ=255)
                                                  gccgccgccgcgTGGCGGGTTTGAAGTCTTGTCGTCTTATCAGCCACAGAAACAGCTTCAGTGCCTGGCAATGGGTATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGGAGCCAGCTGAACCGTCTGGAAAACCTGCGGGGTAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGTTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGc                                                                                                                                                                                           <  2:503472/243‑1 (MQ=255)
                                                         gcgcGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATGCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATAc                                                                                                                                                                                    <  1:44754/251‑1 (MQ=255)
                                                                       tGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCagag                                                                                                                                                                       <  1:368624/250‑1 (MQ=255)
                                                                         aaGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGGTCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGttt                                                                                                                                                                    >  1:62957/1‑251 (MQ=255)
                                                                               tGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGaac                                                                                                                                                              >  1:361764/1‑251 (MQ=255)
                                                                                    tcTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCACTAACGCTAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAAccc                                                                                                                                                         >  2:151104/1‑250 (MQ=255)
                                                                                      ttATCAGCCACAGCAACAGCTTCATTGCCTGGCAATGTTGATTCAATGCGCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCAc                                                                                                                                                        <  2:34009/250‑1 (MQ=255)
                                                                                      ttATCAGCCACAGCAACAGCTTCAGGGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCAc                                                                                                                                                        <  1:442596/250‑1 (MQ=255)
                                                                                       tATCAGCCACAGCCACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTAGCTGGCCCGAGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGTAAAACCTGCTGGGGAATGGCAAGTATCACAGCGACTGGATCTCTACCGCCAACGAATGCCGGGTGATTTACTGGCAAATAACCGGCCCTACAGCCGCAGCTTCCAACTTGGTTCGTCATACGGCTAACCCAGATTTTTCGAACAACCACtt                                                                                                                                                      >  2:488612/1‑251 (MQ=255)
                                                                                          cagcaacagcaaCAGCTTCAGTGACTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAACTATCACAGCGACTGGATCTCTAAAGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGAGATAAAGCCGCCTCGGCCAACTGGTTTCGTCATACGGATAAACCAGAGTTTGCGAACAACCACTTCCt                                                                                                                                                   >  2:394470/1‑251 (MQ=255)
                                                                                                caGCAACAGCTGCAGTGCAGGGCAAAGTTGATTAAAGGCTCGCTGGCCCGTGGTGAGTTAGAAAACCCCCGTGGCCAGATGAACCGTCTGTAAAACCTGCTGGGGAAGGGCAAAGATCATACCGACTGGATCTCTAACGCCAACAAAGGCCGGGGGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCCCTTCCTGCAAgg                                                                                                                                             <  2:241833/251‑1 (MQ=255)
                                                                                                       aGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATgg                                                                                                                                      >  1:209823/1‑251 (MQ=255)
                                                                                                                  cTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGCAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCCATGGCGCAACATTGc                                                                                                                           >  1:345758/1‑251 (MQ=255)
                                                                                                                    tgCAATGTTGATCCAATGCCCGCTGTCCTGTGGTGATTTAGATAACGCTCGTAGCCATCTGAACCGTCTGGAAAACCTGCGGGGGAATGGCAAATATCACAGCGACTGGCTCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAAGGACCGGCGATAAAGCCGCCGATGACAACTGGTTGCGTCATACGGCTAAACCAGAGTTTTCGAACAACCTCTTCCTGCAAGGTCAATGGCGCAACATTGcc                                                                                                                          <  1:11956/249‑1 (MQ=255)
                                                                                                                        attttGATTCAATGCTCGCTGGCACGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCCCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACACCCACTTCCGGCAAGGTCAATGGCGCAACATTGCCCGTGc                                                                                                                     >  2:11326/4‑251 (MQ=255)
                                                                                                                            tGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGcaca                                                                                                                  <  1:441551/250‑1 (MQ=255)
                                                                                                                                  aaTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTg                                                                                                           <  1:475989/251‑1 (MQ=255)
                                                                                                                                  aaTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAAAATTGCCCGTGCCCAAATCTTg                                                                                                           >  2:169887/1‑251 (MQ=255)
                                                                                                                                        cgcgtgccggtgtggATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGc                                                                                                     <  2:103805/237‑1 (MQ=255)
                                                                                                                                                cGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGa                                                                                             <  1:291606/251‑1 (MQ=255)
                                                                                                                                                  tggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGACGCCGCTGCCAACTGGTTGCGTCATACGGCTCAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGACCATGGCGAAACATTGCCCGTGCACAAAACTTGCTGGGCTGGTTTGAAc                                                                                           >  2:323936/1‑251 (MQ=255)
                                                                                                                                                       aTTTAGATAACGCCAGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGc                                                                                       >  1:408387/1‑250 (MQ=255)
                                                                                                                                                         ttAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATTACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGc                                                                                                                                                 <  1:177122/190‑1 (MQ=255)
                                                                                                                                                         ttAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGc                                                                                                                                                 >  2:177122/1‑190 (MQ=255)
                                                                                                                                                           aGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCACTAACGCTAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGaa                                                                                   <  1:151104/250‑1 (MQ=255)
                                                                                                                                                                  gCCCGTAGCCAGCTGAACCGGCTGGAAAACCTGCTGGGGAATGGCAAATATCCCAGCGACTGGAGCCCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACGTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTtct                                                                           <  2:408387/251‑1 (MQ=255)
                                                                                                                                                                    ccGTAGCCAGCGGAACCGTCTGGAAAACCTGCTGGGGAATGGAAAATATCACAGCGACGGGATCTCTAACGCCAACAAATGCCGGTGGATTTACTGGCAAAGGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCGTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCg                                                                         <  2:231076/251‑1 (MQ=255)
                                                                                                                                                                      gTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGATGCCAACTGGTTGCGTCATACGGCTAAACCAGAGATTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCga                                                                        >  1:100584/1‑250 (MQ=255)
                                                                                                                                                                        aGCCAGCTGAACCGTTTGGAAAACCTGCGGGGGAAGGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCAGACGGCGAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCGGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaaga                                                                     <  2:431624/251‑1 (MQ=255)
                                                                                                                                                                         gCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAACTCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaaga                                                                     >  2:133617/1‑250 (MQ=255)
                                                                                                                                                                          ccAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAAc                                                                   <  1:219916/251‑1 (MQ=255)
                                                                                                                                                                               ggacccggcTGGAAAACCGGCGGGGGAATTGAAAAGATCACAGCGACTGGATCTTTAACGCCAACAAATGCCGGGTTTTTTACTGGAAAATGACCGGGAATAAAGCCGACGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGGTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAATAACTCAAt                                                              <  2:314402/243‑1 (MQ=255)
                                                                                                                                                                                gAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCGCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCAGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCaa                                                               <  1:31188/249‑1 (MQ=255)
                                                                                                                                                                                 aaCCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATg                                                             >  2:174033/1‑250 (MQ=255)
                                                                                                                                                                                         ggAAAACCTGCTGGGGAATGGTAAAGATCATAGCGACTGGATCCCTAACGCCAACAAATGCCGGGGGATTTACGGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCACCATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGccc                                                    <  1:488612/251‑1 (MQ=255)
                                                                                                                                                                                          gAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCg                                                   <  1:397295/251‑1 (MQ=255)
                                                                                                                                                                                           aaaaCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGACCTCAATGAAAATGccagg                                                  >  2:189110/1‑248 (MQ=255)
                                                                                                                                                                                               ccTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCTCTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGATATTGTTCTCGAAGAACTCAATGAAAATTCCCGGAGTc                                              <  1:392699/251‑1 (MQ=255)
                                                                                                                                                                                                  gctgGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTg                                            <  1:344076/250‑1 (MQ=255)
                                                                                                                                                                                                            ggCAAATATCACAGCGACTGGATCGCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTTCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtgatga                                  <  1:472855/250‑1 (MQ=255)
                                                                                                                                                                                                                aaTATCACAGCGACTGGATCTCAAACGCCAACAAATGCCGGGGGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGAt                             <  2:361764/251‑1 (MQ=255)
                                                                                                                                                                                                                  aatCAGAGCAACTGGATCTCTAACGCCAACAAATGCCGGGTGATTGACTGGCAAATGACCGGCGAGAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGAtct                           <  2:175810/250‑1 (MQ=255)
                                                                                                                                                                                                                            aCTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACACATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAAc                 >  2:195833/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAGCAACCACTTCATGCAAGGTCAATGGCGCAACATTGCCCGTGCAAAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGCGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTact  >  1:20341/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
CGCATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAAT  >  NC_000913/3554804‑3555289

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