Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,047 T→G 100% F72V (TTT→GTT)  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,0470TG97.2% 593.2 / 3.4 161F72V (TTT→GTT) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (81/75);  total (84/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 168 total aligned reads displayed.

TCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGT  >  NC_000913/2535901‑2536192
                                                                                                                                                  |                                                                                                                                                 
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                                 <  3:1487171/149‑1 (MQ=255)
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                                 <  2:1326737/149‑1 (MQ=255)
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                                 <  4:1094333/149‑1 (MQ=255)
  aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGa                                                                                                                                               <  2:146703/149‑1 (MQ=255)
    tGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGaaa                                                                                                                                             <  4:1199662/149‑1 (MQ=255)
      cTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaacc                                                                                                                                           >  2:427393/1‑149 (MQ=255)
      cTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaacca                                                                                                                                          >  1:1214637/1‑149 (MQ=255)
            tctctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCAc                                                                                                                                     >  3:567444/1‑149 (MQ=255)
              tctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGc                                                                                                                                   >  1:407200/1‑149 (MQ=255)
               ctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCa                                                                                                                                  <  3:729790/149‑1 (MQ=255)
                  gCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATtc                                                                                                                               <  3:707781/149‑1 (MQ=255)
                      gaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTa                                                                                                                           <  3:134429/149‑1 (MQ=255)
                       aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTAt                                                                                                                          <  3:1386458/149‑1 (MQ=255)
                       aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCCGAAAAAATCATTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCAAt                                                                                                                          <  4:515515/149‑1 (MQ=255)
                        tCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATc                                                                                                                         <  1:427393/149‑1 (MQ=255)
                          gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCTGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGa                                                                                                                       >  2:1029158/1‑149 (MQ=255)
                          gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTa                                                                                                                           >  4:617394/1‑145 (MQ=255)
                          gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGa                                                                                                                       >  3:571009/1‑149 (MQ=255)
                            cAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAAt                                                                                                                     >  4:383743/1‑149 (MQ=255)
                                 tCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGAt                                                                                                                >  4:1529012/1‑149 (MQ=255)
                                 tCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGAt                                                                                                                >  4:230514/1‑149 (MQ=255)
                                        cGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcggcg                                                                                                         <  4:1057062/149‑1 (MQ=255)
                                          tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                                                                                       <  1:802799/149‑1 (MQ=255)
                                          tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                                                                                       <  3:230514/149‑1 (MQ=255)
                                               gATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACt                                                                                                  >  1:159942/1‑149 (MQ=255)
                                                   tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc                                                                                              >  1:1105664/1‑149 (MQ=255)
                                                     gTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGt                                                                                            >  4:107194/1‑149 (MQ=255)
                                                     gTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGt                                                                                            >  2:852046/1‑149 (MQ=255)
                                                      tttttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTc                                                                                           >  1:1202146/1‑149 (MQ=255)
                                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTcc                                                                                          >  1:488701/1‑142 (MQ=255)
                                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                                                                                   >  1:341704/1‑149 (MQ=255)
                                                               aaaaaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGt                                                                                  >  1:679681/1‑149 (MQ=255)
                                                                ataaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTa                                                                                 <  2:461502/147‑1 (MQ=255)
                                                                      cGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGaca                                                                           <  2:378295/149‑1 (MQ=255)
                                                                              aTGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg                                                                                    >  3:916672/1‑132 (MQ=255)
                                                                              aTGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg                                                                                    <  4:916672/132‑1 (MQ=255)
                                                                                  tATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGa                                                               <  1:1014243/149‑1 (MQ=255)
                                                                                    ttGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                             <  1:999622/149‑1 (MQ=255)
                                                                                    ttGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                             <  2:1305722/149‑1 (MQ=255)
                                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAg                                                            <  4:588694/149‑1 (MQ=255)
                                                                                       cTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGc                                                          <  4:567444/149‑1 (MQ=255)
                                                                                       cTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGc                                                          >  1:240242/1‑149 (MQ=255)
                                                                                        tATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCg                                                         >  2:1532611/1‑149 (MQ=255)
                                                                                        tATCAAACCAACGGGTAACAAAATGGTCGCACCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGc                                                          >  4:7142/1‑148 (MQ=255)
                                                                                         aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATc                                                                                                                    <  2:974329/89‑1 (MQ=255)
                                                                                         aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATc                                                                                                                    >  1:974329/1‑89 (MQ=255)
                                                                                         aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                        >  1:1277101/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCGGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       >  3:1408405/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       >  2:556348/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       >  3:1551405/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       <  4:839806/149‑1 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       >  3:840116/1‑149 (MQ=255)
                                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc                                                                               <  4:1488695/123‑1 (MQ=255)
                                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc                                                                               >  3:1488695/1‑123 (MQ=255)
                                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAgg                                                     >  2:1214221/1‑149 (MQ=255)
                                                                                             aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                    >  4:673999/1‑149 (MQ=255)
                                                                                               ccAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCtt                                                  <  3:611060/149‑1 (MQ=255)
                                                                                                cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACt                                                                                       <  1:846000/111‑1 (MQ=255)
                                                                                                cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                                 >  1:1451960/1‑149 (MQ=255)
                                                                                                   cGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                              <  3:819343/149‑1 (MQ=255)
                                                                                                    gggTAACACAATGGTCGCGCCAGTAGACGGCACCCTTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCTACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                             <  2:141281/149‑1 (MQ=255)
                                                                                                    gggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                             >  3:1464018/1‑149 (MQ=255)
                                                                                                    gggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                             <  2:889077/149‑1 (MQ=255)
                                                                                                    gggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACAAACCTCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                             >  2:1155486/1‑149 (MQ=255)
                                                                                                      gTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCAt                                                                                                                                 >  3:256714/1‑63 (MQ=255)
                                                                                                      gTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCAt                                                                                                                                 <  4:256714/63‑1 (MQ=255)
                                                                                                      gTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcg                                                                                                            <  3:334949/84‑1 (MQ=255)
                                                                                                       tAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTa                                          <  3:1020243/149‑1 (MQ=255)
                                                                                                        aaCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAt                                         >  4:1303213/1‑149 (MQ=255)
                                                                                                         aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                                                                                       <  2:336176/86‑1 (MQ=255)
                                                                                                         aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                        >  3:884602/1‑149 (MQ=255)
                                                                                                         aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                        <  2:198537/149‑1 (MQ=255)
                                                                                                         aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                        >  4:490604/1‑149 (MQ=255)
                                                                                                           aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc                                                                               >  4:1364183/1‑108 (MQ=255)
                                                                                                           aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTg                                                                     >  2:834828/1‑118 (MQ=255)
                                                                                                            aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCt                                     >  1:25470/1‑149 (MQ=255)
                                                                                                             aaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTg                                    >  3:1161162/1‑149 (MQ=255)
                                                                                                                   cgcgCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                              <  1:675974/149‑1 (MQ=255)
                                                                                                                   cgcgCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                              <  3:1031999/149‑1 (MQ=255)
                                                                                                                   cgcgCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                              <  1:1158001/149‑1 (MQ=255)
                                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                              >  4:1144412/1‑148 (MQ=255)
                                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                              <  3:1144412/148‑1 (MQ=255)
                                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAgg                             >  2:1163870/1‑149 (MQ=255)
                                                                                                                           tAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGaa                                                                                                    >  1:93131/1‑71 (MQ=255)
                                                                                                                             gACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGt                    <  4:1408405/149‑1 (MQ=255)
                                                                                                                             gACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGt                    <  3:673999/149‑1 (MQ=255)
                                                                                                                              aCGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTc                                                                                     >  3:1200852/1‑83 (MQ=255)
                                                                                                                                ggCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaa                 >  4:725502/1‑149 (MQ=255)
                                                                                                                                 gCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaaa                <  4:173223/149‑1 (MQ=255)
                                                                                                                                     cATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTg            >  2:547628/1‑149 (MQ=255)
                                                                                                                                      aTTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAATGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGCCAGCGCGTGAAAGTTgg           >  4:1324181/1‑149 (MQ=255)
                                                                                                                                      aTTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg           >  4:20670/1‑149 (MQ=255)
                                                                                                                                      aTTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg           <  2:1187872/149‑1 (MQ=255)
                                                                                                                                       ttGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGc          <  2:240242/149‑1 (MQ=255)
                                                                                                                                       ttGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGc          <  1:852046/149‑1 (MQ=255)
                                                                                                                                         ggTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                       >  3:130666/1‑102 (MQ=255)
                                                                                                                                            aaaaTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATAc     <  3:812765/149‑1 (MQ=255)
                                                                                                                                             aaaTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACt    >  2:1550456/1‑149 (MQ=255)
                                                                                                                                               aTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGc                                             <  1:1463431/106‑1 (MQ=255)
                                                                                                                                               aTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGt  <  4:1226796/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGT  >  NC_000913/2535901‑2536192

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