Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,049 T→A 5.2% intergenic (‑436/+23) yhbX ← / ← leuU putative EptAB family phosphoethanolamine transferase, inner membrane protein/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,0490TA 5.2% 299.3 / 3.1 101intergenic (‑436/+23)yhbX/leuUputative EptAB family phosphoethanolamine transferase, inner membrane protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base T (46/49);  new base A (4/2);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCG  >  NC_000913/3321902‑3322195
                                                                                                                                                   |                                                                                                                                                  
aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                                  <  2:133710/150‑1 (MQ=255)
  ttATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAg                                                                                                                                                <  1:496946/150‑1 (MQ=255)
    aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                              <  1:750681/150‑1 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATCCCTGACAACCATTCCACAGCCCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                                                                                                                            >  2:2008552/1‑150 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                                                                                                                            >  2:1982361/1‑150 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGta                                                                                                                                            >  2:703525/1‑149 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCAtc                                                                                                                                            <  2:1823752/150‑1 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTATCAGCGtc                                                                                                                                            <  1:1426827/150‑1 (MQ=255)
        cAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtctc                                                                                                                                          >  1:179069/1‑150 (MQ=255)
        cAGTAGAAGGGAAAAATTTATAACTAACGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtctc                                                                                                                                          <  2:1458532/150‑1 (MQ=255)
           tAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCttt                                                                                                                                       <  1:1550644/150‑1 (MQ=255)
            aGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTGCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTGAAAAGCGTCtttttt                                                                                                                                      <  1:1849111/150‑1 (MQ=255)
             gAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTc                                                                                                                                     >  2:1463966/1‑150 (MQ=255)
               aGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGCCCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTg                                                                                                                                   >  2:1369314/1‑150 (MQ=255)
                  gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                                >  2:552194/1‑150 (MQ=255)
                     aaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACCCTTTTCAGCGTCTCTTTACTGGAAttt                                                                                                                             >  1:337017/1‑150 (MQ=255)
                       aTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTgg                                                                                                                           <  1:1634707/150‑1 (MQ=255)
                        tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTTCCAGCTACGAGTAAAGCAACTGGACGAGATACTGATACCTGACAACCATTCCTCAGACCAGGACCCAAACGAAAAACGCCGCTTATCCGCGGCTCGTTTCTGGATGTTGGt                                                                                                                          >  1:1599775/1‑150 (MQ=255)
                           ataACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTAcc                                                                                                                       <  2:294513/150‑1 (MQ=255)
                            taACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCg                                                                                                                      <  2:1881344/150‑1 (MQ=255)
                               cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGATTTTCAGCGTCTCTTTTCTGGAATATGGTACCGAgg                                                                                                                   >  1:1086173/1‑150 (MQ=255)
                                      gTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGAc                                                                                                            <  2:583334/150‑1 (MQ=255)
                                       tACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACt                                                                                                           <  2:369975/150‑1 (MQ=255)
                                         ccGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCGTGTCTGGAATTTGGTACCGAGGACGGGACttt                                                                                                         >  1:509110/1‑149 (MQ=255)
                                         ccGGCACCATCGATTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTATCTTTTCTGGAAATTGGTACCGAGGACGGGACTTg                                                                                                         <  2:945683/150‑1 (MQ=255)
                                           ggCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGaa                                                                                                       >  1:1709208/1‑150 (MQ=255)
                                            gCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCCATCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAc                                                                                                      <  1:288824/150‑1 (MQ=255)
                                             cACCATCGTTTCAAGGTACCAGCTACGAGGAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGAACAGGACCAAAACGATAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGTAcc                                                                                                     <  1:86232/150‑1 (MQ=255)
                                             cACCATAGTTTCAAGTTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGTACTTGAAcc                                                                                                     <  2:1673047/150‑1 (MQ=255)
                                               ccATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCg                                                                                                   <  1:1239709/150‑1 (MQ=255)
                                                 aTCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTa                                                                                                 <  1:252370/150‑1 (MQ=255)
                                                   cGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGACAAACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCGCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAg                                                                                               <  1:1008281/150‑1 (MQ=255)
                                                     tttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTAGAACCCGTAAGcc                                                                                             >  2:617975/1‑150 (MQ=255)
                                                       tCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGccca                                                                                           >  1:1689566/1‑149 (MQ=255)
                                                        cAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTa                                                                                          >  2:932306/1‑150 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                    <  2:913318/150‑1 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCATTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGTCCGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                    >  1:84710/1‑150 (MQ=255)
                                                               ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGCCTTGAACCCGTAAGCCCTATTGGGcc                                                                                   >  2:746838/1‑149 (MQ=255)
                                                               ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCa                                                                                   >  1:1100655/1‑150 (MQ=255)
                                                                cAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCAc                                                                                  <  1:1825913/150‑1 (MQ=255)
                                                                cAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCATCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCAc                                                                                  <  1:486630/150‑1 (MQ=255)
                                                                     aCGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacca                                                                             <  2:599021/150‑1 (MQ=255)
                                                                     aCGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacca                                                                             <  1:798249/150‑1 (MQ=255)
                                                                     aCGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacca                                                                             <  2:1944332/150‑1 (MQ=255)
                                                                        aGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCt                                                                          <  2:1534655/150‑1 (MQ=255)
                                                                         gTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTc                                                                         >  2:1881130/1‑150 (MQ=255)
                                                                          tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                        <  1:52348/150‑1 (MQ=255)
                                                                            aaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAg                                                                      >  1:1614185/1‑150 (MQ=255)
                                                                             aGCAACTGGACGAGATACAGATACCTGACAACCATTCCTAAGACCACGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTATGGAATTAGGTAACGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACAACCACCTCAAgg                                                                     >  1:2002529/1‑150 (MQ=255)
                                                                              gCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGt                                                                    >  2:432956/1‑150 (MQ=255)
                                                                                     gACGAGATACAGATACCTGACAACCCTTCCTCAGACCAGTACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCTAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtgt                                                             <  1:356044/150‑1 (MQ=255)
                                                                                      aCGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTCAGCCCTATTGGGCACTACCACCTCAAGGTAGCGGGTc                                                            >  2:1252575/1‑150 (MQ=255)
                                                                                        gagaTACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTa                                                          <  1:1463966/150‑1 (MQ=255)
                                                                                          tatacagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTAcc                                                        <  1:1855855/149‑1 (MQ=255)
                                                                                          gatacagatacCTGACAACCATTCCTCAGACCAGCACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTAcc                                                        <  1:259120/150‑1 (MQ=255)
                                                                                                tatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTcc                                                  <  1:1525830/149‑1 (MQ=255)
                                                                                                    ccTGACAACCATTCCTCAGACCCGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGTCTTGAACCCGTAAGCCCTATTGGGCCCTACCCCCTCAAGGTAGCGTGTCTACCAATTccacaa                                              >  2:192264/1‑148 (MQ=255)
                                                                                                     cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                             >  2:674686/1‑150 (MQ=255)
                                                                                                      tGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGGGTCTACCAATTACaccacc                                            >  1:903783/1‑150 (MQ=255)
                                                                                                       gACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGACCCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCt                                           >  2:263237/1‑150 (MQ=255)
                                                                                                          aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                        <  1:1359621/150‑1 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                       >  2:1904454/1‑150 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAGTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTAAGCCACCTCGGc                                       >  2:813667/1‑150 (MQ=255)
                                                                                                             cATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCCCCTCAAGGTAGCGTGTCTGCCAATTCCACCACCTCGGCAc                                     >  2:184653/1‑150 (MQ=255)
                                                                                                             cATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAATGGGCACTACCACCTCAAGGTAGCGTGTGTACCAATTCCACCACCTCGGCAc                                     >  1:1863820/1‑150 (MQ=255)
                                                                                                              aTTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACg                                    >  2:645267/1‑150 (MQ=255)
                                                                                                                tCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGa                                  <  1:456990/150‑1 (MQ=255)
                                                                                                                tCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGa                                  <  2:1100655/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTGTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                 <  1:800423/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                 >  2:547096/1‑150 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                 <  1:1432564/150‑1 (MQ=255)
                                                                                                                    cAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtact                              <  1:432956/150‑1 (MQ=255)
                                                                                                                    cAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtact                              >  2:673718/1‑150 (MQ=255)
                                                                                                                      gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                            >  2:192766/1‑150 (MQ=255)
                                                                                                                          aGGACCAAAACGAAAAAAGCCGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCAGCACGGATACTACtttt                        <  1:1369314/150‑1 (MQ=255)
                                                                                                                           ggACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACttttt                       >  2:342934/1‑150 (MQ=255)
                                                                                                                           ggACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACttttt                       <  1:1881130/150‑1 (MQ=255)
                                                                                                                            gACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTa                      <  1:1732733/150‑1 (MQ=255)
                                                                                                                            gACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTa                      <  2:1104021/150‑1 (MQ=255)
                                                                                                                             aCCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTTGAATTTTGTACCGATGACGGGACTTGTACCCGTAATGCCTATTGGGCACTACCTCCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAg                     <  2:2002529/150‑1 (MQ=255)
                                                                                                                             aCCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAg                     <  1:571695/150‑1 (MQ=255)
                                                                                                                             aCCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAg                     <  2:1681346/150‑1 (MQ=255)
                                                                                                                             aCCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAAATCCACCCCCTCGGCACGGATACTACTTTTTAg                     >  2:1793172/1‑150 (MQ=255)
                                                                                                                              ccAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGt                    >  1:1404988/1‑150 (MQ=255)
                                                                                                                               cAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGtt                   >  2:615778/1‑150 (MQ=255)
                                                                                                                               cAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGtt                   >  2:539227/1‑150 (MQ=255)
                                                                                                                                aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCCCCTCGGCACGGATACTACTTTTTAGTTc                  >  2:1857660/1‑150 (MQ=255)
                                                                                                                                 acactaaaaacGACGCTTTTCCTCGTCTCTTTTCTGGAATTTGGTACCGAGGAAGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCg                 <  2:322093/145‑1 (MQ=255)
                                                                                                                                   cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg               >  1:1724615/3‑150 (MQ=255)
                                                                                                                                   cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGGTCggg               >  1:1467002/3‑150 (MQ=255)
                                                                                                                                   cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCACTATTGGGCACTACCACCTCAAGGTAGCGTGTATACCAATTCCACCACCTCGGAACGGATACTACATTTTAGTTCgcg               >  1:1467003/3‑148 (MQ=255)
                                                                                                                                   aCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg               <  2:179069/150‑1 (MQ=255)
                                                                                                                                   aCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg               >  2:1565004/1‑150 (MQ=255)
                                                                                                                                    cGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGa              >  2:1588628/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGGACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGata            >  2:543556/1‑150 (MQ=255)
                                                                                                                                            gACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTTGGGATATCGCTg      >  1:1431473/1‑150 (MQ=255)
                                                                                                                                            gACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTg      <  1:1560498/150‑1 (MQ=255)
                                                                                                                                               gCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTc   <  1:1554582/150‑1 (MQ=255)
                                                                                                                                                cTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCg  <  2:1448153/150‑1 (MQ=255)
                                                                                                                                                cTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCAAGGATACTACTTTTTAGTTCGGGATATCGCTGGTCg  >  2:544444/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCG  >  NC_000913/3321902‑3322195

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: