Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,246 C→A 77.9% P138Q (CCG→CAG)  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,2460CA77.9% 321.0 / 105.8 149P138Q (CCG→CAG) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base A (54/62);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTC  >  NC_000913/2536100‑2536394
                                                                                                                                                  |                                                                                                                                                    
cGTCCACTTCGTTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAATAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCgg                                                                                                                                                    >  4:137100/1‑149 (MQ=255)
 gTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgt                                                                                                                                                   <  1:293681/149‑1 (MQ=255)
 gTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGgt                                                                                                                                                   <  3:540228/149‑1 (MQ=255)
  tCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCACGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGgtt                                                                                                                                                  <  1:1169861/149‑1 (MQ=255)
     aCTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGgttgtt                                                                                                                                               >  1:1498052/1‑149 (MQ=255)
     aCTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGgttgtt                                                                                                                                               <  3:551845/149‑1 (MQ=255)
       ttCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTAt                                                                                                                                             <  2:747329/149‑1 (MQ=255)
       ttCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTAt                                                                                                                                             >  3:291677/1‑149 (MQ=255)
       ttCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTAt                                                                                                                                             >  2:73909/1‑149 (MQ=255)
            tATCGACACCGGTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCa                                                                                                                                        >  2:714066/1‑149 (MQ=255)
              tCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAAc                                                                                                                                      >  1:1032118/1‑149 (MQ=255)
                gACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACAt                                                                                                                                    >  1:671351/1‑149 (MQ=255)
                gACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACAt                                                                                                                                    >  1:1428286/1‑149 (MQ=255)
                 acacCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATg                                                                                                                                   >  4:739965/1‑149 (MQ=255)
                 acacCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATg                                                                                                                                   <  1:64140/149‑1 (MQ=255)
                  cacCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATgg                                                                                                                                  <  3:478278/149‑1 (MQ=255)
                   acCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGa                                                                                                                                 >  3:682136/1‑149 (MQ=255)
                    ccGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGAc                                                                                                                                <  4:106905/149‑1 (MQ=255)
                    ccGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGAc                                                                                                                                <  2:666367/149‑1 (MQ=255)
                     cGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACg                                                                                                                               >  2:248365/1‑149 (MQ=255)
                       ttGAACTGAAAGGCGAAGGCTTCAAGTGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGaa                                                                                                                             >  2:289002/1‑149 (MQ=255)
                        tGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCTAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGaaa                                                                                                                            <  2:841363/149‑1 (MQ=255)
                        tGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGaaa                                                                                                                            <  2:1032118/149‑1 (MQ=255)
                             tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCaaa                                                                                                                       <  2:1028569/149‑1 (MQ=255)
                             tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCaaa                                                                                                                       <  2:534957/149‑1 (MQ=255)
                             tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCaaa                                                                                                                       >  3:539860/1‑149 (MQ=255)
                             tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCaaa                                                                                                                       <  3:140559/149‑1 (MQ=255)
                             tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCaaa                                                                                                                       <  2:1117824/149‑1 (MQ=255)
                                aaGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGaa                                                                                                                    <  2:1136352/149‑1 (MQ=255)
                                aaGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGaa                                                                                                                    <  3:1268262/149‑1 (MQ=255)
                                aaGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGaa                                                                                                                    <  3:739965/149‑1 (MQ=255)
                                  ggCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACt                                                                                                                  <  4:1348751/149‑1 (MQ=255)
                                   gCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTg                                                                                                                 <  4:336335/149‑1 (MQ=255)
                                    cGAAGGCTTCAAGCGTATTGCTTAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGa                                                                                                                <  4:1006493/149‑1 (MQ=255)
                                    cGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGa                                                                                                                <  2:145816/149‑1 (MQ=255)
                                     gAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGgtt                                                                                                                                                  <  1:287898/114‑1 (MQ=255)
                                      aaGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATc                                                                                                              <  1:73909/149‑1 (MQ=255)
                                         gCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCaaa                                                                                                           <  1:463397/149‑1 (MQ=255)
                                          cTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAAc                                                                                                          <  3:566777/149‑1 (MQ=255)
                                           ttCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACt                                                                                                         >  2:814142/1‑149 (MQ=255)
                                           ttCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAAc                                                                                                          <  2:943506/148‑1 (MQ=255)
                                           ttCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAAc                                                                                                          >  1:943506/1‑148 (MQ=255)
                                           ttCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACt                                                                                                         >  1:225597/1‑149 (MQ=255)
                                           ttCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAAAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACt                                                                                                         >  1:857402/1‑149 (MQ=255)
                                              aaGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTc                                                                                                      <  1:827797/149‑1 (MQ=255)
                                               aGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTcc                                                                                                     <  1:597507/149‑1 (MQ=255)
                                                     ttGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGc                                                                                               <  3:137100/149‑1 (MQ=255)
                                                      tGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCg                                                                                              <  4:595035/149‑1 (MQ=255)
                                                       gCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGt                                                                                             >  3:1150721/1‑149 (MQ=255)
                                                         tGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTaa                                                                                           <  3:690709/149‑1 (MQ=255)
                                                         tGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATATGCCGGTGCTGGAAGAGAAAGCCAAGCCTACCCTGACTAAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTaa                                                                                           >  4:316977/1‑149 (MQ=255)
                                                          gaagaaGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAAc                                                                                          >  3:1096547/1‑149 (MQ=255)
                                                           aagaagGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAAcc                                                                                         >  1:232473/1‑149 (MQ=255)
                                                             gaagGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGTCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGt                                                                                       >  2:513003/1‑149 (MQ=255)
                                                              aagGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTcc                                                                                                                                         <  3:727799/98‑1 (MQ=255)
                                                              aagGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTcc                                                                                                                                         >  4:727799/1‑98 (MQ=255)
                                                               agGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTgg                                                                                     <  3:702985/149‑1 (MQ=255)
                                                                 gTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTa                                                                                                 <  4:1285843/135‑1 (MQ=255)
                                                                 gTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTa                                                                                                 >  3:1285843/1‑135 (MQ=255)
                                                                     gcgcGTGAAAGTTGTCGATACTGTCATAGAATTTGATCCGCCGCTACTGCAACCGAAAGCCATGTCTACCCTGACTCAGGATGTTACCTTCAACATGGACGAAATCAAAGAACTGATCAAACGGTCCGGTAGCGTAACCGTGGGTTaac                                                                               <  2:528030/149‑2 (MQ=255)
                                                                     gcgcGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGaaa                                                                               <  3:810565/149‑1 (MQ=255)
                                                                         gTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAAcccc                                                                           <  3:1024472/149‑1 (MQ=255)
                                                                           gAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCgg                                                                         <  1:514307/149‑1 (MQ=255)
                                                                                 tGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATcc                                                                   <  1:459515/149‑1 (MQ=255)
                                                                                          tGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCaaga                                                          >  1:660563/1‑149 (MQ=255)
                                                                                           gTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGAGTGAAACCCCGGTTATCCGCATCaagaa                                                         <  4:1150721/149‑1 (MQ=255)
                                                                                            tCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCaag                                                           <  3:554866/146‑1 (MQ=255)
                                                                                            tCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCaag                                                           >  4:554866/1‑146 (MQ=255)
                                                                                              aTTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAATCCCGGTTATCCGCATCAAGAAGTa                                                      >  3:482720/1‑149 (MQ=255)
                                                                                              aTTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTa                                                      >  3:19772/1‑149 (MQ=255)
                                                                                               ttGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTaa                                                     >  2:1072568/1‑149 (MQ=255)
                                                                                               ttGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTaa                                                     >  3:262576/1‑149 (MQ=255)
                                                                                               ttGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTaa                                                     >  1:1175745/1‑149 (MQ=255)
                                                                                               ttGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTaa                                                     >  4:525132/1‑149 (MQ=255)
                                                                                                    tttGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAAccc                                                                            <  4:648982/121‑1 (MQ=255)
                                                                                                    tttGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAAccc                                                                            >  3:648982/1‑121 (MQ=255)
                                                                                                    tttGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCtt                                                >  1:1386429/1‑149 (MQ=255)
                                                                                                      tGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGc                                              <  4:463172/149‑1 (MQ=255)
                                                                                                      tGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGc                                              <  3:994649/149‑1 (MQ=255)
                                                                                                        aTCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGAGTGAAACCTCGGTTATCCGCATCAAGAAGTAATTCTTGCCg                                            >  4:887958/1‑149 (MQ=255)
                                                                                                         tCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGc                                           >  2:1373457/1‑149 (MQ=255)
                                                                                                            gCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGt                                        <  2:66499/149‑1 (MQ=255)
                                                                                                                  gctgGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGaaaaa                                  <  1:987911/149‑1 (MQ=255)
                                                                                                                   ctgGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAAt                                 <  3:931318/149‑1 (MQ=255)
                                                                                                                    tgGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATg                                <  2:1454145/149‑1 (MQ=255)
                                                                                                                    tgGAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATg                                >  2:166568/1‑149 (MQ=255)
                                                                                                                      gAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGgc                              >  4:1028106/1‑149 (MQ=255)
                                                                                                                      gAAGAGAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGgc                              >  2:844353/1‑149 (MQ=255)
                                                                                                                        agagAAAGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGgcgc                            <  3:579402/149‑1 (MQ=255)
                                                                                                                            aaaGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCAt                        >  1:693732/1‑149 (MQ=255)
                                                                                                                            aaaGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGTCGAAATCAATGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCAt                        >  1:731286/1‑149 (MQ=255)
                                                                                                                            aaaGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCaaga                                                          <  3:1444543/115‑1 (MQ=255)
                                                                                                                            aaaGCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCaaga                                                          >  4:1444543/1‑115 (MQ=255)
                                                                                                                               gCCAAGTCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCg                                                                          >  4:712543/1‑96 (MQ=255)
                                                                                                                                     tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt               <  1:1072568/149‑1 (MQ=255)
                                                                                                                                     tCTACCCTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt               <  4:707301/149‑1 (MQ=255)
                                                                                                                                         cccTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttt           >  2:1374049/1‑149 (MQ=255)
                                                                                                                                           cTGACTCAGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttttt         >  1:1497639/1‑149 (MQ=255)
                                                                                                                                           cTGACTCAGGTTGTTATCGCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGCAGCGTAACCGTGGGTGAAACCTCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCTCCCATCGGAGCCAtttttct         >  2:1142441/1‑147 (MQ=255)
                                                                                                                                                 ccGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCTGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCtt   >  3:1009899/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTC  >  NC_000913/2536100‑2536394

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: ‑