Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,247,747 T→C T198T (ACA→ACG ycgC ← putative PTS system enzyme I/IIA/HPr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,247,7470TC99.5% 626.6 / ‑5.0 200T198T (ACA→ACGycgCputative PTS system enzyme I/IIA/HPr
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (82/117);  total (83/117)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

Only 100 of 211 total aligned reads displayed.

TGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAG  >  NC_000913/1247657‑1247836
                                                                                          |                                                                                         
tGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGgcgt                                                                                         <  1:7417862/93‑1 (MQ=255)
tGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGgcgt                                                                                         <  1:3183431/93‑1 (MQ=255)
 gCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGcgtg                                                                                        <  2:1115322/93‑1 (MQ=255)
  ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGa                                                                                       <  1:6106575/93‑1 (MQ=255)
  ccTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGa                                                                                       >  1:2616932/1‑93 (MQ=255)
   ctctTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGAc                                                                                      <  2:1495742/93‑1 (MQ=255)
   ctctTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGAc                                                                                      <  2:72770/93‑1 (MQ=255)
       tCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAt                                                                                  <  1:1272303/93‑1 (MQ=255)
       tCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAt                                                                                  >  2:506653/1‑93 (MQ=255)
        cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAtt                                                                                 <  2:5933473/93‑1 (MQ=255)
        cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAtt                                                                                 <  1:5163014/93‑1 (MQ=255)
        cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAtt                                                                                 >  1:4143590/1‑93 (MQ=255)
        cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAtt                                                                                 <  2:7125633/93‑1 (MQ=255)
         aGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAttt                                                                                >  2:5541211/1‑93 (MQ=255)
         aGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCAttt                                                                                <  2:5144992/93‑1 (MQ=255)
          gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTg                                                                               >  2:1916603/1‑93 (MQ=255)
           cTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGc                                                                              >  1:2502029/1‑93 (MQ=255)
           cTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGc                                                                              >  2:5764310/1‑93 (MQ=255)
             tCTTGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCGTTTGCCg                                                                            >  2:7592958/1‑93 (MQ=255)
             tCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCg                                                                            <  1:2364758/93‑1 (MQ=255)
                 gCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGtttt                                                                        <  1:6169156/93‑1 (MQ=255)
                 gCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGtttt                                                                        <  1:3256566/93‑1 (MQ=255)
                   ccTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGtttttt                                                                      <  1:699787/93‑1 (MQ=255)
                   ccTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGtttttt                                                                      >  2:7179017/1‑93 (MQ=255)
                     tttCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTc                                                                    >  1:6342181/1‑93 (MQ=255)
                          cAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCa                                                               <  2:7694846/93‑1 (MQ=255)
                          cAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCa                                                               >  2:5685791/1‑93 (MQ=255)
                              cAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACat                                                           <  1:5415961/93‑1 (MQ=255)
                               aGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACata                                                          <  1:6032715/93‑1 (MQ=255)
                                ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACatat                                                         >  1:2913880/1‑93 (MQ=255)
                                 gcgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATc                                                        <  1:4066361/93‑1 (MQ=255)
                                 gcgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATc                                                        <  1:3718756/93‑1 (MQ=255)
                                   gcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCgg                                                      >  2:113941/1‑93 (MQ=255)
                                      gCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCAt                                                   >  2:1970209/1‑93 (MQ=255)
                                        gTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTa                                                 >  2:5994598/1‑93 (MQ=255)
                                         tATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                <  1:7930311/93‑1 (MQ=255)
                                         tATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                >  1:6212368/1‑93 (MQ=255)
                                         tATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATAATCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                >  2:2233927/1‑93 (MQ=255)
                                          aTCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaaa                                               >  2:2923012/1‑93 (MQ=255)
                                            cGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATg                                             >  2:4277/1‑93 (MQ=255)
                                             gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGt                                            <  1:2441695/93‑1 (MQ=255)
                                              ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                           <  1:1939074/93‑1 (MQ=255)
                                              ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                           <  2:1200448/93‑1 (MQ=255)
                                               tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCg                                          <  1:6427960/93‑1 (MQ=255)
                                                ataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGa                                         <  2:378325/93‑1 (MQ=255)
                                                ataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGa                                         >  2:6719901/1‑93 (MQ=255)
                                                  aGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATa                                       <  1:6964629/93‑1 (MQ=255)
                                                  aGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATa                                       <  1:4453569/93‑1 (MQ=255)
                                                  aGCGAATTGGAGTAACCGAATCTGGTTAATCATTTCTTGGGTGGCGGATTTTCCGGTTTTTTCCAGCAACATATCGGCATTAAATGTCGATaa                                      <  1:760184/93‑1 (MQ=255)
                                                   gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATaa                                      <  1:5888985/93‑1 (MQ=255)
                                                   gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATaa                                      <  2:3939979/93‑1 (MQ=255)
                                                    cGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAg                                     >  1:2792410/1‑93 (MQ=255)
                                                     gAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                    >  2:2674883/1‑93 (MQ=255)
                                                      aaCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGt                                   <  2:5475842/93‑1 (MQ=255)
                                                       aCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGta                                  <  2:4641826/93‑1 (MQ=255)
                                                       aCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGta                                  <  2:3210915/93‑1 (MQ=255)
                                                       aCTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGta                                  <  2:4847289/93‑1 (MQ=255)
                                                        cTTGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                 >  2:3330336/1‑93 (MQ=255)
                                                         ttGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtata                                <  2:7092601/93‑1 (MQ=255)
                                                         ttGTAGTAACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtata                                <  1:4571729/93‑1 (MQ=255)
                                                            tagtaACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAAc                             <  2:6580383/93‑1 (MQ=255)
                                                             agtaACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAAcc                            <  1:6591916/93‑1 (MQ=255)
                                                              gtaACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCa                           <  2:4737858/93‑1 (MQ=255)
                                                              gtaACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCa                           <  1:3559880/93‑1 (MQ=255)
                                                               taACGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAg                          >  2:4824224/1‑93 (MQ=255)
                                                                aaCGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGc                         <  1:3345692/93‑1 (MQ=255)
                                                                 aCGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGcc                        <  2:6597653/93‑1 (MQ=255)
                                                                 aCGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGcc                        >  1:1024432/1‑93 (MQ=255)
                                                                  cGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCg                       <  1:959147/93‑1 (MQ=255)
                                                                  cGCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCg                       <  1:5788979/93‑1 (MQ=255)
                                                                   gCAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCgg                      <  1:1800961/93‑1 (MQ=255)
                                                                    cAATCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                     >  2:6894155/1‑93 (MQ=255)
                                                                    cAATCTGGTTAATACCTTTTGGCGTGACCCATTTGCGGTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                    <  2:2375204/93‑1 (MQ=255)
                                                                     aaTCTGGTTAATTCTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                    <  2:2711954/93‑1 (MQ=255)
                                                                     aaTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                    >  1:997096/1‑93 (MQ=255)
                                                                     aaTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGa                    <  1:3957352/93‑1 (MQ=255)
                                                                      aTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                   <  1:114459/93‑1 (MQ=255)
                                                                      aTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                   <  1:2044575/93‑1 (MQ=255)
                                                                      aTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                   <  1:732579/93‑1 (MQ=255)
                                                                      aTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                   >  1:7829953/1‑93 (MQ=255)
                                                                      aTCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                   >  2:1749899/1‑93 (MQ=255)
                                                                       tCTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                  <  2:915302/93‑1 (MQ=255)
                                                                       tCTGGTTAATACTCTCTGGCGGGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAgg                  <  1:4901341/93‑1 (MQ=255)
                                                                        cTGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGc                 >  2:4986016/1‑93 (MQ=255)
                                                                         tGGTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGcc                <  1:7899493/93‑1 (MQ=255)
                                                                          ggTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCg               >  1:3386847/1‑93 (MQ=255)
                                                                           gTTAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCgg              <  1:485714/93‑1 (MQ=255)
                                                                            ttAATGCTCTCTGGCGTAACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa             <  2:3262760/93‑1 (MQ=255)
                                                                            ttAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa             <  2:3092194/93‑1 (MQ=255)
                                                                            ttAATACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa             <  1:2035459/93‑1 (MQ=255)
                                                                               aTACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGt          >  1:3127250/1‑93 (MQ=255)
                                                                                tcctctctGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTa         <  2:6990619/91‑1 (MQ=255)
                                                                                tACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTa         <  1:384926/93‑1 (MQ=255)
                                                                                tACTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTa         >  1:1596842/1‑93 (MQ=255)
                                                                                 aCTCTCTGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc        >  2:6239986/1‑93 (MQ=255)
                                                                                  ctctctGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCTATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa       <  2:7831620/93‑1 (MQ=255)
                                                                                  ctctctGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa       <  1:805584/93‑1 (MQ=255)
                                                                                   tctctGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACAt      >  1:3303669/1‑93 (MQ=255)
                                                                                      ctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCa   >  2:2468864/1‑93 (MQ=255)
                                                                                      ctGGCGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCa   <  1:1206371/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAG  >  NC_000913/1247657‑1247836

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