Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,678 C→A 5.8% intergenic (+49/+105) 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,6780CA 5.8% 465.4 / 8.5 152intergenic (+49/+105)yjcE/yjcFputative cation/proton antiporter/pentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base C (70/72);  new base A (6/4);  total (76/76)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

ACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGCGAC  >  NC_000913/4281529‑4281827
                                                                                                                                                      |                                                                                                                                                     
 ccaccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCt                                                                                                                                                      <  2:1497650/150‑1 (MQ=255)
   cccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTTCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                                    <  1:857280/149‑1 (MQ=255)
      tGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAATTCAGTAAGCCTAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGc                                                                                                                                                 <  2:1148760/150‑1 (MQ=255)
         gcgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAcc                                                                                                                                              <  1:375259/150‑1 (MQ=255)
          cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGAAGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCt                                                                                                                                             >  1:122055/1‑150 (MQ=255)
           gcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                            >  2:192019/1‑150 (MQ=255)
            ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCCCGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                           >  2:1606236/1‑150 (MQ=255)
            ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGAACGAGCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACAAGGCCCGAGGTAATCTGAAAACGGTGCAATCTTTTGCACCttt                                                                                                                                           >  2:883384/1‑150 (MQ=255)
            ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACATCCCGATCTCGATTTGCTTGAAGCGGGGCTAAGAGAGGACAATCAGTAACCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGGGCAATCTTTTGCACCatt                                                                                                                                           >  1:328061/1‑150 (MQ=255)
                acGTGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                       <  1:1594570/150‑1 (MQ=255)
                acGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCGCGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                       <  1:686473/150‑1 (MQ=255)
                acGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACa                                                                                                                                       >  2:788097/1‑150 (MQ=255)
                      taTCAGCAACGAAACGCTGCAAAAATTACTTCACTATCTCGATTTCCTTGAAGCGTTGCTAATTGACGAAAATCACTAAGCCGAAAGCCACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAg                                                                                                                                 <  1:1653513/149‑1 (MQ=255)
                         cAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTc                                                                                                                              >  1:789066/1‑150 (MQ=255)
                           gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                            >  1:760404/1‑150 (MQ=255)
                           gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                            >  1:579970/1‑150 (MQ=255)
                             aaCGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGc                                                                                                                          <  2:568007/150‑1 (MQ=21)
                             aaCGAAACGCTGCAAAAATTACTGCACGAGCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGc                                                                                                                          <  1:208600/150‑1 (MQ=255)
                                 aaaCGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAAc                                                                                                                      <  1:1489214/150‑1 (MQ=255)
                                  aaCGCTGCAAAAATTACTGCACGATGTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACa                                                                                                                     <  2:1506394/150‑1 (MQ=255)
                                  aaCGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACa                                                                                                                     <  1:1567103/150‑1 (MQ=255)
                                      cTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAGAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTTAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCCCGAGCTAAcctgtt                                                                                                                 <  2:1219340/150‑6 (MQ=255)
                                      cTGAAAAAATTAATGCACGAACTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGTCGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCa                                                                                                                 <  1:1645875/150‑1 (MQ=255)
                                       tGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCa                                                                                                                 <  2:1039187/150‑1 (MQ=255)
                                       tGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGATAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                                <  1:1288097/150‑1 (MQ=255)
                                        gCAAAAATTACTGCACTATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGATGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATg                                                                                                               <  1:1111994/150‑1 (MQ=255)
                                        gCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGTAAATCAGTATGCCGAAAGCGACGTAACCAGGCCAGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATg                                                                                                               <  2:878025/150‑1 (MQ=255)
                                         cAAAAATTACTGCACGATCTCGATTTGCTTGACGCGTTGCTAATTGAGGGAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCCCCTTTAACAGGTTAGCTCGGGCTAACATTCATGa                                                                                                              >  2:339047/1‑150 (MQ=255)
                                         cAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGa                                                                                                              >  2:608665/1‑150 (MQ=255)
                                           aaaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGAtt                                                                                                            >  1:558655/1‑150 (MQ=255)
                                             aaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCt                                                                                                          <  2:521712/150‑1 (MQ=255)
                                                tACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAACGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaa                                                                                                       <  1:628190/150‑1 (MQ=255)
                                                 aCTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGGCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATTATTCTaaaa                                                                                                      >  1:1122944/1‑150 (MQ=255)
                                                  cTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAAc                                                                                                     <  2:12826/150‑1 (MQ=255)
                                                    gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                                   <  2:1276845/150‑1 (MQ=255)
                                                       cGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaata                                                                                                >  2:493465/1‑150 (MQ=255)
                                                         aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGTCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                              <  1:24666/150‑1 (MQ=255)
                                                         aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                              >  1:966154/1‑150 (MQ=255)
                                                         aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                              >  2:823993/1‑150 (MQ=255)
                                                           ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                            <  2:1353275/150‑1 (MQ=255)
                                                               aTTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGATGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtct                                                                                        >  1:1593247/1‑150 (MQ=255)
                                                                tttGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtctc                                                                                       <  1:751570/150‑1 (MQ=255)
                                                                  tGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTAtctctc                                                                                     >  1:1000692/1‑150 (MQ=255)
                                                                   gCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATctctct                                                                                    >  1:601507/1‑150 (MQ=255)
                                                                     ttGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCttt                                                                                  >  1:311252/1‑150 (MQ=255)
                                                                     ttGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCttt                                                                                  <  2:1190861/150‑1 (MQ=255)
                                                                        aaGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaa                                                                               >  1:1708905/1‑150 (MQ=255)
                                                                         aGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaaa                                                                              >  2:892344/1‑150 (MQ=255)
                                                                         aGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTATTTTaaa                                                                              >  2:886376/1‑150 (MQ=255)
                                                                          gCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAg                                                                             <  2:322283/150‑1 (MQ=255)
                                                                            gTTGATAATTGAGGAAATTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGATGTAATCTGAAAAAGGTGCTATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                                                           <  2:936872/150‑1 (MQ=255)
                                                                              tGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                                                         >  2:1110370/1‑150 (MQ=255)
                                                                                  aaTTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCa                                                                     >  2:1562796/1‑150 (MQ=255)
                                                                                     tGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTg                                                                  >  1:10129/1‑150 (MQ=255)
                                                                                     tGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTg                                                                  <  2:1682568/150‑1 (MQ=255)
                                                                                      gAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTgg                                                                 >  2:258435/1‑150 (MQ=255)
                                                                                       aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAGAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                                >  2:1533250/1‑150 (MQ=255)
                                                                                       aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                                <  2:766453/150‑1 (MQ=255)
                                                                                       aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                                <  2:689433/150‑1 (MQ=255)
                                                                                        ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAGATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                               <  2:628217/150‑1 (MQ=255)
                                                                                        ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                               >  2:112124/1‑150 (MQ=255)
                                                                                           aaaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAATCAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGtt                                                            >  1:832906/1‑150 (MQ=255)
                                                                                             aTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttct                                                          <  2:1599405/150‑1 (MQ=255)
                                                                                                aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCATTTCATTGGCTGTTCTTCg                                                       >  2:459520/1‑150 (MQ=255)
                                                                                                   aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                                    >  1:1702931/1‑150 (MQ=255)
                                                                                                    aGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACtt                                                   >  1:251214/1‑150 (MQ=255)
                                                                                                     gCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTACTTCGACTTa                                                  <  2:760404/150‑1 (MQ=255)
                                                                                                      ccGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTaa                                                 >  2:1532696/1‑150 (MQ=255)
                                                                                                        gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAcc                                               <  1:1428873/150‑1 (MQ=255)
                                                                                                            gCGACGTAACCAGGCCCGAGGTAATCTGAAAAATGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGTGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCaa                                           <  1:1369290/150‑1 (MQ=255)
                                                                                                             cGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAAt                                          >  2:105051/1‑150 (MQ=255)
                                                                                                               aagTAACCAGGCCCGAGGTAATCTGAAAAATGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCt                                        <  2:1287139/148‑1 (MQ=255)
                                                                                                               aCGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCt                                        >  1:253675/1‑150 (MQ=255)
                                                                                                                cGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTg                                       <  2:754552/150‑1 (MQ=255)
                                                                                                                 ttAACCAGGCCCGAGGTAATATGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGc                                      <  1:1371410/149‑1 (MQ=255)
                                                                                                                 gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGc                                      <  1:1262978/150‑1 (MQ=255)
                                                                                                                     ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATGTCTCTTTTAAAGTCAATTCATTGGATGTTCTGCGACTTAACCTCAATCTGCATCt                                  >  2:282872/1‑150 (MQ=255)
                                                                                                                     ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCt                                  >  1:520753/1‑150 (MQ=255)
                                                                                                                      cAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCtt                                 >  1:445484/1‑150 (MQ=255)
                                                                                                                        ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                               >  1:1536480/1‑150 (MQ=255)
                                                                                                                          cccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCat                             <  2:355787/150‑1 (MQ=255)
                                                                                                                           ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGTCTGTTCTTCGACTTAACCTCAATCTGCATATTGCata                            >  1:1204136/1‑150 (MQ=255)
                                                                                                                           ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                            <  1:947171/150‑1 (MQ=255)
                                                                                                                            cGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCataa                           <  1:451481/150‑1 (MQ=255)
                                                                                                                                gTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATc                       <  1:1158496/150‑1 (MQ=255)
                                                                                                                                 tAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCt                      <  2:1567644/150‑1 (MQ=255)
                                                                                                                                 tAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAACTCAATCTGCATCTTGCATAATATCt                      >  1:1170320/1‑150 (MQ=255)
                                                                                                                                   acctgttaaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                    >  2:559562/8‑150 (MQ=39)
                                                                                                                                   aTCTGAAAAAGGTGCAATCTTTTGCAACTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc                    <  1:788097/150‑1 (MQ=255)
                                                                                                                                     cTGAATAAGGTGCAATCTTTTGCAGCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACATATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAg                  <  2:449265/150‑1 (MQ=255)
                                                                                                                                      tGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt                 >  2:768398/1‑150 (MQ=255)
                                                                                                                                         taaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt              >  2:1016951/2‑150 (MQ=255)
                                                                                                                                              tgTGCAATCTTTTGCACCTTTAACAGTTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTc          <  1:24616/149‑1 (MQ=255)
                                                                                                                                              ggTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTc          <  1:1243935/150‑1 (MQ=255)
                                                                                                                                               gTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCa         <  1:382595/150‑1 (MQ=255)
                                                                                                                                                 gCAATATTTTTCACCTTTAACAGGTTAGCTCGGGCTAAAATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATa       <  2:1695358/150‑1 (MQ=255)
                                                                                                                                                   aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGc     >  1:1081253/1‑150 (MQ=255)
                                                                                                                                                   aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTATCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGc     >  1:514432/1‑150 (MQ=255)
                                                                                                                                                    aTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGCg    <  1:1068154/150‑1 (MQ=255)
                                                                                                                                                      cTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAACACAAAAAAATTATATCTCTTTTAAAGTCAATTCATTTGCTGTTCTTCGCCTTAACCTCAATCTGCATCTTGAATAATCTCTGCAGTAATCTTCATAGCGAc  <  2:328061/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGCGAC  >  NC_000913/4281529‑4281827

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: