Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,419:1 +T 5.6% intergenic (+23/+163) 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,4191.T 5.6% 388.3 / 16.7 90intergenic (+23/+163)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base . (40/45);  new base T (3/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  NC_000913/3330273‑3330566
                                                                                                                                                   |                                                                                                                                                   
aCGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑Gc                                                                                                                                                   >  3:776155/1‑149 (MQ=255)
    aTGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCggg                                                                                                                                               <  3:429294/149‑1 (MQ=255)
    aTGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCggg                                                                                                                                               <  3:840769/149‑1 (MQ=255)
     tGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                              <  2:15089/149‑1 (MQ=255)
      ggCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGc                                                                                                                                             <  2:726150/149‑1 (MQ=255)
       gCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCt                                                                                                                                            <  2:274493/149‑1 (MQ=255)
       gCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCt                                                                                                                                            >  4:367760/1‑149 (MQ=255)
            tgtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGa                                                                                                                                       <  2:738713/149‑1 (MQ=255)
             gtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAGCCCCGGC‑CATGTGCCGGGGTTTctttt                                                                                                                                      <  2:28674/149‑5 (MQ=255)
                cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                                   <  1:361618/149‑1 (MQ=255)
                 aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAGCCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                  <  4:432100/149‑1 (MQ=255)
                 aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                  <  4:641597/149‑1 (MQ=255)
                   tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGTTTCAGATTATTGa                                                                                                                                >  4:213313/1‑149 (MQ=255)
                   tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                                <  3:822348/149‑1 (MQ=255)
                    atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGAc                                                                                                                               >  2:678777/1‑149 (MQ=255)
                    atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGAc                                                                                                                               >  3:158305/1‑149 (MQ=255)
                     tCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACa                                                                                                                              <  1:382006/149‑1 (MQ=255)
                       tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACaaa                                                                                                                            >  3:48999/1‑149 (MQ=255)
                         tCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGt                                                                                                                          <  4:594912/149‑1 (MQ=255)
                                  gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                 <  3:495723/149‑1 (MQ=255)
                                    aGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGttt                                                                                                               >  3:431303/1‑149 (MQ=255)
                                      tAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTAt                                                                                                             <  1:660977/149‑1 (MQ=255)
                                       aGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATg                                                                                                            <  1:275946/149‑1 (MQ=255)
                                        gAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGc                                                                                                           <  4:501082/149‑1 (MQ=255)
                                          aCCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCg                                                                                                         >  4:272838/1‑149 (MQ=255)
                                            cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGa                                                                                                       >  4:447061/1‑149 (MQ=255)
                                              gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                     >  1:574992/1‑149 (MQ=255)
                                                 gATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                  >  1:839380/1‑149 (MQ=255)
                                                  aTCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                 >  1:192006/1‑149 (MQ=255)
                                                    cAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                               <  2:568256/149‑1 (MQ=255)
                                                      gCGTAATCGACGTATTCCGTTAGTGCGCTTTGAAAGCCGCTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTAGTTTACGCCGCATCCGGCATaa                                                                                             <  2:301138/149‑1 (MQ=255)
                                                       cGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTACGCCGCATCCGGCATaaa                                                                                            <  1:678777/149‑1 (MQ=255)
                                                             cGCCGTATTCCGTTGGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                                      >  3:68153/1‑149 (MQ=255)
                                                             cGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                                      >  3:522949/1‑149 (MQ=255)
                                                               ccGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCTGGATGCGGCGTAAACGCCTTAt                                                                                    >  4:33454/1‑149 (MQ=255)
                                                               ccGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                    >  3:868670/1‑149 (MQ=255)
                                                               ccGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                    >  4:655187/1‑149 (MQ=255)
                                                                    ttCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                               >  2:271671/1‑149 (MQ=255)
                                                                    ttCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                               <  4:362887/149‑1 (MQ=255)
                                                                      acGTTAGTGCGATTCAAAAACCGTTTGTATAAAGACATTTATCAGAAAGATTAGTCAAAAAGAAACCCCAGTATATG‑GCCGAGGCGTCAGATTATTAACAAAGCGCGCTTTGTTTATACCGGATGCGGCGTAAACGCCCTATCCGGccg                                                                             <  4:357600/148‑2 (MQ=25)
                                                                       cGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                            >  3:147566/1‑149 (MQ=255)
                                                                              gCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                     <  2:787141/149‑1 (MQ=255)
                                                                                gTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                   <  1:833814/149‑1 (MQ=255)
                                                                                 ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                  >  3:727609/1‑149 (MQ=255)
                                                                                 ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                  <  4:648741/149‑1 (MQ=255)
                                                                                 ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                  <  3:380602/149‑1 (MQ=255)
                                                                                       aaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                            <  1:791360/149‑1 (MQ=255)
                                                                                        aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                           <  4:798277/149‑1 (MQ=255)
                                                                                          ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                         <  1:19198/149‑1 (MQ=255)
                                                                                          ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                         <  4:862587/149‑1 (MQ=255)
                                                                                           cGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAAACCCGGCACATG‑GCCGGGGCTTCAGATTATTAACAAAGTGCGCTTTACTTATCACGGATGCGGCGTAAACGTCTTATCTGTCCTAAAAAATAGTGCAAATCCAAc                                                        <  4:533015/149‑1 (MQ=255)
                                                                                                 tataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                  >  1:432520/1‑149 (MQ=255)
                                                                                                  ataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGATAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATAtttc                                                 <  4:634039/149‑3 (MQ=255)
                                                                                                     aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                              >  1:840409/1‑149 (MQ=255)
                                                                                                      aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                             >  1:117466/1‑149 (MQ=255)
                                                                                                        atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                                <  4:777005/64‑1 (MQ=255)
                                                                                                        atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                                >  3:777005/1‑64 (MQ=255)
                                                                                                        atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                       <  1:551196/105‑1 (MQ=255)
                                                                                                        atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                       >  2:551196/1‑105 (MQ=255)
                                                                                                        atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                           <  4:396740/149‑1 (MQ=255)
                                                                                                         tatTTATCAGAATAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                          <  2:25390/149‑1 (MQ=255)
                                                                                                           tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                        >  3:509125/1‑149 (MQ=255)
                                                                                                           tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                        <  4:48999/149‑1 (MQ=255)
                                                                                                              aTCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCTGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                     <  3:211853/149‑1 (MQ=255)
                                                                                                               tCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                                    <  4:68153/149‑1 (MQ=255)
                                                                                                               tCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                                    >  1:748405/1‑149 (MQ=255)
                                                                                                                 aGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                  >  1:810804/1‑149 (MQ=255)
                                                                                                                 aGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                  >  2:612079/1‑149 (MQ=255)
                                                                                                                  gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                                    >  2:673283/1‑146 (MQ=255)
                                                                                                                  gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAg                                    <  1:673283/146‑1 (MQ=255)
                                                                                                                   aaCAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                <  2:107006/149‑1 (MQ=255)
                                                                                                                   aaCAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                >  1:723733/1‑149 (MQ=255)
                                                                                                                     cAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                              <  4:492847/149‑1 (MQ=255)
                                                                                                                      aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                             >  2:507705/1‑149 (MQ=255)
                                                                                                                              aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  4:470526/1‑149 (MQ=255)
                                                                                                                              aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  2:595068/1‑149 (MQ=255)
                                                                                                                               aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                    <  4:868670/149‑1 (MQ=255)
                                                                                                                               aaaaGAAACCCCGGCACATG‑ACCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                    >  3:344937/1‑149 (MQ=255)
                                                                                                                                tctgaagccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                   >  1:226064/8‑149 (MQ=21)
                                                                                                                                tctgAAGCCCCGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                   >  2:279521/4‑149 (MQ=18)
                                                                                                                                 aaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcgc                  >  1:750389/1‑149 (MQ=255)
                                                                                                                                    aagccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACACAATCGTGCAAATTCAACATATTGCAATTCGCTTGTAGGCCTGATAAGCGTAGCGCATc               >  4:759091/4‑149 (MQ=18)
                                                                                                                                       ccccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  4:38026/1‑149 (MQ=255)
                                                                                                                                        cccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:435737/1‑149 (MQ=255)
                                                                                                                                        cccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc           <  2:840409/149‑1 (MQ=255)
                                                                                                                                           ggCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGa        <  1:612079/149‑1 (MQ=255)
                                                                                                                                             cacaTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt      <  1:410260/149‑1 (MQ=255)
                                                                                                                                             cacaTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt      <  4:147566/149‑1 (MQ=255)
                                                                                                                                                 tG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc  >  4:627095/1‑149 (MQ=255)
                                                                                                                                                 tG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc  >  2:260/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  NC_000913/3330273‑3330566

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