Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,939 T→C 8.7% intergenic (+17/‑150) acrZ → / → modA AcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913794,9390TC 8.7% 178.3 / 6.2 71intergenic (+17/‑150)acrZ/modAAcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base C (2/4);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAC  >  NC_000913/794795‑795086
                                                                                                                                                |                                                                                                                                                   
tGGTATTCCCCGTAATCATGGTACCTTTCGTGATGGCCATCATAATGGGGCAGATTTACGGTCTTGGTGAAGAATTCTCCATCTGATCTGCTGTTGGTAACAAAGCCCAGCGCGGACAAACGCATTGATTCCCTGAATGCGCGCCTAgg                                                                                                                                                 >  2:295695/1‑148 (MQ=255)
   tATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCgg                                                                                                                                              <  2:501282/149‑1 (MQ=255)
     ttCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTTGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGc                                                                                                                                            <  1:39960/149‑1 (MQ=255)
       cGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                          <  1:526817/149‑1 (MQ=255)
       cGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGAATTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                          <  1:532678/149‑1 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtt                                                                                                                                         >  1:979669/1‑149 (MQ=255)
          cGTAATCATGGTACCTTTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttca                                                                                                                                       <  1:116044/149‑3 (MQ=255)
           gTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcag                                                                                                                                      <  2:661513/149‑4 (MQ=255)
           gTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGCAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcag                                                                                                                                      <  1:42774/149‑4 (MQ=255)
             aaTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACCTCTTTTCTGGTGTTGGTAAAAAAGCCCAGCCCGGCCAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctt                                                                                                                                    >  2:797424/1‑149 (MQ=255)
             aaTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctt                                                                                                                                    <  1:418031/149‑1 (MQ=255)
              aTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttt                                                                                                                                   <  2:305795/149‑1 (MQ=255)
              aTCATGGTACCTGTCGGGATGGCCATCAACCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttt                                                                                                                                   <  2:457856/149‑1 (MQ=255)
                 aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                <  2:719675/149‑1 (MQ=255)
                   ggTACCTGTCGTGATGGCCATCATCCTGGGTCTGAATTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTtctc                                                                                                                              <  2:387532/149‑1 (MQ=255)
                       ccTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACt                                                                                                                          >  1:600017/1‑149 (MQ=255)
                        cTGTCGTGATGGCCATAATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTGTCAACtt                                                                                                                         >  2:984163/1‑149 (MQ=255)
                              tGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGCTTCACTGAATGCCCGCTTAGTCGGGCATTTTCTTTTACTCAACTTCCTGCt                                                                                                                   >  1:87740/1‑149 (MQ=255)
                               gATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCtt                                                                                                                  >  2:332722/1‑149 (MQ=255)
                                  ggCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTc                                                                                                               <  1:181330/149‑1 (MQ=255)
                                    ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCt                                                                                                             >  1:490179/1‑149 (MQ=255)
                                          tcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACCTTCCTGATTTTCCTGCCGat                                                                                                       >  2:439415/1‑149 (MQ=255)
                                                 tcggATGTACGGTCGAGGTGCAGTATTCACCCTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACACAATCATTCATTCCCTGACTGCCCGCTTAGTCGGGCATTTTCTATTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                <  1:702287/146‑1 (MQ=255)
                                                 tCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                >  1:546691/1‑149 (MQ=255)
                                                 tCTGATTTACGGTCTTGGGCTAGTATTCAACTTCTTTGCTGGTGTTGGTAAAAAAGACCAGCCGGGACAAACTCCTGGATTTCCTGAATGCCCGCTAAGTGGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                <  2:731622/149‑1 (MQ=255)
                                                 gcTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                <  2:735251/148‑1 (MQ=255)
                                                  ttGATTTACGGTCTTTGTGATGTATACAACATCTTTTCTGGTGTTTGTAAAAAAGACCAGCACGGACAAATTGATTGATTCCCTGAATGCCGGCCTATTCGTGCATTTTCTTTTTCTGAACTTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                               <  2:245704/148‑1 (MQ=255)
                                                     aTTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGCCCAGCCCGGACAAAATCATTGCTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                            >  2:652141/1‑149 (MQ=255)
                                                         aagttCTTGGGGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGTACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTa                                                                                        <  2:354834/145‑1 (MQ=255)
                                                           ggTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGTGCATTTTCTATTTCTAAACTTCATGCTTTTCATGCCGATATTTTTTCTTATCTAcc                                                                                      >  1:751278/1‑149 (MQ=255)
                                                             tCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAACGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCGTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc                                                                                    >  1:502614/1‑149 (MQ=255)
                                                               ttGGTGCAGTATTCAACATCTTTTCTGTTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCGTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTcac                                                                                  <  2:789551/149‑1 (MQ=255)
                                                                     aaGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                            <  1:291344/149‑1 (MQ=255)
                                                                        ttttCATCATCTTTTCTGGTGTTGGTAAAAAAGAACAGCCCGGACAAAATCATTGATTCCCTGAATGCCCTCTTAGTCGTGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg                                                                         <  2:323451/147‑1 (MQ=255)
                                                                        tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCATCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg                                                                         <  2:319664/149‑1 (MQ=255)
                                                                        tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACAAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCGAGATCACAAAAGTTat                                                                         >  1:695949/1‑148 (MQ=255)
                                                                            cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                     >  1:168081/1‑149 (MQ=255)
                                                                               aaTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAc                                                                  <  1:528992/148‑1 (MQ=255)
                                                                                aTCTTTTCTGGTGTTGGTAAAAAAGCCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGAATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCTATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATTAct                                                                 >  1:734767/1‑149 (MQ=255)
                                                                                aTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCATTTAct                                                                 >  1:774250/1‑149 (MQ=255)
                                                                                aTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                 >  2:375403/1‑149 (MQ=255)
                                                                                  cTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgc                                                               <  2:928497/149‑1 (MQ=255)
                                                                                   ttttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct                                                              <  2:600017/149‑1 (MQ=255)
                                                                                   ttttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCAATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct                                                              >  1:885803/1‑149 (MQ=255)
                                                                                              tGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTcc                                                   <  1:81577/149‑1 (MQ=255)
                                                                                               ggaaaaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                  <  1:329614/146‑1 (MQ=255)
                                                                                                 tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCCTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATATACCTCACAAACGTTAGCAATAACTGCTTGGAAAATTACGat                                                >  1:166229/1‑148 (MQ=255)
                                                                                                 tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTAAGTCGGGCATTTTCTTTTTGTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGTAAAAttaggag                                                >  1:697063/1‑144 (MQ=255)
                                                                                                  aaaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCGTAGTCGGGCATTTTCTTTTTATCAAGTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACAGCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGt                                               >  2:684860/1‑149 (MQ=255)
                                                                                                     aaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAg                                            <  2:490179/149‑1 (MQ=255)
                                                                                                                cGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCGTTTTCTCAACTTCCTTCTTTTCCTGCCGATATTTTTTCTTATCTACCACACAAAGGTAAGCAATAACTGGTGGGAATATTCCGAGTTAGTCTTTTTATTg                                 >  2:350739/1‑149 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGGCCGCTTAGTCGGGCATTTTCTTTTTCTCTACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                               <  1:546852/149‑1 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                               <  1:643340/149‑1 (MQ=255)
                                                                                                                    cAAACTCATTGATTGACTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGc                             <  2:119067/149‑1 (MQ=255)
                                                                                                                     aaaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcg                            >  1:649755/1‑148 (MQ=255)
                                                                                                                     aaaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTAACTCACAAAGTTTAGCAATAACTTCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                            >  2:932523/1‑149 (MQ=255)
                                                                                                                        aTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTAc                         <  2:546691/149‑1 (MQ=255)
                                                                                                                             tGATTCCCTTAATGCCCGCTTTGGCGGGCATTTTCTTTTTCCCATCTTCCTGCTTTTCCTGCCTATATTTTTTCTTATCTAACTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAAc                    <  2:572201/149‑1 (MQ=255)
                                                                                                                              gATTCCCTGATTGCCCGCTTAGTCGGGCATTTTCTTTTTCTGAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGGTAGTCGTTATATTGTCGCCTACATAACg                   <  1:888691/149‑1 (MQ=255)
                                                                                                                                ttCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCATCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGtt                 <  1:1428/149‑1 (MQ=255)
                                                                                                                                 tCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTa                <  1:465614/149‑1 (MQ=255)
                                                                                                                                  tccTGAATTCCCGCATAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATATACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc               <  2:751278/148‑1 (MQ=255)
                                                                                                                                    gaaaaTGCCCGACTAAGCGGGCATTTTCATTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGAAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt             >  2:665380/4‑149 (MQ=255)
                                                                                                                                    cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt             >  1:31298/1‑149 (MQ=255)
                                                                                                                                        aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAg         <  1:261300/149‑1 (MQ=255)
                                                                                                                                        aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTAAACTTCCTGCTTTTCCTGCCGATATATTTTCTTATCTACCTCACACAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAg         <  2:942664/149‑1 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATATCGTTACATTAAGGGGt     >  1:608687/1‑149 (MQ=255)
                                                                                                                                             cGCTTTGGCGGGCATTTTCTTTTTCTCAACTTCGTGCTTTTCCTGCCGTTATTTTTTCTTATCTACCTCTCAAAGGTTTTCAATAACTGCGGGGATAATTACGAGTTAGTGGTTATATAGTCGCGTAAATATCGTTAATATAAGCGGtt    >  2:211286/1‑149 (MQ=255)
                                                                                                                                             cGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt    >  1:879706/1‑149 (MQ=255)
                                                                                                                                              gCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTa   <  1:979411/149‑1 (MQ=255)
                                                                                                                                               cTTAGTCGGGCATTTTCTTTTTCTCACCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCTCCTACATAAGTGTACATTAAGGGGGTAc  >  2:674008/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAC  >  NC_000913/794795‑795086

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