Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,047 T→G 85.0% 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,0470TG85.0% 298.7 / 35.7 93F72V (TTT→GTT) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base G (42/37);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACT  >  NC_000913/2535901‑2536190
                                                                                                                                                  |                                                                                                                                               
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                               >  4:508725/1‑149 (MQ=255)
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                               <  2:460263/149‑1 (MQ=255)
tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                               <  1:672769/149‑1 (MQ=255)
  aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGa                                                                                                                                             >  1:267785/1‑149 (MQ=255)
  aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGa                                                                                                                                             <  1:338656/149‑1 (MQ=255)
  aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGa                                                                                                                                             <  1:569211/149‑1 (MQ=255)
   ttGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGaa                                                                                                                                            >  2:258146/1‑149 (MQ=255)
    tGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGaaa                                                                                                                                           <  4:239065/149‑1 (MQ=255)
    tGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTCCAGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGaaa                                                                                                                                           >  2:634175/1‑149 (MQ=255)
    tGCTCCGCTCCCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGaaa                                                                                                                                           >  2:157610/1‑149 (MQ=255)
     gCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaac                                                                                                                                          >  2:227956/1‑149 (MQ=255)
       tCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaacca                                                                                                                                        >  4:659374/1‑149 (MQ=255)
         cGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaac                                                                                                                                      <  3:168399/149‑1 (MQ=255)
         cGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaac                                                                                                                                      >  3:136371/1‑149 (MQ=255)
         cGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaaccaac                                                                                                                                      <  2:253913/149‑1 (MQ=255)
          gCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAaaccaacc                                                                                                                                     <  4:431898/149‑1 (MQ=255)
            tctctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCAc                                                                                                                                   >  2:96871/1‑149 (MQ=255)
                    gagaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATtctc                                                                                                                           >  1:351722/1‑149 (MQ=255)
                                atCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGa                                                                                                               <  3:476756/149‑1 (MQ=255)
                                  cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                                                                                                             <  1:157610/149‑1 (MQ=255)
                                   gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAg                                                                                                            >  2:607133/1‑149 (MQ=255)
                                    aaGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgc                                                                                                           >  1:559878/1‑149 (MQ=255)
                                    aaGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgc                                                                                                           <  3:13713/149‑1 (MQ=255)
                                      gACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcgg                                                                                                         >  2:792493/1‑149 (MQ=255)
                                       aCGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcggc                                                                                                        <  4:156727/149‑1 (MQ=255)
                                              ggATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAAc                                                                                                 <  1:325171/149‑1 (MQ=255)
                                               gATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACt                                                                                                >  1:590062/1‑149 (MQ=255)
                                                   tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc                                                                                            >  3:219892/1‑149 (MQ=255)
                                                     gTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGt                                                                                          >  3:342615/1‑149 (MQ=255)
                                                     gTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGt                                                                                          <  4:438395/149‑1 (MQ=255)
                                                      tttttGCGGCAAAAATCGTTAGACATGGTATTGCTAGCAAGCCAGCGGGTGATAAAATGGTCGCGCCAGTAGACGGCACTATTGGTAAAATCGTTGAAACCAACCACGCATCCTCCATCGAATGTGATAGCGGCGTTGAACTGTTCGTc                                                                                         >  4:143069/1‑149 (MQ=255)
                                                      tttttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATcagatt                                                                                                                  <  4:272872/124‑6 (MQ=255)
                                                      tttttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATcagatt                                                                                                                  >  3:272872/1‑119 (MQ=255)
                                                      tttttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTc                                                                                         >  2:516537/1‑149 (MQ=255)
                                                          tGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACt                                                                                     <  3:22931/149‑1 (MQ=255)
                                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                                                                                 >  3:62258/1‑149 (MQ=255)
                                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                                                                                 >  2:457243/1‑149 (MQ=255)
                                                              gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                                                                                 <  1:227956/149‑1 (MQ=255)
                                                                  aaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc                                                                             <  4:429929/149‑1 (MQ=255)
                                                                   aaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACTAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCg                                                                            >  2:73527/1‑149 (MQ=255)
                                                                    aTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGa                                                                           <  1:792493/149‑1 (MQ=255)
                                                                     tCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCAc                                                                                                                                   <  2:148938/92‑1 (MQ=255)
                                                                     tCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCAc                                                                                                                                   >  1:148938/1‑92 (MQ=255)
                                                                      cGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGaca                                                                         <  4:136371/149‑1 (MQ=255)
                                                                       gTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac                                                                        >  3:276750/1‑149 (MQ=255)
                                                                          ggTGATGGTATTGCTATCAAGCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGt                                                                     <  3:659024/149‑1 (MQ=255)
                                                                          ggTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGt                                                                     >  1:693234/1‑149 (MQ=255)
                                                                           gTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGtt                                                                    >  3:246542/1‑149 (MQ=255)
                                                                               tGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAAc                                                                <  2:680326/149‑1 (MQ=255)
                                                                                  tATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGa                                                             >  4:325773/1‑149 (MQ=255)
                                                                                   aTTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaa                                                            >  2:156565/1‑149 (MQ=255)
                                                                                    ttGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                           <  2:590062/149‑1 (MQ=255)
                                                                                    ttGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                           <  3:508725/149‑1 (MQ=255)
                                                                                     tGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAg                                                          <  4:276750/149‑1 (MQ=255)
                                                                                         aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGCTATCGACACCGTTGAACTGAAAGGCGa                                                      >  1:409274/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTTGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                               >  4:559583/1‑59 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTTGCGCCAGTAGACGGCACCATTGGTAAAATCgtt                                                                                                                                               <  3:559583/59‑1 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                     >  4:395268/1‑149 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                     <  3:395268/149‑1 (MQ=255)
                                                                                          tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                     >  2:707471/1‑149 (MQ=255)
                                                                                           cAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAg                                                    >  3:628675/1‑149 (MQ=255)
                                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGaca                                                                         >  4:87435/1‑127 (MQ=255)
                                                                                            aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGaca                                                                         <  3:87435/127‑1 (MQ=255)
                                                                                             aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  >  4:650990/1‑149 (MQ=255)
                                                                                             aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  >  3:532464/1‑149 (MQ=255)
                                                                                             aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  >  3:431331/1‑149 (MQ=255)
                                                                                             aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  >  1:292763/1‑149 (MQ=255)
                                                                                               ccAACGGGTATCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCtt                                                <  3:80814/149‑1 (MQ=255)
                                                                                                cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                               <  4:179576/149‑1 (MQ=255)
                                                                                                cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                               >  3:747777/1‑149 (MQ=255)
                                                                                                cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                               >  3:634/1‑149 (MQ=255)
                                                                                                 aaCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATc                                                                                                                       >  1:392290/1‑76 (MQ=255)
                                                                                                 aaCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATc                                                                                                                       <  2:392290/76‑1 (MQ=255)
                                                                                                  aCGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCaa                                             <  3:276903/149‑1 (MQ=255)
                                                                                                   cGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAg                                            <  3:680600/149‑1 (MQ=255)
                                                                                                      gTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGt                                         >  2:247190/1‑149 (MQ=255)
                                                                                                       tAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTa                                        >  2:386206/1‑149 (MQ=255)
                                                                                                        aaCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAt                                       <  4:634/149‑1 (MQ=255)
                                                                                                         aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                      >  3:671829/1‑149 (MQ=255)
                                                                                                           aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                     <  2:236889/132‑1 (MQ=255)
                                                                                                           aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                     >  1:236889/1‑132 (MQ=255)
                                                                                                            aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGACCTGGTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCt                                   >  2:581591/1‑149 (MQ=255)
                                                                                                               tGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaa                                <  2:351722/149‑1 (MQ=255)
                                                                                                               tGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCACTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaa                                <  4:62258/149‑1 (MQ=255)
                                                                                                                  tCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGTGTATTGCTgaagaa                             <  4:91030/149‑1 (MQ=255)
                                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  <  4:288653/126‑1 (MQ=255)
                                                                                                                    gcgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                  >  3:288653/1‑126 (MQ=255)
                                                                                                                     cgcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGt                          <  1:797281/149‑1 (MQ=255)
                                                                                                                      gcCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTc                         >  3:719246/1‑149 (MQ=255)
                                                                                                                           tAGACGGCACGATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGGCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAgcgc                    >  2:448643/1‑149 (MQ=255)
                                                                                                                             gACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGt                  <  1:156565/149‑1 (MQ=255)
                                                                                                                             gACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGt                  <  1:457243/149‑1 (MQ=255)
                                                                                                                                     cATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTg          <  1:386206/149‑1 (MQ=255)
                                                                                                                                      aTTGGTAAAATCGTTGAAACTAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg         <  3:659374/149‑1 (MQ=255)
                                                                                                                                      aTTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg         <  1:96871/149‑1 (MQ=255)
                                                                                                                                             aaaTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACt  <  3:325773/149‑1 (MQ=255)
                                                                                                                                              aaTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg         >  4:111453/1‑141 (MQ=255)
                                                                                                                                              aaTCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg         <  3:111453/141‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACT  >  NC_000913/2535901‑2536190

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: