Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,677 T→A 6.4% intergenic (+48/+106) yjcE → / ← yjcF putative cation/proton antiporter/pentapeptide repeats protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,6770TA 6.4% 481.9 / 15.1 168intergenic (+48/+106)yjcE/yjcFputative cation/proton antiporter/pentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base T (86/71);  new base A (6/5);  total (92/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 168 total aligned reads displayed.

ACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGCG  >  NC_000913/4281529‑4281825
                                                                                                                                                    |                                                                                                                                                    
accaccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATc                                                                                                                                                     >  2:31206/1‑150 (MQ=255)
  caccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtt                                                                                                                                                   >  2:1456616/1‑150 (MQ=255)
   accTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                                  >  2:1703652/1‑150 (MQ=255)
    ccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                                 >  1:1546881/1‑150 (MQ=255)
    ccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                                 >  2:392576/1‑150 (MQ=255)
    ccTGAGCGCCACACGTGAGATCAGCAACGAAAGGCTTCAAAAATTACTGCACGATCTCGATTTGATTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTTACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                                 <  2:1469393/150‑1 (MQ=255)
     cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                                >  1:942570/1‑150 (MQ=255)
     cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                                >  2:1272848/1‑150 (MQ=255)
          cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCt                                                                                                                                           >  1:1640653/1‑150 (MQ=255)
          cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCt                                                                                                                                           <  1:1193591/150‑1 (MQ=255)
          cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCt                                                                                                                                           <  2:376623/150‑1 (MQ=38)
           gcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                          >  2:648409/1‑150 (MQ=255)
           gcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtt                                                                                                                                          <  1:1266463/150‑1 (MQ=255)
               cacTGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAAc                                                                                                                                      <  1:1721678/150‑1 (MQ=255)
               cacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAAc                                                                                                                                      >  1:1063702/1‑150 (MQ=255)
                acGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATCGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                     >  1:200273/1‑150 (MQ=255)
                   ggAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGt                                                                                                                                  <  1:67481/150‑1 (MQ=255)
                    gagaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCATTTGCACCTTTAACAGGtt                                                                                                                                 <  2:694853/150‑1 (MQ=255)
                       aacAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGAAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGc                                                                                                                              >  1:1152385/3‑150 (MQ=255)
                          aGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCg                                                                                                                           >  2:1513847/1‑150 (MQ=255)
                          aGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCg                                                                                                                           >  2:832611/1‑150 (MQ=255)
                            cAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCggg                                                                                                                         >  2:604513/1‑150 (MQ=255)
                              aCGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCt                                                                                                                       >  2:660766/1‑150 (MQ=255)
                              aCGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTAGCACCTTTAACAGGTTAGCTCGGGCt                                                                                                                       >  1:1113375/1‑150 (MQ=255)
                                    cGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACAtt                                                                                                                 >  1:774621/1‑150 (MQ=255)
                                      cTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCa                                                                                                               <  1:415162/150‑1 (MQ=255)
                                            aaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTc                                                                                                         <  2:315811/150‑1 (MQ=255)
                                             aaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCt                                                                                                        <  2:893899/150‑1 (MQ=255)
                                               ttACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaa                                                                                                      <  1:979035/150‑1 (MQ=255)
                                                  cTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAAc                                                                                                   >  2:1020217/1‑150 (MQ=255)
                                                  cTGCACGATCTAGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGATCGTGCTAACATGCATGAGTCTAAAAc                                                                                                   >  2:1013920/1‑150 (MQ=255)
                                                   tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACa                                                                                                  >  1:492705/1‑150 (MQ=255)
                                                    gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                 >  2:769370/1‑150 (MQ=255)
                                                     cACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaa                                                                                                <  1:637934/150‑1 (MQ=255)
                                                           ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTACGCCGAAAGCGACGTAACCAGGCCCGAGTTAATCTGAAAAATGTGCAATCTTTTGCACGTTTAACAGGTTAGCTCGGGCTTACATTCATGATTCTAAAACAAATAAAtt                                                                                          <  1:6799/150‑1 (MQ=255)
                                                           ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                          <  2:1519773/150‑1 (MQ=255)
                                                           ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                          <  2:1337096/150‑1 (MQ=255)
                                                               aTTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtct                                                                                      >  1:1655070/1‑150 (MQ=255)
                                                                  tGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtctctc                                                                                   <  2:579300/150‑1 (MQ=255)
                                                                    cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtt                                                                                 >  2:354984/1‑150 (MQ=255)
                                                                     ttGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCttt                                                                                <  1:912948/150‑1 (MQ=255)
                                                                      ttaaGCGTTGCTAATTTAGGAAAAACAGTAAGCCTAAAGCTACGTAACCAGGCCCGAGTTAATATGATAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                               <  1:855572/148‑1 (MQ=255)
                                                                          gCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAg                                                                           <  1:1413153/150‑1 (MQ=255)
                                                                            gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                         >  1:893191/1‑150 (MQ=255)
                                                                            gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                         <  1:1703652/150‑1 (MQ=255)
                                                                                cTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCATTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAAtt                                                                     >  2:611456/1‑150 (MQ=255)
                                                                                    ttGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAtt                                                                 <  2:942570/150‑1 (MQ=255)
                                                                                      gAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTgg                                                               <  2:1261547/150‑1 (MQ=255)
                                                                                        ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTTCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                             >  1:1452142/1‑150 (MQ=255)
                                                                                        ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                             >  2:1323555/1‑150 (MQ=255)
                                                                                        ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                             >  1:581813/1‑150 (MQ=255)
                                                                                           aaaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGtt                                                          >  2:222455/1‑150 (MQ=255)
                                                                                            aaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttc                                                         <  1:1323555/150‑1 (MQ=255)
                                                                                               cactaccccGAACGCGACGTACCCAGGCCCGAGGTAATCTGAAAAAGCTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                      <  2:1654712/143‑1 (MQ=255)
                                                                                                 gTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGa                                                    >  1:1107817/1‑150 (MQ=255)
                                                                                                  tAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAc                                                   <  2:760698/150‑1 (MQ=255)
                                                                                                   aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                                  <  1:765663/150‑1 (MQ=255)
                                                                                                   aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAACTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                                  <  2:1623711/150‑1 (MQ=255)
                                                                                                   aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAACTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                                  <  1:1135871/150‑1 (MQ=255)
                                                                                                      ccGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAAGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTaa                                               >  1:1206218/1‑150 (MQ=255)
                                                                                                        gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAcc                                             <  2:1063702/150‑1 (MQ=255)
                                                                                                        gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAcc                                             <  1:1355632/150‑1 (MQ=255)
                                                                                                          aaGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTc                                           >  2:1548851/1‑150 (MQ=255)
                                                                                                            gCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCNAACANTCATGATTNTAAAACAANTAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCaa                                         <  1:1456616/150‑1 (MQ=255)
                                                                                                              gACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATc                                       <  2:502418/150‑1 (MQ=255)
                                                                                                              gACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATc                                       <  1:1272848/150‑1 (MQ=255)
                                                                                                                cGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTg                                     <  1:1349975/150‑1 (MQ=255)
                                                                                                                cGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTg                                     >  2:458592/1‑150 (MQ=255)
                                                                                                                 gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGc                                    >  1:1024670/1‑150 (MQ=255)
                                                                                                                  tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                                   <  2:1175639/150‑1 (MQ=255)
                                                                                                                  tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                                   <  1:832611/150‑1 (MQ=255)
                                                                                                                   aaCCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCAt                                  >  2:1573091/1‑150 (MQ=255)
                                                                                                                    aCCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATc                                 <  2:1429227/150‑1 (MQ=255)
                                                                                                                     ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTATCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAATTAACCTCAATCTGCAtat                                >  2:852924/1‑148 (MQ=255)
                                                                                                                      cAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCtt                               >  1:660209/1‑150 (MQ=255)
                                                                                                                        ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                             >  2:1657149/1‑150 (MQ=255)
                                                                                                                        ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                             >  2:844368/1‑150 (MQ=255)
                                                                                                                          cccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCat                           >  1:251611/1‑150 (MQ=255)
                                                                                                                             gAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTATTCGACTTAACCTCAATCTTCATCTTGCataat                        >  2:870195/1‑150 (MQ=255)
                                                                                                                                 tAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCt                    <  2:568795/150‑1 (MQ=255)
                                                                                                                                  aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTg                   >  2:1508820/1‑150 (MQ=255)
                                                                                                                                  aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTg                   >  2:1184968/1‑150 (MQ=255)
                                                                                                                                   aTCTGAAAAAGGTGCTATCTTTTGCAACTTTAACTGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAATATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                  <  1:270710/150‑1 (MQ=255)
                                                                                                                                   aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                  <  2:1572631/150‑1 (MQ=255)
                                                                                                                                   aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                  <  1:1573091/150‑1 (MQ=255)
                                                                                                                                   aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                  >  2:911551/1‑150 (MQ=255)
                                                                                                                                     atGAAAAAGGTGCAATCTTTTTCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAg                <  1:627411/149‑1 (MQ=255)
                                                                                                                                       gttaaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAAAAAATAAATTATCTATCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTa              >  1:702294/4‑150 (MQ=255)
                                                                                                                                       gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTa              >  1:494283/1‑150 (MQ=255)
                                                                                                                                         taaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGCCTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt            >  1:1627043/2‑150 (MQ=255)
                                                                                                                                         aaaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt            >  1:1122692/1‑150 (MQ=255)
                                                                                                                                          aaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATc           >  1:1693338/1‑150 (MQ=255)
                                                                                                                                          aaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATc           >  2:1170513/1‑150 (MQ=255)
                                                                                                                                           aaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCt          >  1:843707/1‑150 (MQ=255)
                                                                                                                                            aGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCtt         >  2:986183/1‑150 (MQ=255)
                                                                                                                                             ggTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTc        <  2:1392650/150‑1 (MQ=255)
                                                                                                                                              gTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCa       >  1:409707/1‑150 (MQ=255)
                                                                                                                                               tGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCAt      <  2:774621/150‑1 (MQ=255)
                                                                                                                                               tGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCAt      <  1:458592/150‑1 (MQ=255)
                                                                                                                                                gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATa     >  1:793630/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGCG  >  NC_000913/4281529‑4281825

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