Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,677 T→A 5.4% intergenic (+48/+106) yjcE → / ← yjcF putative transporter YjcE/pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,6770TA 5.4% 868.7 / 22.3 262intergenic (+48/+106)yjcE/yjcFputative transporter YjcE/pentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base T (128/119);  new base A (10/5);  total (138/124)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

CCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGC  >  NC_000913/4281530‑4281824
                                                                                                                                                   |                                                                                                                                                   
    cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                               <  2:786715/150‑1 (MQ=255)
     tGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGc                                                                                                                                              <  1:887374/150‑1 (MQ=255)
     tGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGc                                                                                                                                              >  1:1204163/1‑150 (MQ=255)
           ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                        >  1:1374965/1‑150 (MQ=255)
              cacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAAc                                                                                                                                     >  2:486510/1‑150 (MQ=255)
               acaggAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                    >  1:7150/4‑150 (MQ=255)
               acGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                    <  2:1115834/150‑1 (MQ=255)
                  ggAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTTAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtttttcagat                                                                                                                                 <  1:92533/150‑8 (MQ=39)
                  ggAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtttttcagat                                                                                                                                 <  1:858967/150‑8 (MQ=39)
                   gagaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGtt                                                                                                                                <  2:557394/150‑1 (MQ=255)
                   gagaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGtt                                                                                                                                <  1:806876/150‑1 (MQ=255)
                    agaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTa                                                                                                                               <  1:947630/150‑1 (MQ=255)
                     gaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAg                                                                                                                              <  1:486510/150‑1 (MQ=255)
                      aTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGc                                                                                                                             <  1:337747/150‑1 (MQ=255)
                        cAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTc                                                                                                                           <  2:213652/150‑1 (MQ=255)
                        cAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTc                                                                                                                           >  1:112641/1‑150 (MQ=255)
                            aaCGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGc                                                                                                                       <  1:515202/150‑1 (MQ=255)
                                     cTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCa                                                                                                              >  1:1252798/1‑150 (MQ=255)
                                      tGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                             <  1:228306/150‑1 (MQ=255)
                                      tGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                             <  1:612062/150‑1 (MQ=255)
                                      tGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                             >  1:983208/1‑150 (MQ=255)
                                       gCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATg                                                                                                            >  2:504687/1‑150 (MQ=255)
                                        cAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGa                                                                                                           >  1:577126/1‑150 (MQ=255)
                                        cAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGa                                                                                                           >  1:945627/1‑150 (MQ=255)
                                           aaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGAttg                                                                                                        <  1:994170/150‑2 (MQ=255)
                                           aaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTc                                                                                                        <  1:999087/150‑1 (MQ=255)
                                           aaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTc                                                                                                        <  1:1389588/150‑1 (MQ=255)
                                               tACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaa                                                                                                    <  2:1407382/150‑1 (MQ=255)
                                                aCTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGATGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaaa                                                                                                   <  2:746327/150‑1 (MQ=255)
                                                 cTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAAc                                                                                                  <  2:640768/150‑1 (MQ=255)
                                                  tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACa                                                                                                 >  2:924663/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  2:892500/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  2:1345362/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  1:793569/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  2:81395/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  1:1011700/1‑150 (MQ=255)
                                                   gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                >  1:30001/1‑150 (MQ=255)
                                                             gATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtc                                                                                      >  1:646132/1‑150 (MQ=255)
                                                              aTTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtct                                                                                     <  2:1217020/150‑1 (MQ=255)
                                                              aTTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtct                                                                                     >  1:55991/1‑150 (MQ=255)
                                                                    ttGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCttt                                                                               <  2:705734/150‑1 (MQ=255)
                                                                     tGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                              >  1:166406/1‑150 (MQ=255)
                                                                      gAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTa                                                                             >  2:1148414/1‑150 (MQ=255)
                                                                         gCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAg                                                                          <  2:735380/150‑1 (MQ=255)
                                                                           gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                        >  2:326359/1‑150 (MQ=255)
                                                                           gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                        <  2:240587/150‑1 (MQ=255)
                                                                           gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                        <  1:53261/150‑1 (MQ=255)
                                                                             tGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                        <  2:803078/148‑1 (MQ=255)
                                                                             tGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                                                      >  1:184476/1‑150 (MQ=255)
                                                                             tGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                                                      >  1:1198436/1‑150 (MQ=255)
                                                                              gCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAAt                                                                     >  2:1306268/1‑150 (MQ=255)
                                                                                 aaTTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCa                                                                  >  2:1375326/1‑150 (MQ=255)
                                                                                 aaTTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCa                                                                  >  2:33397/1‑150 (MQ=255)
                                                                                   ttGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAtt                                                                >  2:296388/1‑150 (MQ=255)
                                                                                   ttGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAtt                                                                >  2:618693/1‑150 (MQ=255)
                                                                                      aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                             >  1:1281491/1‑150 (MQ=255)
                                                                                      aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                             <  1:1239701/150‑1 (MQ=255)
                                                                                       ggAAAATCAGTAAGCCGAAAGCTACGTAACCAGGCCCGAGGTAATATGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            >  2:1045433/1‑150 (MQ=255)
                                                                                       ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            >  1:1087249/1‑150 (MQ=255)
                                                                                       ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            <  2:1024380/150‑1 (MQ=255)
                                                                                        gAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTg                                                           <  2:166406/150‑1 (MQ=255)
                                                                                        gAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTGTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTg                                                           <  1:849061/150‑1 (MQ=255)
                                                                                           aaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGTGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTTGCTGttc                                                        <  1:1150738/150‑1 (MQ=255)
                                                                                           aaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttc                                                        <  1:1417919/150‑1 (MQ=255)
                                                                                              cAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                     >  1:1118365/1‑150 (MQ=255)
                                                                                               aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                                    <  2:17957/150‑1 (MQ=255)
                                                                                               aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                                    >  2:1246504/1‑150 (MQ=255)
                                                                                                gTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGa                                                   >  2:1161228/1‑150 (MQ=255)
                                                                                                   aGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACtt                                                >  2:651063/1‑150 (MQ=255)
                                                                                                   aGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACtt                                                >  1:375104/1‑150 (MQ=255)
                                                                                                    gCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTa                                               >  2:426949/1‑150 (MQ=255)
                                                                                                    gCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTa                                               >  1:653305/1‑150 (MQ=255)
                                                                                                    gCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTa                                               <  2:771626/150‑1 (MQ=255)
                                                                                                           gCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCaa                                        >  1:70635/1‑150 (MQ=255)
                                                                                                           gCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCaa                                        <  1:878221/150‑1 (MQ=255)
                                                                                                            cGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAAt                                       >  2:566775/1‑150 (MQ=255)
                                                                                                              aCGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCt                                     <  1:491255/150‑1 (MQ=255)
                                                                                                                gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGc                                   >  2:763/1‑150 (MQ=255)
                                                                                                                    ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGCCTTAACCTCAATCTGCATCt                               >  2:1138444/1‑150 (MQ=255)
                                                                                                                     ttagCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCtt                              >  1:1364736/4‑150 (MQ=21)
                                                                                                                      aGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTg                             <  2:577126/150‑1 (MQ=255)
                                                                                                                         cccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCat                          >  2:1432583/1‑150 (MQ=255)
                                                                                                                          ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTGTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                         >  1:462927/1‑150 (MQ=255)
                                                                                                                              ggTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAatat                     >  2:280560/1‑150 (MQ=255)
                                                                                                                                tAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCt                   >  1:397385/1‑150 (MQ=255)
                                                                                                                                 aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                         >  1:768563/1‑143 (MQ=255)
                                                                                                                                  aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                 <  1:1157427/150‑1 (MQ=255)
                                                                                                                                   tCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCa                >  1:6345/1‑150 (MQ=255)
                                                                                                                                    cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                         <  1:779216/140‑1 (MQ=255)
                                                                                                                                    cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                         >  2:779216/1‑140 (MQ=255)
                                                                                                                                     tGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt              <  2:1204163/150‑1 (MQ=255)
                                                                                                                                      gttaaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAACTCAATCTGCATCTTGCATAATATCTGCAGTa             >  2:929781/4‑150 (MQ=255)
                                                                                                                                      gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTa             >  1:32847/1‑150 (MQ=255)
                                                                                                                                        aaaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt           >  1:1329220/1‑150 (MQ=255)
                                                                                                                                         aaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATc          <  2:154236/150‑1 (MQ=255)
                                                                                                                                          aaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCt         >  2:442962/1‑150 (MQ=255)
                                                                                                                                               gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATa    >  2:954064/1‑150 (MQ=255)
                                                                                                                                               gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATa    >  2:52634/1‑150 (MQ=255)
                                                                                                                                                cAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAg   <  2:710493/150‑1 (MQ=255)
                                                                                                                                                 aaTCTTTTGCNCCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAgg  >  1:1053086/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGC  >  NC_000913/4281530‑4281824

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