Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,063 C→A 22.0% A77E (GCA→GAA)  crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,0630CA22.0% 271.3 / 67.0 109A77E (GCA→GAA) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base A (13/11);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-02

Only 100 of 109 total aligned reads displayed.

TGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGC  >  NC_000913/2535917‑2536209
                                                                                                                                                  |                                                                                                                                                  
tGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCAt                                                                                                                                                  >  4:703125/1‑149 (MQ=255)
  gCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCAtt                                                                                                                                                 <  1:166878/148‑1 (MQ=255)
   cGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtc                                                                                                                                                >  1:855480/1‑148 (MQ=255)
   cGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtc                                                                                                                                                >  1:12625/1‑148 (MQ=255)
   cGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATtc                                                                                                                                                <  2:664330/148‑1 (MQ=255)
      gaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTctct                                                                                                                                             <  1:72616/148‑1 (MQ=255)
       aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTa                                                                                                                                            <  1:158283/148‑1 (MQ=255)
        tCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATc                                                                                                                                          <  4:353932/149‑1 (MQ=255)
          gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                                                                                                        >  3:559006/1‑149 (MQ=255)
           tCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                                                                                                        >  2:416059/1‑148 (MQ=255)
            cAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGaa                                                                                                                                       >  1:562007/1‑148 (MQ=255)
             aaTATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAAt                                                                                                                                      <  1:573980/148‑1 (MQ=255)
             aaTATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAAt                                                                                                                                      >  1:638028/1‑148 (MQ=255)
             aaTATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAAt                                                                                                                                      <  1:676063/148‑1 (MQ=255)
             aaTATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATc                                                                                                                                     <  3:787982/149‑1 (MQ=255)
             aaTATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAAt                                                                                                                                      >  1:667372/1‑148 (MQ=255)
              atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATc                                                                                                                                     <  2:908127/148‑1 (MQ=255)
              atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCt                                                                                                                                    >  3:884571/1‑149 (MQ=255)
              atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCt                                                                                                                                    >  4:50891/1‑149 (MQ=255)
                   gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAg                                                                                                                               >  4:707615/1‑149 (MQ=255)
                   gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAg                                                                                                                               >  3:125703/1‑149 (MQ=255)
                          tGCCGGATGTCGTTTTTGCTGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGt                                                                                                                         <  2:45165/148‑1 (MQ=255)
                          tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGt                                                                                                                         <  2:487362/148‑1 (MQ=255)
                          tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGt                                                                                                                         <  1:162934/148‑1 (MQ=255)
                           gCCTGATGTCGTTTTTGCGGATAAAATCGTTGGTGATGGTATTGCTATCAATCCAACGGGTAACAAAATGGTCGCGGCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTg                                                                                                                       >  4:102188/1‑149 (MQ=255)
                             cGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGaa                                                                                                                     <  3:125168/149‑1 (MQ=255)
                             cGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGa                                                                                                                      <  2:12625/148‑1 (MQ=255)
                                aTGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACt                                                                                                                   >  1:304725/1‑148 (MQ=255)
                                  gTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATATGATAGCGGCGTTGAACTGtt                                                                                                                >  4:691580/1‑149 (MQ=255)
                                   tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGtg                                                                                                                <  2:434077/148‑2 (MQ=255)
                                            cGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACtt                                                                                                       <  1:962392/148‑1 (MQ=255)
                                             ggAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg                                                                                                     >  4:534497/1‑149 (MQ=255)
                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                                                                                                    >  3:446731/1‑149 (MQ=255)
                                                  aaaGCGTTGGTGATAGTATTGCTATCAAACCAACGGGTAACAATATGGTCGCGCCAGTAGACGGCACCATTGTTAAACTCTTTGGAACCATCCACGCATTGTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc                                                                                                <  4:690310/149‑1 (MQ=255)
                                                     tCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGac                                                                                             >  3:486312/1‑149 (MQ=255)
                                                       gTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac                                                                                           >  3:213308/1‑149 (MQ=255)
                                                             gATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGa                                                                                     >  3:662514/1‑149 (MQ=255)
                                                              aTGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGa                                                                                     <  1:416059/148‑1 (MQ=255)
                                                              aTGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGaa                                                                                    <  4:239043/149‑1 (MQ=255)
                                                               tGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAAc                                                                                   <  4:213308/149‑1 (MQ=255)
                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGGATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAg                                                                             >  4:809082/1‑149 (MQ=255)
                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                                              >  1:668642/1‑148 (MQ=255)
                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                                              <  1:58812/148‑1 (MQ=255)
                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAg                                                                             <  3:72919/149‑1 (MQ=255)
                                                                         aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                                         >  3:332409/1‑149 (MQ=255)
                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                                        <  4:265586/149‑1 (MQ=255)
                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                                        <  3:50891/149‑1 (MQ=255)
                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                                         <  2:667372/148‑1 (MQ=255)
                                                                           cAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                                        >  2:510383/1‑148 (MQ=255)
                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAg                                                                       <  1:909563/148‑1 (MQ=255)
                                                                              aCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                                     >  1:113135/1‑148 (MQ=255)
                                                                                 aaCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGt                                                                                        >  2:433247/1‑126 (MQ=255)
                                                                                 aaCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGt                                                                                        <  1:433247/126‑1 (MQ=255)
                                                                                  aCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCa                                                                 <  2:638028/148‑1 (MQ=255)
                                                                                    gggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                                               <  1:510383/148‑1 (MQ=255)
                                                                                    gggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                                               >  1:341422/1‑148 (MQ=255)
                                                                                     ggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                                              >  1:570811/1‑148 (MQ=255)
                                                                                       tAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                                               >  2:927453/1‑145 (MQ=255)
                                                                                       tAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                                               <  1:927453/145‑1 (MQ=255)
                                                                                          taaacTGGCCGCCCCAGTAGTCATCACCGCTAGTAAACTCTATGAAACCAACCACGCATTCCCTCTCGAATCTGATACAGGCGTTTAACTGTTCATCCACCTCGATATCGACACCGTTGAACTGAAAGCCGAAGGCTTCAAGCGTAtt                                                         >  2:510440/2‑148 (MQ=255)
                                                                                          cAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                                         >  1:96093/1‑148 (MQ=255)
                                                                                          cAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                                         >  2:917719/1‑148 (MQ=255)
                                                                                          cAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTCGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTg                                                        >  3:218978/1‑149 (MQ=255)
                                                                                           aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGc                                                       >  4:180527/1‑149 (MQ=255)
                                                                                            aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGc                                                       >  1:154436/1‑148 (MQ=255)
                                                                                            aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCt                                                      <  4:662514/149‑1 (MQ=255)
                                                                                            aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTCGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCt                                                      <  4:218978/149‑1 (MQ=255)
                                                                                                 gTCGCGCCAGTGGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaaga                                                 >  3:556995/1‑149 (MQ=255)
                                                                                                  tCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaagaa                                                <  4:339020/149‑1 (MQ=255)
                                                                                                   cgcgCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                                               <  3:951627/149‑1 (MQ=255)
                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAgg                                              >  3:341784/1‑149 (MQ=255)
                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAgg                                              >  4:511397/1‑149 (MQ=255)
                                                                                                     cgcCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAgg                                              >  1:128343/1‑148 (MQ=255)
                                                                                                     cgcCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAgg                                              <  2:154436/148‑1 (MQ=255)
                                                                                                       ccAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATc                                                                                                >  1:101960/1‑96 (MQ=255)
                                                                                                       ccAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATc                                                                                                <  2:101960/96‑1 (MQ=255)
                                                                                                       ccAGTAGACGGCACCATTGGTAAAATCTATGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCa                                           <  3:314944/149‑1 (MQ=255)
                                                                                                            aGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAgcgc                                       >  2:214836/1‑148 (MQ=255)
                                                                                                               cGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGa                                   <  3:102188/149‑1 (MQ=255)
                                                                                                                 gCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaaa                                 <  4:332409/149‑1 (MQ=255)
                                                                                                                     cATTGGTAAAATCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTg                             <  4:58456/149‑1 (MQ=255)
                                                                                                                      aTTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGAGAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg                            <  3:703125/149‑1 (MQ=255)
                                                                                                                       ttGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg                            >  2:920020/1‑148 (MQ=255)
                                                                                                                           tAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATa                       >  4:665943/1‑149 (MQ=255)
                                                                                                                            aaaaTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATa                       >  2:443809/1‑148 (MQ=255)
                                                                                                                             aaaTCTTTGAAACCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATAc                      >  2:583605/1‑148 (MQ=255)
                                                                                                                                tCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGt                   >  2:300921/1‑148 (MQ=255)
                                                                                                                                     gAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCa                                                                 <  3:622550/97‑1 (MQ=255)
                                                                                                                                     gAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCa                                                                 >  4:622550/1‑97 (MQ=255)
                                                                                                                                      aaaCCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaa                                                                               >  3:80464/1‑82 (MQ=255)
                                                                                                                                      aaaCCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaa                                                                               <  4:80464/82‑1 (MQ=255)
                                                                                                                                      aaaCCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGaa            <  3:534497/149‑1 (MQ=255)
                                                                                                                                      aaaCCAACCACGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGa             <  1:722999/148‑1 (MQ=255)
                                                                                                                                       aaccaaccACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAt           <  4:765773/149‑1 (MQ=255)
                                                                                                                                        accaaccaCGCGTTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAtt          >  3:617661/1‑149 (MQ=255)
                                                                                                                                         ccaaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAtt          >  1:818492/1‑148 (MQ=255)
                                                                                                                                           aaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTGCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTg        >  2:745083/1‑148 (MQ=255)
                                                                                                                                           aaccaCGAATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTg        <  2:562007/148‑1 (MQ=255)
                                                                                                                                            accaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGAt      >  3:728835/1‑149 (MQ=255)
                                                                                                                                                 gCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGc  >  1:603936/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGC  >  NC_000913/2535917‑2536209

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: