Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,968,180 T→G I85I (ATA→ATC ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,968,1800TG98.5% 471.6 / ‑3.0 136I85I (ATA→ATCptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (59/75);  total (60/76)
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.16e-01

Only 100 of 136 total aligned reads displayed.

GTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGC  >  NC_000913/2968083‑2968275
                                                                                                 |                                                                                               
 ttGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATc                                                                                              <  1:1144453/100‑1 (MQ=255)
 ttGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATc                                                                                              <  1:3055908/100‑1 (MQ=255)
 ttGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATc                                                                                              <  1:165661/100‑1 (MQ=255)
  tGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCg                                                                                             <  2:3293349/100‑1 (MQ=255)
  tGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCg                                                                                             >  1:1531385/1‑100 (MQ=255)
    aataatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGt                                                                                           >  2:764365/1‑100 (MQ=255)
    aataatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGt                                                                                           >  2:2492317/1‑100 (MQ=255)
      taatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTc                                                                                         <  1:1993090/100‑1 (MQ=255)
       aatTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTcc                                                                                        <  1:1853517/100‑1 (MQ=255)
            gTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAATCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAg                                                                                   >  2:533819/1‑100 (MQ=255)
            gTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAg                                                                                   >  1:4060174/1‑100 (MQ=255)
             tACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGc                                                                                  >  1:3217811/1‑100 (MQ=255)
               cGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTAATTGAAGCTGGGGTGCTTTTCCGCATCTGCAGGGTTGATCGGTTCCGCCAGCCt                                                                                >  2:1917110/1‑100 (MQ=255)
                gCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTg                                                                               >  1:3165769/1‑100 (MQ=255)
                gCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTg                                                                               >  2:3757553/1‑100 (MQ=255)
                 ccTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGc                                                                              <  2:1041589/100‑1 (MQ=255)
                  cTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGcc                                                                             >  2:649377/1‑100 (MQ=255)
                   tAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCa                                                                            <  1:1731398/100‑1 (MQ=255)
                   tAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCa                                                                            <  2:1544814/100‑1 (MQ=255)
                   tAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCa                                                                            >  2:2343737/1‑100 (MQ=255)
                       aaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCCAGGTTGATCGGTTCCGCCAGCCTGCCAACCa                                                                        >  2:3110373/1‑100 (MQ=255)
                       aaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCa                                                                        >  1:1068941/1‑100 (MQ=255)
                       aaCGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCa                                                                        <  1:2705789/100‑1 (MQ=255)
                         cgcgcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAgg                                                                      >  1:1853106/1‑100 (MQ=255)
                          gcgcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGc                                                                     >  1:1294138/1‑100 (MQ=255)
                            gcgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCcg                                                                   >  2:1363673/1‑100 (MQ=255)
                             cgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCcga                                                                  >  1:1306377/1‑100 (MQ=255)
                             cgGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCcga                                                                  >  1:32696/1‑100 (MQ=255)
                              ggAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCcgac                                                                 <  1:3979349/100‑1 (MQ=255)
                               gAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCcgacg                                                                <  2:1209420/100‑1 (MQ=255)
                                  aCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATc                                                             <  2:1920736/100‑1 (MQ=255)
                                  aCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATc                                                             <  2:688542/100‑1 (MQ=255)
                                  aCGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATc                                                             <  1:1523439/100‑1 (MQ=255)
                                   cGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATcc                                                            >  1:17887/1‑100 (MQ=255)
                                   cGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATcc                                                            <  2:787430/100‑1 (MQ=255)
                                   cGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATcc                                                            <  1:2988525/100‑1 (MQ=255)
                                      tcttctTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCACCCAGGCCGACGATCCCtt                                                         >  1:3827205/1‑100 (MQ=255)
                                      tcttctTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCtt                                                         <  2:3561686/100‑1 (MQ=255)
                                        ttcttTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTtca                                                       <  1:2908758/100‑1 (MQ=255)
                                         tcttTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTtcat                                                      >  1:3322806/1‑100 (MQ=255)
                                            tttACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCaa                                                   >  1:2692497/1‑100 (MQ=255)
                                             ttACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCaaa                                                  <  1:13486/100‑1 (MQ=255)
                                              tACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCaaan                                                 <  1:1650598/100‑2 (MQ=255)
                                               aCGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcg                                                >  2:1675386/1‑100 (MQ=255)
                                                cGGAGGGGATTTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGGTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgc                                               >  2:902467/1‑100 (MQ=255)
                                                cGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgc                                               >  1:2484109/1‑100 (MQ=255)
                                                cGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgc                                               >  1:1993029/1‑100 (MQ=255)
                                                cGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgc                                               >  1:3789849/1‑100 (MQ=255)
                                                 ggAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAAcgcg                                              <  2:3552736/100‑1 (MQ=255)
                                                   aGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAg                                            <  2:1324089/100‑1 (MQ=255)
                                                   aGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAg                                            <  2:3089798/100‑1 (MQ=255)
                                                    ggggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGc                                           >  2:2866715/1‑100 (MQ=255)
                                                    ggggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGc                                           <  2:213025/100‑1 (MQ=255)
                                                     gggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCTAGCg                                          <  1:2684812/100‑1 (MQ=255)
                                                     gggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCg                                          <  1:680915/100‑1 (MQ=255)
                                                        aTGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTa                                       <  1:74816/100‑1 (MQ=255)
                                                         tGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTAc                                      <  2:1304701/100‑1 (MQ=255)
                                                         ggTATTTGAAGCGGGGGTGCTTTTGCGCATCAGCAAGTTGGATAGGAGCCCCCAGACTGCCTACCAGGCCAACACTCCCTCCATCACACGCAAGCggtgt                                      >  2:3959760/2‑96 (MQ=255)
                                                          gTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACa                                     >  2:3703872/1‑100 (MQ=255)
                                                            aTTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGt                                   >  2:3606143/1‑100 (MQ=255)
                                                             tttGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTg                                  >  2:151842/1‑100 (MQ=255)
                                                             tttGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTg                                  <  2:232397/100‑1 (MQ=255)
                                                               tGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCg                                <  2:2838989/100‑1 (MQ=255)
                                                                gAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGa                               >  2:2454493/1‑100 (MQ=255)
                                                                  aGCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGAcc                             >  2:1637634/1‑100 (MQ=255)
                                                                   gCTGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcg                            <  1:2023768/100‑1 (MQ=255)
                                                                     tGGGGTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACcgcg                          >  2:3627409/1‑100 (MQ=255)
                                                                      ggggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGt                         >  1:1242335/1‑100 (MQ=255)
                                                                       gggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTg                        <  1:3500175/100‑1 (MQ=255)
                                                                       gggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTg                        >  1:535450/1‑100 (MQ=255)
                                                                        ggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTgg                       <  2:190002/100‑1 (MQ=255)
                                                                        ggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTgg                       <  1:570629/100‑1 (MQ=255)
                                                                        ggTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCAGCCAGCCTGCCAACCAGGACGACGATCCCTTCATCAAACACGAGCGTTACAGTGCGACCGCGTgg                       >  1:3790502/1‑100 (MQ=255)
                                                                         gTGCTTTTGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGt                      >  1:2325775/1‑100 (MQ=255)
                                                                            cTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGtttt                   >  2:3428867/1‑100 (MQ=255)
                                                                             ttttGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGttttt                  <  2:786508/100‑1 (MQ=255)
                                                                             ttttGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGttttt                  <  2:3180136/100‑1 (MQ=255)
                                                                             ttttGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGttttt                  <  1:1848900/100‑1 (MQ=255)
                                                                               ttGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGGTTTTTc                >  2:1388965/1‑100 (MQ=255)
                                                                               ttGCGCATCTGCAAGGTTGATCGGGTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTc                <  1:1942306/100‑1 (MQ=255)
                                                                                tGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCa               <  1:1119639/100‑1 (MQ=255)
                                                                                tGCGCATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCa               <  1:4002092/100‑1 (MQ=255)
                                                                                accgcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCa               <  2:1722447/98‑1 (MQ=255)
                                                                                 gcgcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAg              <  2:2626374/100‑1 (MQ=255)
                                                                                 gcgcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAg              >  2:3881861/1‑100 (MQ=255)
                                                                                  cgcATCTGGAAGGTTGGTCGGTTGTGCCAGCCTGCCCACCAGGCCGACGATCCCCTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGc             <  1:3309306/100‑1 (MQ=255)
                                                                                  cgcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGc             >  1:4025989/1‑100 (MQ=255)
                                                                                  cgcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGc             <  1:499646/100‑1 (MQ=255)
                                                                                   gcATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGcc            <  1:2399763/100‑1 (MQ=255)
                                                                                    cATCTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccc           <  2:839911/100‑1 (MQ=255)
                                                                                    cATCTGCAAGGTTGATCGGTGCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTCACAGTGCGACCGCGTGGTCTTTTCAGccc           >  2:3576207/1‑100 (MQ=255)
                                                                                       cTGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCg        <  1:2513441/100‑1 (MQ=255)
                                                                                        tGCAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGTTCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCgg       <  2:1025372/100‑1 (MQ=255)
                                                                                        tGCAAGGTTGATAGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGTGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCgg       >  2:3849260/1‑100 (MQ=255)
                                                                                          cAAGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGt     <  1:3889692/100‑1 (MQ=255)
                                                                                           aaGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTc    <  2:2934527/100‑1 (MQ=255)
                                                                                           aaGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTc    <  2:3293014/100‑1 (MQ=255)
                                                                                            aGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCg   >  2:139876/1‑100 (MQ=255)
                                                                                            aGGTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCg   <  1:2402423/100‑1 (MQ=255)
                                                                                             ggTTGATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGc  >  2:1069938/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGC  >  NC_000913/2968083‑2968275

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