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

CGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAA  >  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)
                                          gAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGtt                                                                                                                                                                                            <  1:524190/251‑1 (MQ=255)
                                                 cGGCGCGTAGCGGGATGGAAGTCTTTTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGCCCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCa                                                                                                                                                                                      <  1:247411/250‑1 (MQ=255)
                                                        tAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCt                                                                                                                                                                               >  1:457226/1‑250 (MQ=255)
                                                             ggATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGCTTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGAACTCTCACGCCAACAAAGGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCTCCACTGC‑CAACTGGTTACGTAATACGGCTAAAcc                                                                                                                                                                          >  2:180443/1‑250 (MQ=255)
                                                              gATTGAAGTCGTGTCGGTTTATCAGCCCCGGCACCTGCTTCAGGGCCAGGCAATGTGGATTAAAGGATGGCTGGCCCGTGTGGATTTAGATAACGCCCGTTGCCAGCTGAACCGTCGGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATTCCGGGTGAGTTTCGGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCag                                                                                                                                                                        <  2:509981/251‑1 (MQ=255)
                                                                  gaaggtttttctctttATCACCAACAGTAATAGCTTCAGTGCCTTGCAATGTTGATTCAATTCTCGCTGGCCCGTGGTGATCTAGATAACTCCCGTAGCCACCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGt                                                                                                                                                                     <  1:275753/236‑1 (MQ=255)
                                                                            cgtcTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGGAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGaacaac                                                                                                                                                          <  1:219228/251‑1 (MQ=255)
                                                                               cTTATAAGCCAGAGCAACAGCTTCAGAGCATGGCAATGTTGATTCAATGCTCACTGGACAGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGCAAACCTGCTGGGACATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTAATGGCAACTGACAGGCGACAAAGCCGCCGATGCACAA‑TGGTTGCGTCATCAGGCTAAACCACAGTTAGCGacaacaca                                                                                                                                                        >  2:436795/1‑243 (MQ=255)
                                                                               cTTATAAGCCACAGCAACAGCATCAGTGTCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACAGGTTCTCTAACGACAACAAATGCCGGGTGATTTACGGGCAGATGACCGGAGATAAACCCGCCGCTGC‑CCACTGGTTGCGTCATACGGCTAAACCCGAGTTAGCGAACAACCAc                                                                                                                                                       >  2:632120/1‑251 (MQ=255)
                                                                                     aGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGATTGCGAACAACCACTTCCTg                                                                                                                                                 >  1:625314/1‑251 (MQ=255)
                                                                                             ccacatcttcatggcttgtcAATGTTGAGTCAATGCTCGTTGGCCCGTGTTGTTTTAGATAACGTCCGTAGCCAGCTTAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGACCTCTAACTCCAACAAATGCCGGGTGATTTACTGGCAAATGTCCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTGATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCTAGGTCa                                                                                                                                         <  2:127271/232‑1 (MQ=255)
                                                                                              aaCAGCTTCCGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAAAGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGTCTGGCTCTCTAACGCCAACAAATGCCGGGTGCTTTACTGGCAAATGACCGGCGATAAAGCCGCCGCAGC‑CAACTGGTTGCGTCATACGGCTAAAACAGAGTTTGCGAACAACCAATTCCTTCAAGGTcgg                                                                                                                                        >  2:18822/1‑249 (MQ=255)
                                                                                                 aGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGACGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATgg                                                                                                                                     >  1:78804/1‑251 (MQ=255)
                                                                                                      cAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTATCCAGCTGAACCGTCTGGAAAACCTGCTGTGGAATGGCAAATATCCCAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATCACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCaa                                                                                                                                <  1:635106/251‑1 (MQ=255)
                                                                                                           ccTGGCAATGTTGATTCAGTGCTCGCTGGACTGTGGTGATTTAGATAAGGCCCGTAGCCATCTGAGCCGGCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACTCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTg                                                                                                                           <  1:284884/251‑1 (MQ=255)
                                                                                                            cTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGCACAACCACTTCCTGCGAGGTCAATGGCGCAACATTGc                                                                                                                          >  2:307505/1‑251 (MQ=255)
                                                                                                             tGGCAATGTTGATTCAATGCTAGCTGGCCCGTGGTGATTTAGATAACGACCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATAACCGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGACGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCCGAGTTTGCGACCAAACACTTCCTGCAAGGTCAATGgcgaaaacttgcc                                                                                                                         >  2:639489/1‑241 (MQ=255)
                                                                                                                     gtaattaaaagtccGCTGGCCGGTGGTGTTTTAGAGAACGCCCTTAGCCAAATAAACCGGCTTCAAAACCTAAGGTGGAATGGCAAAGATCACAGCGACTGGACCTCTAACAACAACAAAAGCAGGGTGATTGACTGGCAAATGACCGGCGATAAAGCCGCCGCTAC‑CCACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGTAAGGCCAATGGCGCAACAtataccttgcaca                                                                                                                 <  2:665889/238‑13 (MQ=255)
                                                                                                                           cAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCtt                                                                                                           >  1:161621/1‑251 (MQ=255)
                                                                                                                            aaTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTg                                                                                                          <  1:646650/251‑1 (MQ=255)
                                                                                                                             aTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGc                                                                                                         <  1:632120/251‑1 (MQ=255)
                                                                                                                             aTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGc                                                                                                         >  1:254299/1‑251 (MQ=255)
                                                                                                                                  cGCTGGCCCGTGGGTTTTTAGATAACGCCCATAGCCAGCTTAACCGTCTGGAAAACCGTCTGGGGAAGGGCAAATATCACAGCGTCTGGATCTCTAACGCCAACAAATGCCGTGTAATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGACCGTGCACAAATCTTGCtt                                                                                                       <  2:380862/248‑2 (MQ=255)
                                                                                                                                          cgtgggaaTTTAGATAACTCCTGTAGCCATCTGAACCGTCTGTAAAACCTGCTTGGGAATGTCAAAGATCACAGCGACGGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGTCGATAAAGCCGCCGCTGC‑CAACTGGTTTCGGCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGa                                                                                            <  1:661936/244‑1 (MQ=255)
                                                                                                                                          cGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGa                                                                                            <  1:141978/251‑1 (MQ=255)
                                                                                                                                                 aTTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGa                                                                                                                                                                                                                                           <  1:30303/102‑1 (MQ=255)
                                                                                                                                                 aTTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGa                                                                                                                                                                                                                                           >  2:30303/1‑102 (MQ=255)
                                                                                                                                                            gCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCtt                                                                                                           <  1:479813/218‑1 (MQ=255)
                                                                                                                                                            gCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTtct                                                                          <  1:639489/251‑1 (MQ=255)
                                                                                                                                                               cGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCg                                                                        >  1:389617/1‑250 (MQ=255)
                                                                                                                                                                    ccAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAAc                                                                  >  1:430241/1‑251 (MQ=255)
                                                                                                                                                                     cAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGACCAACAACTTCCTGCAAGGTCACTGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaagacct                                                                 >  2:501964/1‑248 (MQ=255)
                                                                                                                                                                       gCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTc                                                                >  1:595262/1‑250 (MQ=255)
                                                                                                                                                                            aCCGTCGGGAAAACCTGCTGGGGAAGGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGGCCGGGTGATTTACTGGCAAATGACCGGCGATAAATCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGaa                                                          <  2:286301/251‑1 (MQ=255)
                                                                                                                                                                             ccGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTCCTGGCCAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACCAATCTTGCTTGGCGCGTTTGAACCGGCAGAAATTGTTCTCgaag                                                                     >  2:530532/1‑239 (MQ=255)
                                                                                                                                                                             ccGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCgaag                                                                     <  1:530532/239‑1 (MQ=255)
                                                                                                                                                                              cgcctagAAAACCAGCTGGGGAATGGCAAATATCACATCGCCTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGCGTTTGCGACCAACCACTTCCTGCATGGCCAAAGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGTAATTGGTCTCGACGAAATCaacaaaa                                                         >  2:248332/7‑245 (MQ=255)
                                                                                                                                                                              cGTCTGGAAAACCTGCGGGGGAATGGCAAATATCACAGCGATTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGaaaa                                                        <  2:132803/251‑1 (MQ=255)
                                                                                                                                                                               gTCTGGAAAACCTGCTGGGGAAGGGTAAATAAGACAGCGACTGGATCTCTAACGCCAACAAAGGCCGGGTGATTTACGGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAAt                                                       <  2:208306/251‑1 (MQ=255)
                                                                                                                                                                                     aaaaCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCgg                                                 >  2:158786/1‑251 (MQ=255)
                                                                                                                                                                                               ggggAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCAGGAGTCTGCGGt                                       >  2:53989/1‑251 (MQ=255)
                                                                                                                                                                                                 ggAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtg                                     <  1:158786/251‑1 (MQ=255)
                                                                                                                                                                                                               cacaGCGACTGGATCTCTACCGCCAACAGATGCCGGGTGATTATCTGGAAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGGTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCa                        <  2:305313/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCg                    >  1:346889/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACACCCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAAcc               >  2:592389/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGCGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGtt           >  2:143858/1‑251 (MQ=255)
                                                                                                                                                                                                                               ctAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTtgctg       >  1:96355/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCCAACAAAGCCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTaa  <  2:96355/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                 
CGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGC‑CAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGCTTAA  >  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: ‑