Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,700 C→T 6.1% intergenic (‑404/+19) psuK ← / ← fruA pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,7000CT 6.1% 255.3 / 3.9 82intergenic (‑404/+19)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATCAGCA  >  NC_000913/2259552‑2259849
                                                                                                                                                    |                                                                                                                                                     
taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGcc                                                                                                                                                      <  1:762342/150‑1 (MQ=255)
  atCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCATTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCt                                                                                                                                                    >  1:1081689/1‑150 (MQ=255)
    cTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATTTTACAGGACAGGAAATTTCTGCCCTGt                                                                                                                                                  >  1:193199/1‑150 (MQ=255)
     tGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTa                                                                                                                                                 <  2:130706/150‑1 (MQ=255)
      gCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTaa                                                                                                                                                >  1:1142673/1‑150 (MQ=255)
        ccTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAaca                                                                                                                                              <  2:363546/150‑1 (MQ=255)
         cTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacat                                                                                                                                             <  1:600729/150‑2 (MQ=255)
           tGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTTCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacattt                                                                                                                                           <  2:794324/150‑4 (MQ=39)
           tGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAacacac                                                                                                                                           <  1:495615/150‑1 (MQ=255)
           tGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAacacac                                                                                                                                           >  2:4821/1‑150 (MQ=255)
           tGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATAACATAAACAAACGGTTCATACTCTTTGCCATCATCAAATTTTACAGGACAGTAAATTAATGCCATTTAacacac                                                                                                                                           >  2:302695/1‑150 (MQ=255)
            gAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGTACAGGAAATTTCTGCCCTGTAACACAcc                                                                                                                                          <  2:537990/150‑1 (MQ=255)
              ggAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACAtt                                                                                                                                        >  1:489818/1‑150 (MQ=255)
              ggAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTAATGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACTAACCGGTAATTCTCTATGCGCTCTTCAAATGTTACAGGACAGGAAATTTATGAACTTTAAAACAACtt                                                                                                                                        >  1:1198258/1‑150 (MQ=255)
                  aTAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttat                                                                                                                                    <  1:709617/150‑1 (MQ=255)
                     aCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTAc                                                                                                                                 >  2:62348/1‑150 (MQ=255)
                     aCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTAc                                                                                                                                 <  1:790264/150‑1 (MQ=255)
                               ttGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                                                                                       <  1:270238/150‑1 (MQ=255)
                                tgatgaACGTTTGAAGTTTCTACTGAAAGATATAACCATTTAATTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGc                                                                                                                      <  1:173066/150‑1 (MQ=255)
                                tgatgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTTTAACACACCTTTTATTACGCTGCTTTCGc                                                                                                                      >  2:28795/1‑150 (MQ=255)
                                tgatgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCATGTAACACACCTTTTATTACGCTGCTTTCGc                                                                                                                      >  2:1000799/1‑150 (MQ=255)
                                 gatgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCt                                                                                                                     >  2:197225/1‑150 (MQ=255)
                                  atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTATGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCACTGTAACCCACCTTTTATTACGCTGCATTCGCTa                                                                                                                    >  1:756062/1‑150 (MQ=255)
                                   tgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTGTCGCTAc                                                                                                                   >  2:1117645/1‑150 (MQ=255)
                                     aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                                                                                 <  1:475596/150‑1 (MQ=255)
                                          ttGAAGTTTATACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCATTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                                            <  1:207788/150‑1 (MQ=255)
                                           tGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCa                                                                                                           >  1:692880/1‑150 (MQ=255)
                                                  tCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCgg                                                                                                    >  1:93623/1‑150 (MQ=255)
                                                   cTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTAATACGCTGCTTTCGCTACTGCGTCCACTTCCggc                                                                                                   >  2:1143618/1‑149 (MQ=255)
                                                    tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTATGCGCTCATCAAATGTTACAGGACAGGAAAATTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGACCACTTCCGGAc                                                                                                  >  2:1008460/1‑150 (MQ=255)
                                                     aCTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCGTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACg                                                                                                 >  2:352883/1‑150 (MQ=255)
                                                      cTGAAAGATAAAAACATTTAAATTCATAAAATTACATTCACAAACCCGTCATTCTCATTGCGCTCATCAAATGTTACAGGACACGAAATATCTGCCCTGTTACACACCTTTTAATACGCTGCTTTCGCTACTGCTTCCACTTCCGGACGt                                                                                                >  2:683429/1‑150 (MQ=255)
                                                        taaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGttt                                                                                              <  1:1301827/149‑1 (MQ=255)
                                                           atataataCCATTTAAATTCATAAAATTACATTCACATACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTATCACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                                                                                           <  1:1160281/148‑1 (MQ=255)
                                                             aTAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGa                                                                                         >  1:274729/1‑150 (MQ=255)
                                                                aatccATTTAAATTCATATAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACCGGACAGGAAATTTCTGCCCTGTAAGTCACCTTTTATTACGCTTATTTCGGTACTGCGTCCTGTTCCAGAGGTTTCAGGatgt                                                                                      <  2:924724/147‑4 (MQ=255)
                                                                aaaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                                                                                      >  1:909931/1‑150 (MQ=255)
                                                                      tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                                                                                <  2:1157699/150‑1 (MQ=255)
                                                                       ttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGc                                                                               <  2:1027080/150‑1 (MQ=255)
                                                                        tAAATTCATAAAATTTCATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGcc                                                                              >  1:1045442/1‑150 (MQ=255)
                                                                         aaaTTCTTAAATTTACATTCAAAAAGCGGTGATTCTCTTTGCGCTCATCATATGTTTCAGGACAGGAAATTTCTGCCTTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCa                                                                             <  2:641609/150‑1 (MQ=255)
                                                                          aaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                                            <  1:94160/150‑1 (MQ=255)
                                                                          aaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                                            >  1:1125394/1‑150 (MQ=255)
                                                                            ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAACTTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                                                                          <  1:86317/150‑1 (MQ=255)
                                                                              cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc                                                                        >  2:358373/1‑150 (MQ=255)
                                                                                 aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                                                                     <  1:925953/150‑1 (MQ=255)
                                                                                        cATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                                              >  1:797721/1‑150 (MQ=255)
                                                                                           tCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTAc                                                           >  2:583200/1‑150 (MQ=255)
                                                                                             acaAACCGTTGATTCTCATTGCGCGCATCATATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTAGCTACTGCGTCAACTTGCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCAGCAGCGTAcct                                                         <  2:878736/150‑2 (MQ=255)
                                                                                             acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCg                                                         >  2:339476/1‑150 (MQ=255)
                                                                                              caAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCgg                                                        <  2:907028/150‑1 (MQ=255)
                                                                                                  ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaat                                                    >  2:1040674/1‑150 (MQ=255)
                                                                                                    ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaataa                                                  <  1:1295800/150‑1 (MQ=255)
                                                                                                    ggTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGATGCTTTCGCTACTGCGTACACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCCGCTGTCCGGCaataa                                                  >  2:811308/1‑150 (MQ=255)
                                                                                                     gTCATTCTCTTTGCGCTCATGAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaataat                                                 <  2:756062/150‑1 (MQ=255)
                                                                                                       cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTg                                               <  2:699938/150‑1 (MQ=255)
                                                                                                         ttCTCTTTGCGCTCATCAAATGTTACAGGACAGGAACTTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCt                                             <  1:133264/150‑1 (MQ=255)
                                                                                                             ctTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATAGCTACCa                                         >  2:1028055/1‑150 (MQ=255)
                                                                                                                      catcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCa                                >  2:1279923/1‑150 (MQ=255)
                                                                                                                      catcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCa                                <  1:607605/150‑1 (MQ=255)
                                                                                                                      catcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCa                                >  2:1268137/1‑150 (MQ=255)
                                                                                                                          gtgtGTTACAGGGCAGAAATTTCCTGTCCTGTAACATACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATAc                            >  1:181464/4‑150 (MQ=39)
                                                                                                                           aaTGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACTCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATAcc                           <  1:674449/150‑1 (MQ=255)
                                                                                                                           aaTGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACTCAATAcc                           <  1:1127632/150‑1 (MQ=255)
                                                                                                                           aaTGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATAcc                           <  1:318691/150‑1 (MQ=255)
                                                                                                                              gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTTCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc                        >  2:639377/1‑150 (MQ=255)
                                                                                                                              gTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCGCTTCCGGACGTTTCAGGAAGGCATAGGCCACACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTACCCCAATACCGGc                        >  1:886805/1‑150 (MQ=255)
                                                                                                                              gTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCTGACTTTTCAGGAAGGCATAGGCCAAACACGACACCAGCGTAACGGCAATTAATGCTACCAGGTAACCAAATACCGGc                        >  1:887418/1‑150 (MQ=255)
                                                                                                                               ttACATGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCg                       <  1:447473/150‑1 (MQ=255)
                                                                                                                               ttACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCg                       <  1:250663/150‑1 (MQ=255)
                                                                                                                                   aggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACAAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAt                   >  1:1300976/1‑150 (MQ=255)
                                                                                                                                    ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATa                  >  1:477934/1‑150 (MQ=255)
                                                                                                                                    ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATa                  >  2:262939/1‑150 (MQ=255)
                                                                                                                                     gacaggaTATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTATCCCAATACCGGCGTAATAg                 <  1:302695/150‑1 (MQ=255)
                                                                                                                                     gacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAg                 <  1:4821/150‑1 (MQ=255)
                                                                                                                                     gacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAg                 <  2:909931/150‑1 (MQ=255)
                                                                                                                                         ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTACCCCAATACCGGCGTAATAGCGcc             >  2:963669/1‑150 (MQ=255)
                                                                                                                                         ggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGcc             <  1:798108/150‑1 (MQ=255)
                                                                                                                                             aTTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAggg         <  2:1302426/150‑1 (MQ=255)
                                                                                                                                              tttCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGa        >  2:1033225/1‑150 (MQ=255)
                                                                                                                                              tttCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGGAATAGCGCCAGGGa        >  2:1302769/1‑150 (MQ=255)
                                                                                                                                                tCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATc      <  2:692880/150‑1 (MQ=255)
                                                                                                                                                  tGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGTCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATcag    <  2:1142673/150‑1 (MQ=255)
                                                                                                                                                    cccTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATcagca  >  2:478884/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATCAGCA  >  NC_000913/2259552‑2259849

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