Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,854,024 C→A 7.5% intergenic (‑45/+72) ydjE ← / ← ydjF putative transporter YdjE/putative DNA‑binding transcriptional regulator YdjF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,854,0240CA 7.5% 349.1 / 11.1 105intergenic (‑45/+72)ydjE/ydjFputative transporter YdjE/putative DNA‑binding transcriptional regulator YdjF
Reads supporting (aligned to +/- strand):  ref base C (54/43);  new base A (4/4);  total (58/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 107 total aligned reads displayed.

CCCTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGTA  >  NC_000913/1853875‑1854173
                                                                                                                                                     |                                                                                                                                                     
  cTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCACCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                     >  2:150930/1‑150 (MQ=255)
    gTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACtt                                                                                                                                                   <  1:511098/150‑1 (MQ=255)
      taataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGc                                                                                                                                                 <  1:1048320/150‑1 (MQ=255)
        ataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCg                                                                                                                                               >  2:1963160/1‑150 (MQ=255)
        ataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCg                                                                                                                                               >  1:175459/1‑150 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCAt                                                                                                                                            >  2:1676950/1‑150 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCAt                                                                                                                                            >  1:1375783/1‑150 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCAt                                                                                                                                            >  1:32117/1‑150 (MQ=255)
            cAGACTAAAGCTTATAATACCAAATATACGATAATGAAAGCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                           <  1:661093/150‑3 (MQ=255)
            cAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                           <  2:1526585/150‑3 (MQ=255)
            cAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATc                                                                                                                                           >  1:81762/1‑150 (MQ=255)
            cAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATc                                                                                                                                           <  1:586432/150‑1 (MQ=255)
              gACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCtt                                                                                                                                         >  2:1251081/1‑150 (MQ=255)
                  aaaGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTatc                                                                                                                                     <  2:1030908/150‑1 (MQ=255)
                     gCTTAAAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTATCAATTGCGGCAAATTCACTTGCCGCATCTTTatcatc                                                                                                                                  >  2:1238100/1‑150 (MQ=255)
                      cTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTAtcatca                                                                                                                                 <  1:1558766/150‑1 (MQ=255)
                              aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                         >  2:1159389/1‑150 (MQ=255)
                               ccAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTc                                                                                                                        >  1:93592/1‑150 (MQ=255)
                               ccAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTc                                                                                                                        <  2:1778455/150‑1 (MQ=255)
                                cAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCa                                                                                                                       >  2:1809722/1‑150 (MQ=255)
                                cAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCa                                                                                                                       <  2:1928859/150‑1 (MQ=255)
                                   atatACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGTGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGTTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAgg                                                                                                                    <  2:269089/150‑1 (MQ=255)
                                       aCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCa                                                                                                                >  1:716230/1‑150 (MQ=255)
                                       aCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCa                                                                                                                >  2:1896901/1‑150 (MQ=255)
                                        cGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAg                                                                                                               <  1:781648/150‑1 (MQ=255)
                                            aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgcg                                                                                                                  <  1:751744/143‑1 (MQ=255)
                                            aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgcg                                                                                                                  >  2:751744/1‑143 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          >  2:1463779/1‑150 (MQ=255)
                                                aaaCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaat                                                                                                       <  1:1968916/150‑1 (MQ=255)
                                                   ccGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaataat                                                                                                    >  2:562247/1‑150 (MQ=255)
                                                     gggCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAatat                                                                                                  <  1:1676950/150‑1 (MQ=255)
                                                      ggCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATAtt                                                                                                 <  2:282331/150‑1 (MQ=255)
                                                         cAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGAc                                                                                              >  2:214767/1‑150 (MQ=255)
                                                         cAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGAc                                                                                              >  1:999408/1‑150 (MQ=255)
                                                             ggCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACg                                                                                          <  2:81762/150‑1 (MQ=255)
                                                              gCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgt                                                                                         >  1:767270/1‑150 (MQ=255)
                                                               cAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCGAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgtg                                                                                        >  1:645738/1‑150 (MQ=255)
                                                               cAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgtg                                                                                        <  1:376741/150‑1 (MQ=255)
                                                                aaGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgtgt                                                                                       >  2:1323973/1‑150 (MQ=255)
                                                                 aGCGGTCCATTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                      <  2:1717350/150‑1 (MQ=255)
                                                                 aGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                      >  1:318583/1‑150 (MQ=255)
                                                                  gCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGc                                                                                     >  1:421266/1‑150 (MQ=255)
                                                                     gTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTa                                                                                  >  2:1261425/1‑150 (MQ=255)
                                                                         aggcTTGGGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGTATTTGCTAAAAAGGTAGCAATTGCTGCAAATTCACTTGCGGCATATTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTTACAACGTGTGCTTAAAtt                                                                              <  1:1238100/147‑1 (MQ=255)
                                                                         aGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTAGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCTTCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAAtt                                                                              >  2:1287141/1‑150 (MQ=255)
                                                                          gTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTa                                                                             <  2:1216648/150‑1 (MQ=255)
                                                                           tCTTGCACCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTAc                                                                            <  1:150930/150‑1 (MQ=255)
                                                                            cTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTAcc                                                                           <  2:93592/150‑1 (MQ=255)
                                                                                 gcCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACAcc                                                                      <  1:599461/150‑1 (MQ=255)
                                                                                  ccAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGACGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCa                                                                     >  1:1230190/1‑150 (MQ=255)
                                                                                     aTTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCt                                                                  <  1:1416578/150‑1 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     >  1:512486/1‑150 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     >  2:870450/1‑150 (MQ=255)
                                                                                                    tCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCtt                                                   >  2:1056795/1‑150 (MQ=255)
                                                                                                    tCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCtt                                                   >  2:148520/1‑150 (MQ=255)
                                                                                                     ccATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaata                                                                                                      <  1:1067752/98‑1 (MQ=255)
                                                                                                           ggTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCaaaa                                            >  2:1382497/1‑150 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGc                                         >  1:1128756/1‑150 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       >  1:1213441/1‑150 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       <  1:379438/150‑1 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       >  1:1638590/1‑150 (MQ=255)
                                                                                                                 tCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTAt                                      >  1:1645199/1‑150 (MQ=255)
                                                                                                                      tttGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCAct                                 >  1:1676996/1‑150 (MQ=255)
                                                                                                                        tGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGATATCCActct                               >  2:387765/1‑150 (MQ=255)
                                                                                                                         gCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTg                              <  2:175459/150‑1 (MQ=255)
                                                                                                                         gCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTg                              >  1:1790383/1‑150 (MQ=255)
                                                                                                                            aaaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGACCAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAc                           >  1:1281044/1‑150 (MQ=255)
                                                                                                                            aaaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAc                           <  2:716230/150‑1 (MQ=255)
                                                                                                                            aaaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAc                           >  2:356109/1‑150 (MQ=255)
                                                                                                                            aaaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAc                           >  1:1059735/1‑150 (MQ=255)
                                                                                                                              aaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAccc                         >  1:1044434/1‑150 (MQ=255)
                                                                                                                               aaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCg                        <  1:96640/150‑1 (MQ=255)
                                                                                                                                 ggTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGt                      <  1:573922/150‑1 (MQ=255)
                                                                                                                                   tAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGa                    >  2:259852/1‑150 (MQ=255)
                                                                                                                                    taaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCTAGGCGCAGGAACAAATAATATTGAAAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   >  2:94062/7‑150 (MQ=255)
                                                                                                                                    taaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCTAGGCGCAGGAACAAATAATATTGAAAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   <  1:94062/144‑1 (MQ=255)
                                                                                                                                    aGCAATTGCGGCGAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   <  2:645738/150‑1 (MQ=255)
                                                                                                                                     aaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAtt                  >  2:595220/6‑150 (MQ=255)
                                                                                                                                        gatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   >  1:440697/3‑146 (MQ=255)
                                                                                                                                        gatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   <  2:440697/144‑1 (MQ=255)
                                                                                                                                        gatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTAt               >  2:1287336/3‑150 (MQ=255)
                                                                                                                                         ttGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATc              <  1:771853/150‑1 (MQ=255)
                                                                                                                                         ttGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATc              <  1:723174/150‑1 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAg            >  1:844943/1‑150 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAg            >  1:1591537/1‑150 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAg            >  2:720358/1‑150 (MQ=255)
                                                                                                                                            cGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGt           <  1:1261425/150‑1 (MQ=255)
                                                                                                                                            cGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGt           <  1:1463779/150‑1 (MQ=255)
                                                                                                                                            cGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGt           <  2:1059735/150‑1 (MQ=255)
                                                                                                                                             ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt          >  1:1568296/1‑150 (MQ=255)
                                                                                                                                             ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACAAGCTGGATATTATTGTCTTTGCAAAATTCTATCCAATCTGCACCCGGTGATTTATAAGtt          >  2:847446/1‑150 (MQ=255)
                                                                                                                                              gCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtta         >  1:1241888/1‑150 (MQ=255)
                                                                                                                                                aaaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatt       <  1:1963160/150‑1 (MQ=255)
                                                                                                                                                aaaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatt       <  2:844943/150‑1 (MQ=255)
                                                                                                                                                aaaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatt       <  2:1591537/150‑1 (MQ=255)
                                                                                                                                                 aaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatta      >  1:1526189/1‑150 (MQ=255)
                                                                                                                                                 aaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatta      <  2:1790383/150‑1 (MQ=255)
                                                                                                                                                  aTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTtattat     <  1:214767/150‑1 (MQ=255)
                                                                                                                                                  aTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTtattat     <  1:1382497/150‑1 (MQ=255)
                                                                                                                                                     cACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGta  <  1:404257/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGTA  >  NC_000913/1853875‑1854173

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