Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,529 A→C 14.1% G302G (GGT→GGG)  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,967,5290AC14.1% 99.6 / 2.0 60G302G (GGT→GGG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (27/22);  new base C (7/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCC  >  NC_000913/2967381‑2967677
                                                                                                                                                    |                                                                                                                                                    
cgcCATCTCGCACGACAACACCGACTAACCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGa                                                                                                                                                      <  2:1158619/149‑1 (MQ=255)
  ccATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGAcc                                                                                                                                                    <  2:958463/149‑1 (MQ=255)
         gCACGACAACACCGACTAAGCGCTCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTGCATCTGCCACCAGAATGCAACGTTCCGGCGAGGCGTTCGGCCATTGATTAGCGTCATGAAGATGAATCCGCAATCGCTGGCCCAgcgcc                                                                                                                                             >  2:672563/1‑148 (MQ=255)
         gCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGCCCCAGCGCa                                                                                                                                             >  1:1054019/1‑149 (MQ=255)
         gCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCa                                                                                                                                             <  2:823617/149‑1 (MQ=255)
              acaacaCCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTaa                                                                                                                                        <  1:285149/149‑1 (MQ=255)
               caacaCCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTaaa                                                                                                                                       >  1:23908/1‑149 (MQ=255)
                   acCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGc                                                                                                                                   <  1:370397/149‑1 (MQ=255)
                      gACTACGCGATCCTGGGGCAGCGCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAg                                                                                                                                <  2:82940/149‑1 (MQ=255)
                       aCTAAGCGATCCTGGGGCAGCTCAGCAAGAGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGGTCCGGCCCGGCGTTCGGCCCTGGAGTAGCGTCATCACGATGCAACAGCCATCGCTGCCCCAGCGCACGTACATCGCCCGc                                                                                                                               >  1:698331/1‑149 (MQ=255)
                         tAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGCTTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGccc                                                                                                                             >  1:797257/1‑149 (MQ=255)
                           aGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGc                                                                                                                           >  2:742827/1‑149 (MQ=255)
                           aGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCACTTGATTAGCGTCATCAACATGAAACAGCAATCGCTGACCCCGCGCACGTAAATCGCCAGCCCGc                                                                                                                           >  2:207598/1‑149 (MQ=255)
                              gATCCTGGGGCCGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGAATAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCAGctct                                                                                                                        >  2:149444/1‑149 (MQ=255)
                                tCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCttt                                                                                                                      <  1:1077600/149‑1 (MQ=255)
                                     gggCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCATGCGTTCGGCCCTTGATTAGCGTCAACAAGATGAAACAGCAATCGCTGACCCAGCTCACGTAAATCGCCAGCACGCTCGTTgagtt                                                                                                                 >  2:97741/1‑147 (MQ=255)
                                         aGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGtt                                                                                                             <  1:828576/149‑1 (MQ=255)
                                              aGCACGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCCGCGCACGTAATTCGCCAGCCCGCTCTGTGAGATAGTTATCAc                                                                                                        >  1:23631/1‑149 (MQ=255)
                                                cAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATTGTTATCACtt                                                                                                      >  1:1071553/1‑149 (MQ=255)
                                                      ttGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCCTCAAGATGAAACAGCAATCGCTGACCCCGCGCACGTAAATCGCCAGCCCGGTCTTTGAGATAGTTATCACTTAGcgcg                                                                                                >  1:939978/1‑149 (MQ=255)
                                                                 aGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTc                                                                                     >  1:412866/1‑149 (MQ=255)
                                                                        cTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAt                                                                              >  2:336850/1‑149 (MQ=255)
                                                                         tGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAtt                                                                             >  1:613865/1‑149 (MQ=255)
                                                                            caccaGATTGAAACGTTACGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGTTGAAACAGCATTCGCTGACCAAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAttttt                                                                          <  2:752276/149‑1 (MQ=255)
                                                                            caccaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAATATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGTTAGTTATCACTTCGCGCGGCAAACTGTTCGGCAAAttttt                                                                          >  1:135153/1‑149 (MQ=255)
                                                                             accaGAATGAACCGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGGTATCACTTAGCGCGGCAAACTGGTCGGCAAATTTTTc                                                                         >  2:865182/1‑149 (MQ=255)
                                                                             accaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGCCCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCATACTGTTCGGCAAATTTTTc                                                                         >  1:1157667/1‑149 (MQ=255)
                                                                             accaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCGTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCATATTTTTc                                                                         >  1:1022798/1‑149 (MQ=255)
                                                                              ccaTAATGAAACGTTCCGGCCATGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCTCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCa                                                                        <  2:542161/149‑1 (MQ=255)
                                                                              ccaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCa                                                                        <  2:745197/149‑1 (MQ=255)
                                                                               caGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGAAAACGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTGAGCGCGGCGAACTGTTCGGCAAATTTTTCaa                                                                       >  1:1018091/1‑149 (MQ=255)
                                                                                 gAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCATTGAGATAGTTATAACTTATCGCGGCAAACTGTTCGGCAAATTTTTCAATg                                                                     >  1:475471/1‑149 (MQ=255)
                                                                                  aaTGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGCCCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTGATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGa                                                                    <  2:1054019/149‑1 (MQ=255)
                                                                                  aaTGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGa                                                                    <  1:22125/149‑1 (MQ=255)
                                                                                    tGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGCCCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGAcc                                                                  >  1:561740/1‑149 (MQ=255)
                                                                                    tGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGAcc                                                                  >  1:9700/1‑149 (MQ=255)
                                                                                               ggCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTATTCAATGACCGTTTTTACCGc                                                       >  1:278854/1‑149 (MQ=255)
                                                                                                 ccAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGccc                                                     <  1:5328/149‑1 (MQ=255)
                                                                                                   aGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGCCCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCAc                                                   >  1:562530/1‑149 (MQ=255)
                                                                                                           ggCCCTTGAGTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAc                                           <  1:1013022/149‑1 (MQ=255)
                                                                                                                ttGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc                                      <  1:152898/149‑1 (MQ=255)
                                                                                                                ttGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc                                      >  1:1159078/1‑149 (MQ=255)
                                                                                                                  gATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCt                                    <  1:406913/149‑1 (MQ=255)
                                                                                                                      cgCGGCATCAAGATGAAACAGCAATCGCGGCCCCAGCGCACGTTAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTCGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTAACGCCCACTCTGCCACCGAGCCTTTAt                                <  1:82067/148‑1 (MQ=255)
                                                                                                                      aGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGCGCCTTTAt                                >  2:1005809/1‑149 (MQ=255)
                                                                                                                      aGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTAt                                <  2:23908/149‑1 (MQ=255)
                                                                                                                        cGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCa                              >  2:372261/1‑149 (MQ=255)
                                                                                                                        cGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCa                              <  2:797257/149‑1 (MQ=255)
                                                                                                                         gTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCGACTATGAGACCGAGCCGTTATCaa                             >  2:930648/1‑149 (MQ=255)
                                                                                                                            atcaAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCt                          <  2:728444/149‑1 (MQ=255)
                                                                                                                              caAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCg                        <  1:515353/149‑1 (MQ=255)
                                                                                                                               aaGATGAAACAGCAATCGCTGCCCCAGCGCAAGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCTGCAAAATGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACACTGCCACCGAGCCTTAATCAACCTCgg                       >  2:299243/1‑149 (MQ=255)
                                                                                                                               aaGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCgg                       <  1:258147/149‑1 (MQ=255)
                                                                                                                                     aaaCAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACa                 >  2:1138000/1‑149 (MQ=255)
                                                                                                                                       aCAGCAATCGCTGCCCCACCGCACGTAAATCGCCAGCCCGCTCATTGAGATTGTTAACACAAAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCCTTTTTACAGCACACTCTGCCACCGAGCCTTTATCTACCTCGGCAAACAAt               >  2:541551/1‑149 (MQ=255)
                                                                                                                                        cAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACATCGGCAAACAAtt              >  1:959040/1‑149 (MQ=255)
                                                                                                                                          gCAATCGCTGACCCAGCGCTCGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTcg            <  2:539629/149‑1 (MQ=255)
                                                                                                                                          gCAAGCGCTGTCCCATCGCACGGAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGACAACTGTTCGGCAAATTTTTCCATGAACGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTcg            <  1:826929/149‑1 (MQ=255)
                                                                                                                                                gCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCACTGCCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGc      >  1:627517/1‑149 (MQ=255)
                                                                                                                                                gCTGACCCAGCGCACGGAAATCGCCAGCCCGCGCTTTGAGATAGTTATCACTTACCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGc      <  2:67213/149‑1 (MQ=255)
                                                                                                                                                  tGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCCACCTGTTCGGCAAAGTTTTCAATGCCCGTTTTTCCCGCCCCCTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGCCGCAg    >  2:984067/1‑149 (MQ=255)
                                                                                                                                                    aCCCAGCGCACGTAAATCGCCAGCCCGCTCTGTGAGATAGTTATCAATAAGCGCGGCAAACTGTTCGGAAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCATTTATCAACCTCGGCAAACAATTCGCGACGAAGcc  >  1:822656/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCC  >  NC_000913/2967381‑2967677

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