Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,425,968 T→G 6.0% intergenic (‑52/+38) hisQ ← / ← hisJ histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,425,9680TG 6.0% 241.7 / 6.5 84intergenic (‑52/+38)hisQ/hisJhistidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/39);  new base G (3/2);  total (43/41)
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.80e-01

ATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTA  >  NC_000913/2425819‑2426114
                                                                                                                                                     |                                                                                                                                                  
aTTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGt                                                                                                                                                    <  1:79295/150‑1 (MQ=255)
 ttAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGtt                                                                                                                                                   <  2:1486984/150‑1 (MQ=255)
 ttAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGtt                                                                                                                                                   <  1:1141534/150‑1 (MQ=255)
     cccGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTAccc                                                                                                                                               <  1:1647967/149‑1 (MQ=255)
     aCCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTAtcc                                                                                                                                               <  2:873343/150‑1 (MQ=255)
      ccGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGGCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCg                                                                                                                                              >  1:135747/1‑150 (MQ=255)
        gatgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                            <  2:937842/150‑1 (MQ=255)
          tgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAgg                                                                                                                                          >  1:1218815/1‑150 (MQ=255)
           gcttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGAGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAggg                                                                                                                                         <  1:595385/148‑1 (MQ=255)
            atTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGa                                                                                                                                        <  2:1487015/150‑1 (MQ=255)
             ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAg                                                                                                                                       <  2:774998/150‑1 (MQ=255)
              tACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTTCGGGTAACCggagggaga                                                                                                                                      <  1:1081023/150‑2 (MQ=255)
               aCAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGATAACCCATACAACACGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCg                                                                                                                                     <  1:107373/150‑1 (MQ=255)
                cAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTTCGGGTAACCggagggagaca                                                                                                                                    <  2:502053/150‑4 (MQ=255)
                   cGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGct                                                                                                                                 <  2:1258040/150‑1 (MQ=255)
                   cGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGct                                                                                                                                 >  2:969062/1‑150 (MQ=255)
                      gCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACAAGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCATGTCTCCCTCCGGTTACCCGAAGGGTGCGGGctttc                                                                                                                              >  1:1016552/1‑150 (MQ=255)
                        aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttctt                                                                                                                            <  1:380492/150‑1 (MQ=255)
                         cTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                           >  2:1641720/1‑150 (MQ=255)
                         cTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                           >  2:1519653/1‑150 (MQ=255)
                         cTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttcttt                                                                                                                           <  2:629014/150‑1 (MQ=255)
                          tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                          <  1:3485/150‑1 (MQ=255)
                          tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                          <  2:964783/150‑1 (MQ=255)
                            cAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAc                                                                                                                        <  2:893427/150‑1 (MQ=255)
                               agCTGATAGCCAGCTCCAGCGTGAGGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACACCATGCTGCCTGTCCTGTCGTGTGTGGGGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGa                                                                                                                     <  1:71864/150‑1 (MQ=255)
                               agCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGa                                                                                                                     <  1:696313/150‑1 (MQ=255)
                                         cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCGAGTCGTGTGTGGTGAATTACCGCCTCTTGGCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAg                                                                                                           <  1:320047/150‑1 (MQ=255)
                                         cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAg                                                                                                           >  2:124740/1‑150 (MQ=255)
                                                 gCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGGTGGCTAATCGccatcggt                                                                                                   <  2:1372281/150‑6 (MQ=255)
                                                 gCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATa                                                                                                   <  2:135747/150‑1 (MQ=255)
                                                     gACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCGTTCACCGATGGCGATTAGCCACCATAAACa                                                                                               >  2:1413395/1‑150 (MQ=255)
                                                        gAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACAtaa                                                                                            >  1:144386/1‑148 (MQ=255)
                                                            gcACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAAt                                                                                        <  1:857881/150‑1 (MQ=255)
                                                                 cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGGCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                   <  1:103508/150‑1 (MQ=255)
                                                                  tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTATCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGt                                                                                  >  1:728158/1‑150 (MQ=255)
                                                                         tAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCACTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACttttt                                                                           >  1:1594620/1‑150 (MQ=255)
                                                                           acacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCg                                                                         <  1:1631977/150‑1 (MQ=255)
                                                                           acacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCg                                                                         >  2:1082497/1‑150 (MQ=255)
                                                                             acCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCt                                                                       >  2:847132/1‑150 (MQ=255)
                                                                              ccTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCACCGGTTACCCGAAGGGAGCGGGCTTTCTTTAACCGATGGAGATTAGCCACCATAAACATCAAAATCGAAGAACTTTTTCGATa                                                                      >  1:178321/1‑150 (MQ=255)
                                                                                 gAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAtt                                                                   <  1:1366348/150‑1 (MQ=255)
                                                                                   aaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtc                                                                 <  1:969062/150‑1 (MQ=255)
                                                                                     aCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                               <  2:1639935/150‑1 (MQ=255)
                                                                                      cccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                              >  1:1024709/1‑150 (MQ=255)
                                                                                      cccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCACCCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                              >  1:5146/1‑150 (MQ=255)
                                                                                       ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGt                                                             >  1:879219/1‑150 (MQ=255)
                                                                                        cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCCCTAATTTCTCGTa                                                            >  1:1078873/1‑150 (MQ=255)
                                                                                          tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAg                                                          >  1:516488/1‑150 (MQ=255)
                                                                                          tACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAACGATGGTGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTATAGTaat                                                          >  1:884067/1‑149 (MQ=255)
                                                                                          tACAACATGCTGACTGTCCTGTCGGGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAg                                                          <  2:1150607/150‑1 (MQ=255)
                                                                                            caacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTa                                                        >  1:604883/1‑150 (MQ=255)
                                                                                                  gctgcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAATGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCa                                                  <  1:173060/150‑1 (MQ=255)
                                                                                                     gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTAcctgacgcg                                               >  2:876667/1‑143 (MQ=255)
                                                                                                      cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCa                                                  >  1:1523403/1‑146 (MQ=255)
                                                                                                      cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCa                                                  <  2:1523403/146‑1 (MQ=255)
                                                                                                      cctgtcctgtCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgc                                              <  2:1654311/150‑1 (MQ=255)
                                                                                                         gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                             >  2:268466/1‑148 (MQ=255)
                                                                                                         gtcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                             <  1:268466/148‑1 (MQ=255)
                                                                                                           cctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTACCTCCGGTTACCCGAAGGGAGCGGGCGTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTATCGTAAGTACCGTCAGCGCGCAtt                                         >  1:1427635/1‑150 (MQ=255)
                                                                                                            ctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAttt                                        >  1:553243/1‑150 (MQ=255)
                                                                                                            ctgtcGTGTGTGGTGAATTACAGCCTATTGTCTCCATACGGTTACCCGAAGGGAGAGTGCTTTCTTTCACCGATGGCGATTTGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTAGCGTCAGCGCGCAttt                                        >  2:215354/1‑150 (MQ=255)
                                                                                                             tgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                       <  2:158462/150‑1 (MQ=255)
                                                                                                             tgtcGTGTGTGGTGAATTACCGCCACTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                       >  2:304891/1‑150 (MQ=255)
                                                                                                             agtcGTGTGTGGTGAATTACAGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGATCGGGCTTTCTTTCACCGATGGCTATTAGACACCATAAAAATCAAAATCGAAGTACTTATTCGCAAATTTCTCGTAAGTACCGTCAGCGCGCATATc                                       >  1:16701/2‑150 (MQ=255)
                                                                                                               tcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                     <  1:1005964/150‑1 (MQ=255)
                                                                                                               tcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                     <  1:847132/150‑1 (MQ=255)
                                                                                                                    tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                >  2:573233/1‑150 (MQ=255)
                                                                                                                       ggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                             >  1:679842/1‑150 (MQ=255)
                                                                                                                       ggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                             >  2:1002765/1‑150 (MQ=255)
                                                                                                                       ggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACACGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                             >  1:778190/1‑150 (MQ=255)
                                                                                                                        gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                            <  1:1519653/150‑1 (MQ=255)
                                                                                                                        gtgAATTACCGCCTATTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCtt                            >  1:167635/1‑150 (MQ=255)
                                                                                                                         tgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCCCCGATGGCGATTAGCCACCATAAACCTCACAATCGATCTACGTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGGAAAGGCttt                           >  2:1230589/1‑150 (MQ=255)
                                                                                                                             ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTc                       <  2:516488/150‑1 (MQ=255)
                                                                                                                                ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGt                    <  1:124740/150‑1 (MQ=255)
                                                                                                                                 cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTg                   >  1:883470/1‑150 (MQ=255)
                                                                                                                                  gCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTATCACCGATGGCGATTAGCCACCATAACCATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGAAAAGGCTTTGTTCAGTgg                  >  1:94862/1‑149 (MQ=255)
                                                                                                                                      ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcg              >  2:4240/1‑150 (MQ=255)
                                                                                                                                         gTCTCCCTCCGGTTACCCGAAGGGAGCGGGATTTCTTTCACCGATGGCGATTAGCCACCATAAATATCAAAATCTTAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCgcgc           <  2:62407/150‑1 (MQ=255)
                                                                                                                                         cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCgcgc           >  2:739449/3‑150 (MQ=255)
                                                                                                                                          gctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCa          >  2:725394/2‑150 (MQ=255)
                                                                                                                                           ctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAg         >  2:511076/1‑150 (MQ=255)
                                                                                                                                             ccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTTCTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGtt       <  1:1003004/150‑1 (MQ=255)
                                                                                                                                               ctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg     >  1:1588898/1‑150 (MQ=255)
                                                                                                                                               ctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg     >  1:376999/1‑150 (MQ=255)
                                                                                                                                                  cGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTa  <  2:883470/150‑1 (MQ=255)
                                                                                                                                                  cGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTa  <  2:1687708/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTA  >  NC_000913/2425819‑2426114

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