Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,677 T→A 5.2% 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 5.2% 480.3 / 7.6 159intergenic (+48/+106)yjcE/yjcFputative cation/proton antiporter/pentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base T (83/67);  new base A (4/5);  total (87/72)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGC  >  NC_000913/4281528‑4281824
                                                                                                                                                     |                                                                                                                                                   
tACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAt                                                                                                                                                     >  2:1211872/1‑150 (MQ=255)
  ccaccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCt                                                                                                                                                   <  2:1400128/150‑1 (MQ=255)
  ccaccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCt                                                                                                                                                   >  1:1419450/1‑150 (MQ=255)
    accTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                                 >  1:885463/1‑150 (MQ=255)
      cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                               >  2:6829/1‑150 (MQ=255)
         cgcgcCACACGGGAGATCAGCAACGAAATGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAc                                                                                                                                            >  2:530734/1‑150 (MQ=255)
           cgcCACCCGGGAGATCAGAAACGAAACGATGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGAGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCt                                                                                                                                          <  2:865530/150‑1 (MQ=255)
           cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCt                                                                                                                                          <  2:1309376/150‑1 (MQ=38)
             ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                        <  2:921620/150‑1 (MQ=255)
             ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                        >  2:590726/1‑150 (MQ=255)
             ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                        >  1:1380235/1‑150 (MQ=255)
             ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCttt                                                                                                                                        >  1:686631/1‑150 (MQ=255)
              cacaCGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTa                                                                                                                                       <  2:46417/150‑1 (MQ=255)
               acacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTaa                                                                                                                                      >  2:1023994/1‑150 (MQ=255)
                cacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtttttc                                                                                                                                     <  2:896049/150‑4 (MQ=255)
                            gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGTGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                         >  1:580125/1‑150 (MQ=255)
                            gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                         >  1:1466982/1‑150 (MQ=255)
                            gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                         <  2:1257181/150‑1 (MQ=255)
                            gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCgg                                                                                                                         <  2:410730/150‑1 (MQ=25)
                              aaCGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGc                                                                                                                       <  1:1152736/150‑1 (MQ=255)
                              aaCGAAACGCTGCAAAAATTACTGCACGATCGCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCATAAAGCGACGAAACCTGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGc                                                                                                                       <  1:1556061/150‑1 (MQ=255)
                                   aaCGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTa                                                                                                                               <  2:1139583/137‑1 (MQ=255)
                                   aaCGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTa                                                                                                                               >  1:1139583/1‑137 (MQ=255)
                                   aaCGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACa                                                                                                                  >  1:844705/1‑150 (MQ=255)
                                       cTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCa                                                                                                              >  1:1093749/1‑150 (MQ=255)
                                           aaaaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGAt                                                                                                          >  1:918336/1‑150 (MQ=255)
                                            aaaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGAtt                                                                                                         >  1:1491115/1‑150 (MQ=255)
                                            aaaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGAtt                                                                                                         >  2:493679/1‑150 (MQ=255)
                                                 tACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaa                                                                                                    <  2:1107865/150‑1 (MQ=255)
                                                    tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACa                                                                                                 <  2:1586168/150‑1 (MQ=255)
                                                    tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGAGACGTAACCAGGACCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACa                                                                                                 <  2:1579753/150‑1 (MQ=255)
                                                        cGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaaga                                                                                             >  2:294760/1‑148 (MQ=255)
                                                           tctcGATTTGCTTGAAGCGTTGCTAATTGAGGAATATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                         <  2:1700895/150‑1 (MQ=255)
                                                            ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAttt                                                                                         >  1:1506288/1‑150 (MQ=255)
                                                                    gCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATctctct                                                                                 >  1:1458939/1‑150 (MQ=255)
                                                                     cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaata                                                                                             <  1:990065/137‑1 (MQ=255)
                                                                     cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaata                                                                                             >  2:990065/1‑137 (MQ=255)
                                                                     cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtt                                                                                <  2:1515350/150‑1 (MQ=255)
                                                                        gAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTa                                                                             >  1:1198789/1‑150 (MQ=255)
                                                                         aaGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaa                                                                            <  2:1093749/150‑1 (MQ=255)
                                                                             gTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGta                                                                        >  1:624574/1‑149 (MQ=255)
                                                                               tGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                                                      >  1:73627/1‑150 (MQ=255)
                                                                               agCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                                                      >  1:73631/2‑150 (MQ=255)
                                                                                  tAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTc                                                                   >  1:1361673/1‑150 (MQ=255)
                                                                                   aaTTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCa                                                                  >  2:597300/1‑150 (MQ=255)
                                                                                      tGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTg                                                               >  1:384426/1‑150 (MQ=255)
                                                                                       gAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAACGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTgg                                                              <  1:1210237/150‑1 (MQ=255)
                                                                                        aGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGc                                                             <  1:403144/150‑1 (MQ=255)
                                                                                         ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            <  1:1023994/150‑1 (MQ=255)
                                                                                         ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            >  2:524019/1‑150 (MQ=255)
                                                                                         ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            >  2:288621/1‑150 (MQ=255)
                                                                                         ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                            >  2:1569116/1‑150 (MQ=255)
                                                                                          gAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTg                                                           <  1:757508/150‑1 (MQ=255)
                                                                                          gAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTg                                                           >  1:1144945/1‑150 (MQ=255)
                                                                                           aaaaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAGATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGt                                                          <  2:816143/150‑1 (MQ=255)
                                                                                            aaaTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGtt                                                         <  1:1618624/150‑1 (MQ=255)
                                                                                                atttAAGCCGATAGCGACGTAACCAGTCACGAGGTAATCTGAATAAGTTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCTTGATTCTAAAACATATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttt                                                     <  1:1238165/147‑2 (MQ=255)
                                                                                                 aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                                    <  1:524019/150‑1 (MQ=255)
                                                                                                 aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                                    <  2:518003/150‑1 (MQ=255)
                                                                                                 aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                                    >  2:1476183/1‑150 (MQ=255)
                                                                                                 aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATATCTCTTTTAAATTCAATTCATTGGCTGttcttct                                                    >  2:1572315/1‑149 (MQ=255)
                                                                                                    aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                                 >  1:349711/1‑150 (MQ=255)
                                                                                                    aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGATCTTCGAca                                                 >  1:344374/1‑149 (MQ=255)
                                                                                                     aGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACtt                                                >  1:1071876/1‑150 (MQ=255)
                                                                                                       ccGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTaa                                              >  2:878039/1‑150 (MQ=255)
                                                                                                        cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                             <  2:1466982/150‑1 (MQ=255)
                                                                                                        cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                             <  1:1184250/150‑1 (MQ=255)
                                                                                                          aaaGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAcat                                           >  2:1522404/1‑148 (MQ=255)
                                                                                                              cGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAAt                                       >  2:822854/1‑150 (MQ=255)
                                                                                                              cGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAAt                                       >  2:311436/1‑150 (MQ=255)
                                                                                                               gACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATc                                      >  1:1672231/1‑150 (MQ=255)
                                                                                                                aCGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCt                                     <  1:1095401/150‑1 (MQ=255)
                                                                                                                   tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                                  >  1:556469/1‑150 (MQ=255)
                                                                                                                   tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                                  <  1:878039/150‑1 (MQ=255)
                                                                                                                     aCCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCAGCAta                                >  1:612367/1‑149 (MQ=255)
                                                                                                                       cAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCtt                              <  2:1198789/150‑1 (MQ=255)
                                                                                                                        aGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTg                             <  1:1596303/150‑1 (MQ=255)
                                                                                                                        aGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTA‑‑TCTCATTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCa                           <  2:957945/150‑1 (MQ=255)
                                                                                                                         ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                            <  1:1581926/150‑1 (MQ=255)
                                                                                                                         ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                            >  1:26585/1‑150 (MQ=255)
                                                                                                                           cccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCat                          >  1:1335081/1‑150 (MQ=255)
                                                                                                                            ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                         <  2:1353965/150‑1 (MQ=255)
                                                                                                                               aGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCataata                      >  1:924089/1‑150 (MQ=255)
                                                                                                                                     tCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCa                <  2:1143175/150‑1 (MQ=255)
                                                                                                                                      cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCat               >  1:1534411/1‑149 (MQ=255)
                                                                                                                                      cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAg               <  2:73631/150‑1 (MQ=255)
                                                                                                                                      cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAg               <  1:960475/150‑1 (MQ=255)
                                                                                                                                      cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAg               <  2:73627/150‑1 (MQ=255)
                                                                                                                                       tGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt              >  1:1707883/1‑150 (MQ=255)
                                                                                                                                        gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTa             >  2:681381/1‑150 (MQ=255)
                                                                                                                                        gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTa             >  2:222447/1‑150 (MQ=255)
                                                                                                                                          taaaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTTACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt           >  2:1621766/2‑150 (MQ=255)
                                                                                                                                            aaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAAAATCTGCAGTAATCt         >  1:186645/1‑150 (MQ=255)
                                                                                                                                            aaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCt         >  2:105586/1‑150 (MQ=255)
                                                                                                                                             aGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTTATCtt        >  1:252580/1‑150 (MQ=255)
                                                                                                                                                tGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCAt     >  2:703844/1‑150 (MQ=255)
                                                                                                                                                 gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATa    <  1:311436/150‑1 (MQ=255)
                                                                                                                                                 acAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTGAATCTGCATGTTGCATAATATCTGCAGTAATCTTCATa    >  2:249035/2‑150 (MQ=255)
                                                                                                                                                  cAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAg   <  2:580125/150‑1 (MQ=255)
                                                                                                                                                   aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGc  <  1:795546/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATAGC  >  NC_000913/4281528‑4281824

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: