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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,854,0240CA 6.0% 315.8 / 8.2 109intergenic (‑45/+72)ydjE/ydjFputative MFS sugar transporter, membrane protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (49/53);  new base A (3/4);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 109 total aligned reads displayed.

CTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGT  >  NC_000913/1853877‑1854172
                                                                                                                                                   |                                                                                                                                                    
cTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                    >  1:385291/1‑150 (MQ=255)
cTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                    >  2:1181272/1‑150 (MQ=255)
 tGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGTGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACt                                                                                                                                                   >  2:498395/1‑150 (MQ=255)
   ttAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTg                                                                                                                                                 <  1:1617526/150‑1 (MQ=255)
    taataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCAATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGAGGCAAATTCACTTGc                                                                                                                                                >  1:1906203/1‑150 (MQ=255)
       taaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGc                                                                                                                                             >  2:809541/1‑150 (MQ=255)
        aaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCa                                                                                                                                            >  2:1897855/1‑150 (MQ=255)
         aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                    >  2:917442/1‑141 (MQ=255)
         aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                    <  1:917442/141‑1 (MQ=255)
          cacaCTAAAGCTTATAATACCAAATATACGATAAAGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                          <  1:521156/147‑3 (MQ=255)
          cAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCACAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                          <  1:528687/150‑3 (MQ=255)
           aGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaatt                                                                                                                                         <  2:1917048/150‑4 (MQ=255)
           aGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaatt                                                                                                                                         <  2:647539/150‑4 (MQ=255)
             aCTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCttt                                                                                                                                       >  1:128999/1‑150 (MQ=255)
             aCTAAAGCTATTAATAACAAATATCCGATAATGAAACCGGGCCAAAGTCAAGCGGTCCAGGCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCttt                                                                                                                                       <  2:290397/150‑1 (MQ=255)
              cTACAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGAGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGAAGCAATTGCGGCAAATTAACTTGCCGCATCATTa                                                                                                                                      >  1:1255442/1‑150 (MQ=255)
              cTAAAGCTTATAATACCAAATATACGATAATGAAACCGGTCCAAAGGCAAGCGGTCCAGACATCAGCCAATTTGATCATATAGTTCCATTAGGTGCCTCGCATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                      >  1:1255441/1‑150 (MQ=255)
              cTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                      <  2:1844552/150‑1 (MQ=255)
              atAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                      <  2:1006377/149‑1 (MQ=255)
                ataGCTTATAATACCATATATACGATTATGAATCCGGGGCAATGTCTAGCGGTCCAGTCTTGTGCGAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAATGGTAGCAATTGCGGCATATTCACTTGCCGCATCTTTatc                                                                                                                                    <  2:1631565/148‑1 (MQ=255)
                aaaGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTatc                                                                                                                                    <  2:1101555/150‑1 (MQ=255)
                      tataATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCAtt                                                                                                                              >  2:569391/1‑150 (MQ=255)
                       ataataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCAttt                                                                                                                             >  2:1842791/1‑150 (MQ=255)
                         aataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTa                                                                                                                           >  1:459175/1‑150 (MQ=255)
                         aataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTa                                                                                                                           <  1:1181272/150‑1 (MQ=255)
                            aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                        <  1:1526978/150‑1 (MQ=255)
                            aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                        <  2:1866691/150‑1 (MQ=255)
                             ccAAATATACGATAATGAAACCGGTCCAAATGCAAGCGGTCCAGTCTTGCGCCAATTTGATCTTATTGTTCCATTATGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTc                                                                                                                       <  2:643600/150‑1 (MQ=255)
                             ccAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTc                                                                                                                       >  1:1027582/1‑150 (MQ=255)
                              cAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCa                                                                                                                      >  2:1883342/1‑150 (MQ=255)
                              cAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCa                                                                                                                      >  1:926653/1‑150 (MQ=255)
                               aaaTATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCaa                                                                                                                     >  1:808449/1‑150 (MQ=255)
                                 atatatGATAATGAAACCTGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATAAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCCAgg                                                                                                                   <  1:159042/150‑1 (MQ=255)
                                  tataCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgc                                                                                                                  >  1:1804683/1‑150 (MQ=255)
                                  tataCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgc                                                                                                                  <  2:1586155/150‑1 (MQ=255)
                                    taCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCACCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgcgc                                                                                                                >  1:1457629/1‑150 (MQ=255)
                                     aCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCa                                                                                                               <  2:623117/150‑1 (MQ=255)
                                     aCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGGTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCa                                                                                                               <  1:379233/150‑1 (MQ=255)
                                      cGATAATGCAACCGGGCCAAAGGCAAGCGGTCCAGTCTTTCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAg                                                                                                              <  2:1941380/150‑1 (MQ=255)
                                      cGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGGCTCGGATTTGCTTAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCCTCTTGATCATCATTTTATTTCAAGGCGCAg                                                                                                              >  2:1269403/1‑150 (MQ=255)
                                       gATTATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAgg                                                                                                             <  2:128999/150‑1 (MQ=255)
                                          aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAAc                                                                                                          <  2:1377556/150‑1 (MQ=255)
                                              aaaCAGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAATGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaat                                                                                                      <  1:254064/150‑1 (MQ=255)
                                                   gggCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAatat                                                                                                 <  2:411985/150‑1 (MQ=255)
                                                    ggCCAAATGCAAGCGGTCCAGTCTTGCTCCAATTTGATCATATTGTTCCATTATGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATAtt                                                                                                <  2:368409/150‑1 (MQ=255)
                                                     gCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTg                                                                                               <  1:1316752/150‑1 (MQ=255)
                                                        aaaGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACa                                                                                            >  2:73997/1‑150 (MQ=255)
                                                             cAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgtg                                                                                       >  2:744062/1‑150 (MQ=255)
                                                              aaGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTATATTTCAAGGCGCAGGAACAAATAATATTGACAACgtgt                                                                                      <  1:1887598/150‑1 (MQ=255)
                                                               aGCGGTCCAGTCTTGCGCCAATTTGATCTTATTGTTCCATTAGGTGCCTCGGATTTGCTAAATAGGTAGCAATTGCGTCAAATTCACTTTCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                     <  2:1118814/150‑1 (MQ=255)
                                                               aGCGGTCCAGTCTTGCGCCAATTTGATCATATTGATCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                     <  1:1840272/150‑1 (MQ=255)
                                                                  ggTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCAGAATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAAAGTGTGCtt                                                                                  >  2:359542/1‑150 (MQ=255)
                                                                    tCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTaa                                                                                >  1:1766829/1‑150 (MQ=255)
                                                                    tCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTaa                                                                                >  1:526472/1‑150 (MQ=255)
                                                                      cAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAAt                                                                              >  1:1981054/1‑150 (MQ=255)
                                                                         tCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTAc                                                                           <  1:1641059/150‑1 (MQ=255)
                                                                          attGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTAcc                                                                          <  2:807115/149‑1 (MQ=255)
                                                                            tGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACcac                                                                        <  2:546718/150‑1 (MQ=255)
                                                                              cgcCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCacac                                                                      >  2:319277/1‑150 (MQ=255)
                                                                                 cAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                   <  1:498395/150‑1 (MQ=255)
                                                                                 cAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                   <  1:1121604/150‑1 (MQ=255)
                                                                                   aTTTGATCATATTGTTCCATTTGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCt                                                                 >  1:1431643/1‑150 (MQ=255)
                                                                                    tttGATCATATTTTTCCATTAGGTGCCTAGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTg                                                                <  2:292170/150‑1 (MQ=255)
                                                                                    tttGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTg                                                                <  1:1598355/150‑1 (MQ=255)
                                                                                      tGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGa                                                              <  2:1667667/150‑1 (MQ=255)
                                                                                         tCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGatat                                                           <  1:1061712/150‑1 (MQ=255)
                                                                                                gTTCCATTAGGTGCGTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACAACACCAGCTGGATATTATTGTc                                                    >  1:1394708/1‑150 (MQ=255)
                                                                                                gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTATATTTCAAGGCGCAGGAACAAATAATATTGACACCGTGTGCTTAAATTACCACACCAGCGGGATATTATTGtg                                                    >  2:143402/1‑149 (MQ=255)
                                                                                                    cATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTg                                                <  2:1838750/150‑1 (MQ=255)
                                                                                                     aTTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCTGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGc                                               <  1:809541/150‑1 (MQ=255)
                                                                                                       tAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCaa                                             >  1:2007046/1‑150 (MQ=255)
                                                                                                        aGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCAAACCAGCTGGATATTATTGTCTTTGCaaa                                            >  2:1003881/1‑150 (MQ=255)
                                                                                                          gTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAAt                                          <  2:794545/150‑1 (MQ=255)
                                                                                                           ttccTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATg                                         <  1:1450484/148‑1 (MQ=255)
                                                                                                           tGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATg                                         <  1:1601154/150‑1 (MQ=255)
                                                                                                              cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                      >  1:1305917/1‑150 (MQ=255)
                                                                                                                 ggATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATcc                                   <  1:1883342/150‑1 (MQ=255)
                                                                                                                 ggATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATcc                                   <  1:1764668/150‑1 (MQ=255)
                                                                                                                   aTTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCAc                                 >  1:945898/1‑150 (MQ=255)
                                                                                                                        cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                            <  2:297580/150‑1 (MQ=255)
                                                                                                                        cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                            >  2:848427/1‑150 (MQ=255)
                                                                                                                           aaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATTTTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAcc                         <  1:1808488/150‑1 (MQ=255)
                                                                                                                           aaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTAACATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAcc                         >  1:397206/1‑150 (MQ=255)
                                                                                                                              aGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCgg                      <  1:1589188/150‑1 (MQ=255)
                                                                                                                               ggTAGCAATTGCGGCAAATTCACTTGCCGGATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGt                     >  1:1228178/1‑150 (MQ=255)
                                                                                                                               ggTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGt                     <  1:2003074/150‑1 (MQ=255)
                                                                                                                                 tAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGa                   <  2:808449/150‑1 (MQ=255)
                                                                                                                                  aGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                  >  1:737870/1‑150 (MQ=255)
                                                                                                                                    cAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAttt                <  2:1457629/150‑1 (MQ=255)
                                                                                                                                     aaTTGCGGCAAATTCACTTGCCGCATCTTTATCATAATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTa               >  2:1679939/1‑150 (MQ=255)
                                                                                                                                        tGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa            >  1:987758/1‑150 (MQ=255)
                                                                                                                                        tGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa            >  2:149365/1‑150 (MQ=255)
                                                                                                                                        tGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa            <  1:992046/150‑1 (MQ=255)
                                                                                                                                         gCGGCAAATTCAATTGCCGCATATTTATAATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACAAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTTTCAg           <  1:1214060/150‑1 (MQ=255)
                                                                                                                                          cGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCTGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGt          >  2:1501994/1‑150 (MQ=255)
                                                                                                                                           ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt         >  1:1248829/1‑150 (MQ=255)
                                                                                                                                           ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt         >  1:1889410/1‑150 (MQ=255)
                                                                                                                                            gCAAATTCACTTGCCGCATCTTTATCATAATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTTTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtta        >  1:1899209/1‑150 (MQ=255)
                                                                                                                                                 ttCACTTGCCGCATCTGTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCCCCCGGTGATTTATCAGTTATTATg   >  2:356830/1‑150 (MQ=255)
                                                                                                                                                  tCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGt  <  1:744062/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGT  >  NC_000913/1853877‑1854172

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: