Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,854,024 C→A 5.2% 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 5.2% 313.1 / 5.3 106intergenic (‑45/+72)ydjE/ydjFputative MFS sugar transporter, membrane protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (55/45);  new base A (3/3);  total (58/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 109 total aligned reads displayed.

CCCTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGTA  >  NC_000913/1853875‑1854173
                                                                                                                                                     |                                                                                                                                                     
cccTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTc                                                                                                                                                       >  2:661342/1‑150 (MQ=255)
  cTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCAc                                                                                                                                                     >  1:1492069/1‑150 (MQ=255)
      taataaCAGACTAACGCTTATAATACCAAATATACGATAATGAAACCGCGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTAGGATTTGGTATAAAGGTAGCATTTGCGGCAAATTCACTTGc                                                                                                                                                 >  2:444741/1‑150 (MQ=255)
      taataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGc                                                                                                                                                 >  1:817639/1‑150 (MQ=255)
       aataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGcc                                                                                                                                                <  2:560212/150‑1 (MQ=255)
       aataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGcc                                                                                                                                                >  1:710858/1‑150 (MQ=255)
       aataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGcc                                                                                                                                                >  2:1320612/1‑150 (MQ=255)
          aaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCa                                                                                                                                             <  2:801355/150‑1 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCAt                                                                                                                                            >  2:901876/1‑150 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCAt                                                                                                                                            >  1:726009/1‑150 (MQ=255)
            aaGACTAAAGCTAATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGCTGCCTCGGATTTGCTAAAAAGGGAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                           <  1:449735/149‑3 (MQ=255)
             aGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaatt                                                                                                                                          <  1:661445/150‑4 (MQ=255)
               aCTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaattgc                                                                                                                                        <  2:377574/150‑6 (MQ=255)
                atAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                       <  2:705873/149‑1 (MQ=255)
                 aaaaGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTat                                                                                                                                      <  1:772440/149‑1 (MQ=255)
                  aaaGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTatc                                                                                                                                     >  1:1465437/1‑150 (MQ=255)
                  aaaGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTata                                                                                                                                     >  2:617291/1‑149 (MQ=255)
                      cTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTAtcatca                                                                                                                                 <  2:217933/150‑1 (MQ=255)
                      cTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTAtcatca                                                                                                                                 <  1:1232957/150‑1 (MQ=255)
                      cTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTAtcatca                                                                                                                                 >  2:1000194/1‑150 (MQ=255)
                        tataATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATAATGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCAtt                                                                                                                               <  1:1517596/150‑1 (MQ=255)
                          taataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCAtttt                                                                                                                             <  1:1313161/150‑1 (MQ=255)
                           aataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTa                                                                                                                            >  1:807497/1‑150 (MQ=255)
                           aataCCAAATATACGATAATGAAAACGGGCCAAAGGCATTCGGTCCAGTCTTGCGCCAATTTGATCTTATTGTTCCATTAGGTGCGTCGGATTTGCTAAAAAGGTAGCAATTGCTGCAAATTCACTTGCGGCATCTTTATAATCAttatt                                                                                                                            >  1:1570520/1‑147 (MQ=255)
                             taCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAtt                                                                                                                          <  1:810248/150‑1 (MQ=255)
                              aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                         <  2:667031/150‑1 (MQ=255)
                              aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                         >  2:1119094/1‑150 (MQ=255)
                                      taCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgcgc                                                                                                                 >  2:583970/1‑150 (MQ=255)
                                      taCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCAATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgcgc                                                                                                                 >  1:372194/1‑150 (MQ=255)
                                        cGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAg                                                                                                               >  1:1583367/1‑150 (MQ=255)
                                          aTAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGa                                                                                                             >  2:252812/1‑150 (MQ=255)
                                           tAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGaa                                                                                                            >  1:178667/1‑150 (MQ=255)
                                            aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGaaa                                                                                                           >  2:326537/1‑149 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          <  1:1161030/150‑1 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          >  1:1094498/1‑150 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          >  2:1326156/1‑150 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGACAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAACGTAGCAATTGCGGCATATTCACTTGCCGCATCTTTATCATCATTTTATTTAAAGGCGCAGGAACa                                                                                                          >  1:729408/1‑150 (MQ=255)
                                              tGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACaa                                                                                                         >  1:763096/1‑150 (MQ=255)
                                                aaaCCGGGCCAAAGGCAAGTGGGCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaat                                                                                                       <  2:252261/150‑1 (MQ=255)
                                                  aCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          >  1:1313864/1‑145 (MQ=255)
                                                  aCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaataa                                                                                                     >  2:1020067/1‑150 (MQ=255)
                                                  aCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaataa                                                                                                     >  2:855270/1‑150 (MQ=255)
                                                   ccGGGCCAAAGGCAAGCGGTACAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaataat                                                                                                    <  1:1102742/150‑1 (MQ=255)
                                                       gCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTg                                                                                                <  2:1465437/150‑1 (MQ=255)
                                                        ccAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGa                                                                                               >  1:620115/1‑150 (MQ=255)
                                                         cAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGAc                                                                                              <  2:1444054/150‑1 (MQ=255)
                                                                   cGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCt                                                                                    <  2:1492069/150‑1 (MQ=255)
                                                                    ggTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCtt                                                                                   <  2:710858/150‑1 (MQ=255)
                                                                    ggTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCtt                                                                                   <  2:1212597/150‑1 (MQ=255)
                                                                              tGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACcac                                                                         >  1:874873/1‑150 (MQ=255)
                                                                                cgcCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCacac                                                                       >  1:621887/1‑150 (MQ=255)
                                                                                 gcCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACAcc                                                                      >  2:788858/1‑150 (MQ=255)
                                                                                   cAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                    <  1:544498/150‑1 (MQ=255)
                                                                                    aaTTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGc                                                                   >  1:649289/1‑150 (MQ=255)
                                                                                        tGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGa                                                               >  1:1011136/1‑150 (MQ=255)
                                                                                          aTCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGata                                                             >  2:1407162/1‑150 (MQ=255)
                                                                                           tCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGatat                                                            <  2:1492294/150‑1 (MQ=255)
                                                                                                ttGTTCCATTAGGTGCCTCTGATTTGCTAAAAATGTATCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGAGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTg                                                       <  1:125683/150‑1 (MQ=255)
                                                                                                ttGTTCCATTAGGTGCCTCGGATTTACTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTg                                                       <  1:1607629/150‑1 (MQ=255)
                                                                                                atGTTCCATTAGGTGACTCGGATTTGCTAAAAAGGTAGCAAATGCGCCAAATTCACTTGCCGAATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATattattt                                                       >  2:1193119/2‑149 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     <  2:399531/150‑1 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     <  1:1020067/150‑1 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     <  1:855270/150‑1 (MQ=255)
                                                                                                   ttCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCt                                                    <  2:172611/150‑1 (MQ=255)
                                                                                                      cATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTTAATTACCACACCAGCTGGATATTATTGTCTTTg                                                 >  2:1156049/1‑150 (MQ=255)
                                                                                                           ggTGACTAGGAATTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAATTCAAGGCGCAGGAACAAATAATATTGACAACGTGAGCTTAAATTACAACACCAGCTGGATATAATTCTCTTTGAaaaa                                            >  2:1353764/1‑150 (MQ=255)
                                                                                                             tGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATg                                          <  1:326537/150‑1 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGc                                         <  2:1311408/150‑1 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCACGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCGCCCCATCGGGATATTATTGGCTTTGCAAAATGCTa                                       >  2:870808/1‑150 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       >  2:1569182/1‑150 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       >  2:1287823/1‑150 (MQ=255)
                                                                                                                cTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTa                                       >  2:1280874/1‑150 (MQ=255)
                                                                                                                     ttttGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCAc                                  <  2:51749/149‑1 (MQ=255)
                                                                                                                        ttcTAAAAAGGTAGCAATTGCGGAAAATTCACTTTCAGAATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCActct                               <  1:195883/148‑1 (MQ=255)
                                                                                                                        tGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCActct                               <  1:74338/150‑1 (MQ=255)
                                                                                                                        tGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCActct                               <  2:807497/150‑1 (MQ=255)
                                                                                                                        tGCTAAAAAGGTAGCAATTGCGGCAAAATCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACAACACAAGCTGGATATTATTGTCTTTGCAAAATGCTATCCActct                               >  1:1416905/1‑150 (MQ=255)
                                                                                                                         gCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACAACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTg                              >  2:783082/1‑150 (MQ=255)
                                                                                                                          cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                             >  2:96807/1‑150 (MQ=255)
                                                                                                                          cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATAGTCTTTGCAAAATGCTATCCACTCTGc                             >  2:87148/1‑150 (MQ=255)
                                                                                                                          cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGATTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                             >  1:923073/1‑150 (MQ=255)
                                                                                                                           tAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGcc                            >  2:79979/1‑149 (MQ=255)
                                                                                                                            aaaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAc                           <  2:562647/150‑1 (MQ=255)
                                                                                                                                 ggTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGt                      >  1:876064/1‑150 (MQ=255)
                                                                                                                                   tAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGTCAACGTGTGCTTAAATTACCACACCAGCTGGATAATATTGTCTTTGCAAAAGGCTATCCACTCTGCACCCGGTGa                    <  1:103856/150‑1 (MQ=255)
                                                                                                                                   ataaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGATATCCACTCTGCACCCGGTGa                    >  2:759891/8‑150 (MQ=255)
                                                                                                                                     gCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAtt                  >  2:304486/1‑150 (MQ=255)
                                                                                                                                       agatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCTCACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTa                >  1:519170/4‑150 (MQ=255)
                                                                                                                                          tGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa             <  2:1282728/150‑1 (MQ=255)
                                                                                                                                          tGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa             <  1:1119094/150‑1 (MQ=255)
                                                                                                                                           gCGGCAAGTGAATTTGCCGCATCATTATCATCATATTATTTCAATGCGCAGTAACAAATAATATTGACACCGTGTGCTTAAATTACAACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATACACTCTGCACACGTTGATTTATCAg            >  2:195125/1‑150 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTTTTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAGTTATCAg            >  2:645323/1‑150 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAg            >  2:374313/1‑150 (MQ=255)
                                                                                                                                             ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTCCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt          >  1:1630755/1‑150 (MQ=255)
                                                                                                                                               cAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttat        <  2:391629/150‑1 (MQ=255)
                                                                                                                                                aaaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatt       <  2:1583367/150‑1 (MQ=255)
                                                                                                                                                   ttCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTTTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATg    <  1:87234/150‑1 (MQ=255)
                                                                                                                                                   ttCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATg    >  2:1015732/1‑150 (MQ=255)
                                                                                                                                                     cACTTGCCGCATCTTTATCATCATTTTATTTGTAGGCGCAGGAACAATTAATATTGAGAACGTGTGCTTAAATTACCTCTCCTGCTGGTTATTATTGTCTTTGAAAAATGCTATCAACTATGCACCCTGTGATTTATCAGTTATTATGta  <  2:1506272/150‑1 (MQ=255)
                                                                                                                                                     cACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGta  <  1:901876/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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: