Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,415 Δ1 bp 7.3% intergenic (+19/+167) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,330,4150A. 7.3% 616.3 / 38.0 179intergenic (+19/+167)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (102/64);  new base . (4/9);  total (106/73)
Fisher's exact test for biased strand distribution p-value = 4.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 179 total aligned reads displayed.

GCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGA  >  NC_000913/3330270‑3330560
                                                                                                                                                 |                                                                                                                                                  
gCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACAt                                                                                                                                                 >  3:697600/1‑149 (MQ=255)
gCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACAt                                                                                                                                                 >  2:629773/1‑149 (MQ=255)
  aaCGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGTCACATg‑g                                                                                                                                              >  4:378996/1‑149 (MQ=255)
      aaTGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCgg                                                                                                                                          <  2:921150/149‑1 (MQ=255)
        tGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                        <  4:1076759/149‑1 (MQ=255)
        tGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                        >  3:1076759/1‑149 (MQ=255)
        tGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATAAGTCAAAAAGAAACCCCGGCACATG‑GCCgggg                                                                                                                                        <  2:1234288/149‑1 (MQ=255)
         ggCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGc                                                                                                                                       <  1:29231/149‑1 (MQ=255)
            gTTTGTGCAATATCTTTCTGGCTATGCGGTAGAAGCTGCGGATCAGCGTAATCGGCGTATTCCGTTAGTGCGTGTTGAAAGCGGTTTGTATAAAGGTATTTATCAGCACAATTAGTGAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTc                                                                                                                                    >  2:78594/1‑149 (MQ=255)
               tgtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCGGGGTTTcttt                                                                                                                                 <  1:1023497/149‑4 (MQ=255)
               tgtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCGGGGTTTcttt                                                                                                                                 <  1:340471/149‑4 (MQ=255)
                  gCAATATCTTTCTGGCTATGCAGTAGTACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCTTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAGCCCCGGC‑CATGTGCCGGGGTTTctttttg                                                                                                                              <  4:781972/149‑7 (MQ=255)
                   cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                             <  2:1094799/149‑1 (MQ=255)
                    aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCATCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                            <  4:691386/149‑1 (MQ=255)
                      tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                          >  3:1359063/1‑149 (MQ=255)
                      tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                          >  3:832071/1‑149 (MQ=255)
                      tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                        <  1:150533/135‑1 (MQ=255)
                       atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGAc                                                                                                                         >  4:1283080/1‑149 (MQ=255)
                          tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGGCAAAAAAAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACaaa                                                                                                                      >  3:437188/1‑149 (MQ=255)
                                gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                >  4:266397/1‑149 (MQ=255)
                                 cTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                                               >  3:1280140/1‑149 (MQ=255)
                                  tATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                        >  1:163557/1‑123 (MQ=255)
                                   aTGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACGATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCttt                                                                                                             <  3:755520/149‑1 (MQ=255)
                                     gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                           >  2:1144851/1‑149 (MQ=255)
                                        gTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                        >  2:888875/1‑149 (MQ=255)
                                        gTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                        <  3:544743/149‑1 (MQ=255)
                                         tAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTAt                                                                                                       >  3:1291805/1‑149 (MQ=255)
                                           gAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGc                                                                                                     <  3:266397/149‑1 (MQ=255)
                                               cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                           >  2:1114114/1‑139 (MQ=255)
                                               cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                           <  1:1114114/139‑1 (MQ=255)
                                               cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGa                                                                                                 <  1:719255/149‑1 (MQ=255)
                                                     aTCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGc                                                                                                     >  4:72324/1‑139 (MQ=255)
                                                     aTCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGc                                                                                                     <  3:72324/139‑1 (MQ=255)
                                                          cGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                      <  2:1251369/149‑1 (MQ=255)
                                                           gTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTACGCCGCATCCGGCATaaa                                                                                      <  3:626537/148‑1 (MQ=255)
                                                                   cGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                                             <  2:39170/149‑1 (MQ=255)
                                                                    gTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATcc                                                                            <  4:832071/149‑1 (MQ=255)
                                                                    gTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATcc                                                                            >  1:471633/1‑149 (MQ=255)
                                                                    gTATTCCGTTAGTGCG‑TTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                           <  1:235692/117‑1 (MQ=255)
                                                                    gTATTCCGTTAGTGCG‑TTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                           >  2:235692/1‑117 (MQ=255)
                                                                       ttCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAATCCCGGCACATG‑GCCGGGGCCTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCCGATGCGGCGTAAACGCCTTATCAGGc                                                                         <  4:1265428/149‑1 (MQ=255)
                                                                        tCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                        >  1:514105/1‑149 (MQ=255)
                                                                          cGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATcc                                                                            <  1:608560/143‑1 (MQ=255)
                                                                            ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                          >  4:983800/1‑111 (MQ=255)
                                                                            ttAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                          <  3:983800/111‑1 (MQ=255)
                                                                               gTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                       >  4:131925/1‑143 (MQ=255)
                                                                               gTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                       <  3:131925/143‑1 (MQ=255)
                                                                                tGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAt                                                                >  1:330091/1‑149 (MQ=255)
                                                                                 gCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                               <  3:1409887/149‑1 (MQ=255)
                                                                                    ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                            >  4:198970/1‑149 (MQ=255)
                                                                                    ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                            >  2:159277/1‑149 (MQ=255)
                                                                                    ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                            >  3:236766/1‑149 (MQ=255)
                                                                                         aaaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGcc                                                                                                    <  4:503514/104‑1 (MQ=255)
                                                                                         aaaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGcc                                                                                                    >  3:503514/1‑104 (MQ=255)
                                                                                             ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                   <  1:921482/149‑1 (MQ=255)
                                                                                                    tataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                            >  2:666754/1‑149 (MQ=255)
                                                                                                     ataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                                               >  1:827644/1‑81 (MQ=255)
                                                                                                     ataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                           >  2:1143184/1‑149 (MQ=255)
                                                                                                        aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                        >  4:190221/1‑149 (MQ=255)
                                                                                                         aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                       >  2:712905/1‑149 (MQ=255)
                                                                                                         aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                       >  2:72791/1‑149 (MQ=255)
                                                                                                          gATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                      >  3:298308/1‑149 (MQ=255)
                                                                                                           atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                     >  3:95623/1‑149 (MQ=255)
                                                                                                            tatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                    <  4:437188/149‑1 (MQ=255)
                                                                                                             atTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                   <  2:1392168/149‑1 (MQ=255)
                                                                                                              tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                  <  3:190221/149‑1 (MQ=255)
                                                                                                                tATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCAGATGGGGCGTAAACGCCTCATCCGGCCTACAAAATCGTGCAAATTCTACATATTGCAATGCTCTTGt                                >  1:1164770/1‑149 (MQ=37)
                                                                                                                 aTCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                               <  2:1514114/149‑1 (MQ=255)
                                                                                                                  tCAGAACAATTAGTCGAAAAGAAACCCCGGAACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTGTGTTTGTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCGTGGAAATTCAACATATTGCAATTCTCTTGTAg                              >  4:917621/1‑149 (MQ=25)
                                                                                                                  tCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                              >  2:1503747/1‑149 (MQ=255)
                                                                                                                  tCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                              >  1:625079/1‑149 (MQ=255)
                                                                                                                  tCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                              >  4:1453916/1‑149 (MQ=255)
                                                                                                                     gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                   <  1:191859/141‑1 (MQ=255)
                                                                                                                     gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTctct                                   >  2:191859/1‑141 (MQ=255)
                                                                                                                     gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                           >  4:102385/1‑149 (MQ=255)
                                                                                                                     gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGAATATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                           >  2:1033536/1‑149 (MQ=255)
                                                                                                                       aCAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                               <  3:1467230/127‑1 (MQ=255)
                                                                                                                       aCAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                               >  4:1467230/1‑127 (MQ=255)
                                                                                                                         aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                  <  2:175202/138‑1 (MQ=255)
                                                                                                                          aTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATa                      >  1:1264772/1‑149 (MQ=255)
                                                                                                                           ttAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATaa                     >  2:888052/1‑149 (MQ=255)
                                                                                                                              gTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCg                  >  1:973916/1‑149 (MQ=255)
                                                                                                                               tCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGt                 >  1:1414825/1‑149 (MQ=255)
                                                                                                                                cAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTa                <  2:281175/149‑1 (MQ=255)
                                                                                                                                 aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACTTATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg               >  4:11246/1‑149 (MQ=255)
                                                                                                                                 aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg               >  3:1047818/1‑149 (MQ=255)
                                                                                                                                 aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg               >  2:1132575/1‑149 (MQ=255)
                                                                                                                                 aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg               >  1:1001976/1‑149 (MQ=255)
                                                                                                                                  aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                   >  2:884610/1‑112 (MQ=255)
                                                                                                                                  aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                   <  1:884610/112‑1 (MQ=255)
                                                                                                                                  aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc              <  3:198970/149‑1 (MQ=255)
                                                                                                                                  aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc              >  2:1502032/1‑149 (MQ=255)
                                                                                                                                  aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc              <  2:1036162/149‑1 (MQ=255)
                                                                                                                                   aaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg             >  3:1159890/1‑149 (MQ=255)
                                                                                                                                   aaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATGATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg             >  1:1054022/1‑149 (MQ=255)
                                                                                                                                    aaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAACTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcgc            >  4:1180031/1‑149 (MQ=255)
                                                                                                                                        agccccGGC‑CATGTGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATACGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa        >  4:1413724/3‑149 (MQ=18)
                                                                                                                                         aCCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAg       >  3:497156/1‑149 (MQ=255)
                                                                                                                                          ccccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg      <  4:236766/149‑1 (MQ=255)
                                                                                                                                              ggCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATAATGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGa  <  4:1429226/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGA  >  NC_000913/3330270‑3330560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑