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,0470GT84.6% 107.2 / ‑2.2 53V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base T (22/22);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,555,0480TG92.3% 136.2 / ‑1.6 53V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base G (23/25);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAA  >  NC_000913/3554810‑3555288
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                 
ggggcgtagttgtgaggggCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTTGTGTCGGCTTATCCGCCACAGCATCAGCTTGAGGGCCTGGCAATGTTGATGCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGGCCGGGTGATTTa                                                                                                                                                                                                                                       <  2:662366/235‑1 (MQ=255)
           gTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAAt                                                                                                                                                                                                                              <  1:311328/249‑1 (MQ=255)
                gggCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCgg                                                                                                                                                                                                                        <  1:483007/250‑1 (MQ=255)
                      ggggtcgGCTGGACGAAGCCGATGCGTGGGCGCGTAGCGGGATTGAGGTCTTTGCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTTATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATa                                                                                                                                                                                                                   <  2:189927/244‑1 (MQ=255)
                            ggcttgaggAAGCTGAAGCGTGGGGGCGTAGCGGGATTGAAGTCTTGTCGTATTAGCAGCCACAGCAACAGCTTCAGTGCCTGTCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGTCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAgccgcc                                                                                                                                                                                                           <  2:603388/243‑1 (MQ=255)
                                          gAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGtt                                                                                                                                                                                            <  1:524190/251‑1 (MQ=255)
                                                 cGGCGCGTAGCGGGATGGAAGTCTTTTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGCCCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCa                                                                                                                                                                                      <  1:247411/250‑1 (MQ=255)
                                                        tAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCt                                                                                                                                                                               >  1:457226/1‑250 (MQ=255)
                                                             ggATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGCTTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGAACTCTCACGCCAACAAAGGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCTCCACTGCCAACTGGTTACGTAATACGGCTAAAcc                                                                                                                                                                          >  2:180443/1‑250 (MQ=255)
                                                              gATTGAAGTCGTGTCGGTTTATCAGCCCCGGCACCTGCTTCAGGGCCAGGCAATGTGGATTAAAGGATGGCTGGCCCGTGTGGATTTAGATAACGCCCGTTGCCAGCTGAACCGTCGGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATTCCGGGTGAGTTTCGGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCag                                                                                                                                                                        <  2:509981/251‑1 (MQ=255)
                                                                  gaaggtttttctctttATCACCAACAGTAATAGCTTCAGTGCCTTGCAATGTTGATTCAATTCTCGCTGGCCCGTGGTGATCTAGATAACTCCCGTAGCCACCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGt                                                                                                                                                                     <  1:275753/236‑1 (MQ=255)
                                                                            cgtcTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGGAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGaacaac                                                                                                                                                          <  1:219228/251‑1 (MQ=255)
                                                                               cTTATAAGCCAGAGCAACAGCTTCAGAGCATGGCAATGTTGATTCAATGCTCACTGGACAGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGCAAACCTGCTGGGACATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTAATGGCAACTGACAGGCGACAAAGCCGCCGATGCACAATGGTTGCGTCATCAGGCTAAACCACAGTTAGCGacaacaca                                                                                                                                                        >  2:436795/1‑243 (MQ=255)
                                                                               cTTATAAGCCACAGCAACAGCATCAGTGTCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACAGGTTCTCTAACGACAACAAATGCCGGGTGATTTACGGGCAGATGACCGGAGATAAACCCGCCGCTGCCCACTGGTTGCGTCATACGGCTAAACCCGAGTTAGCGAACAACCAc                                                                                                                                                       >  2:632120/1‑251 (MQ=255)
                                                                                     aGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGATTGCGAACAACCACTTCCTg                                                                                                                                                 >  1:625314/1‑251 (MQ=255)
                                                                                             ccacatcttcatggcttgtcAATGTTGAGTCAATGCTCGTTGGCCCGTGTTGTTTTAGATAACGTCCGTAGCCAGCTTAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGACCTCTAACTCCAACAAATGCCGGGTGATTTACTGGCAAATGTCCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTGATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCTAGGTCa                                                                                                                                         <  2:127271/232‑1 (MQ=255)
                                                                                              aaCAGCTTCCGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAAAGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGTCTGGCTCTCTAACGCCAACAAATGCCGGGTGCTTTACTGGCAAATGACCGGCGATAAAGCCGCCGCAGCCAACTGGTTGCGTCATACGGCTAAAACAGAGTTTGCGAACAACCAATTCCTTCAAGGTcgg                                                                                                                                        >  2:18822/1‑249 (MQ=255)
                                                                                                 aGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGACGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATgg                                                                                                                                     >  1:78804/1‑251 (MQ=255)
                                                                                                      cAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTATCCAGCTGAACCGTCTGGAAAACCTGCTGTGGAATGGCAAATATCCCAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATCACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCaa                                                                                                                                <  1:635106/251‑1 (MQ=255)
                                                                                                           ccTGGCAATGTTGATTCAGTGCTCGCTGGACTGTGGTGATTTAGATAAGGCCCGTAGCCATCTGAGCCGGCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACTCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTg                                                                                                                           <  1:284884/251‑1 (MQ=255)
                                                                                                            cTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGCACAACCACTTCCTGCGAGGTCAATGGCGCAACATTGc                                                                                                                          >  2:307505/1‑251 (MQ=255)
                                                                                                             tGGCAATGTTGATTCAATGCTAGCTGGCCCGTGGTGATTTAGATAACGACCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATAACCGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGACGCTGCCAACTGGTTGCGTCATACGGCTAAACCCGAGTTTGCGACCAAACACTTCCTGCAAGGTCAATGgcgaaaacttgcc                                                                                                                         >  2:639489/1‑241 (MQ=255)
                                                                                                                     gtaattaaaagtccGCTGGCCGGTGGTGTTTTAGAGAACGCCCTTAGCCAAATAAACCGGCTTCAAAACCTAAGGTGGAATGGCAAAGATCACAGCGACTGGACCTCTAACAACAACAAAAGCAGGGTGATTGACTGGCAAATGACCGGCGATAAAGCCGCCGCTACCCACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGTAAGGCCAATGGCGCAACAtataccttgcaca                                                                                                                 <  2:665889/238‑13 (MQ=255)
                                                                                                                           cAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCtt                                                                                                           >  1:161621/1‑251 (MQ=255)
                                                                                                                            aaTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTg                                                                                                          <  1:646650/251‑1 (MQ=255)
                                                                                                                             aTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGc                                                                                                         <  1:632120/251‑1 (MQ=255)
                                                                                                                             aTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGc                                                                                                         >  1:254299/1‑251 (MQ=255)
                                                                                                                                  cGCTGGCCCGTGGGTTTTTAGATAACGCCCATAGCCAGCTTAACCGTCTGGAAAACCGTCTGGGGAAGGGCAAATATCACAGCGTCTGGATCTCTAACGCCAACAAATGCCGTGTAATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGACCGTGCACAAATCTTGCtt                                                                                                       <  2:380862/248‑2 (MQ=255)
                                                                                                                                          cgtgggaaTTTAGATAACTCCTGTAGCCATCTGAACCGTCTGTAAAACCTGCTTGGGAATGTCAAAGATCACAGCGACGGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGTCGATAAAGCCGCCGCTGCCAACTGGTTTCGGCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGa                                                                                            <  1:661936/244‑1 (MQ=255)
                                                                                                                                          cGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGa                                                                                            <  1:141978/251‑1 (MQ=255)
                                                                                                                                                 aTTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGa                                                                                                                                                                                                                                           <  1:30303/102‑1 (MQ=255)
                                                                                                                                                 aTTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGa                                                                                                                                                                                                                                           >  2:30303/1‑102 (MQ=255)
                                                                                                                                                            gCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCtt                                                                                                           <  1:479813/218‑1 (MQ=255)
                                                                                                                                                            gCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTtct                                                                          <  1:639489/251‑1 (MQ=255)
                                                                                                                                                               cGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCg                                                                        >  1:389617/1‑250 (MQ=255)
                                                                                                                                                                    ccAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAAc                                                                  >  1:430241/1‑251 (MQ=255)
                                                                                                                                                                     cAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGACCAACAACTTCCTGCAAGGTCACTGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaagacct                                                                 >  2:501964/1‑248 (MQ=255)
                                                                                                                                                                       gCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTc                                                                >  1:595262/1‑250 (MQ=255)
                                                                                                                                                                            aCCGTCGGGAAAACCTGCTGGGGAAGGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGGCCGGGTGATTTACTGGCAAATGACCGGCGATAAATCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGaa                                                          <  2:286301/251‑1 (MQ=255)
                                                                                                                                                                             ccGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTCCTGGCCAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACCAATCTTGCTTGGCGCGTTTGAACCGGCAGAAATTGTTCTCgaag                                                                     >  2:530532/1‑239 (MQ=255)
                                                                                                                                                                             ccGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaag                                                                     <  1:530532/239‑1 (MQ=255)
                                                                                                                                                                              cgcctagAAAACCAGCTGGGGAATGGCAAATATCACATCGCCTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGCGTTTGCGACCAACCACTTCCTGCATGGCCAAAGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGTAATTGGTCTCGACGAAATCaacaaaa                                                         >  2:248332/7‑245 (MQ=255)
                                                                                                                                                                              cGTCTGGAAAACCTGCGGGGGAATGGCAAATATCACAGCGATTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGaaaa                                                        <  2:132803/251‑1 (MQ=255)
                                                                                                                                                                               gTCTGGAAAACCTGCTGGGGAAGGGTAAATAAGACAGCGACTGGATCTCTAACGCCAACAAAGGCCGGGTGATTTACGGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAAt                                                       <  2:208306/251‑1 (MQ=255)
                                                                                                                                                                                     aaaaCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCgg                                                 >  2:158786/1‑251 (MQ=255)
                                                                                                                                                                                               ggggAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCAGGAGTCTGCGGt                                       >  2:53989/1‑251 (MQ=255)
                                                                                                                                                                                                 ggAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtg                                     <  1:158786/251‑1 (MQ=255)
                                                                                                                                                                                                               cacaGCGACTGGATCTCTACCGCCAACAGATGCCGGGTGATTATCTGGAAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGGTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCa                        <  2:305313/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCg                    >  1:346889/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACACCCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAAcc               >  2:592389/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGCGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGtt           >  2:143858/1‑251 (MQ=255)
                                                                                                                                                                                                                               ctAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTtgctg       >  1:96355/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCCAACAAAGCCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTaa  <  2:96355/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                 
CGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAA  >  NC_000913/3554810‑3555288

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