Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,527,603 A→G L70P (CTA→CCA)  sgcC ← putative PTS system enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,527,6030AG83.2% 692.4 / 143.6 322L70P (CTA→CCA) sgcCputative PTS system enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (149/119);  total (176/146)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 324 total aligned reads displayed.

AAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCA  >  NC_000913/4527513‑4527695
                                                                                          |                                                                                            
aaaaTAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTggg                                                                                            >  1:2810817/1‑93 (MQ=255)
 aaaTAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGc                                                                                           >  2:6681338/1‑93 (MQ=255)
    tAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGa                                                                                        <  1:1764539/93‑1 (MQ=255)
       cTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAAc                                                                                     >  1:734920/1‑93 (MQ=255)
       cTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAAc                                                                                     >  2:6151523/1‑93 (MQ=255)
       cTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAAc                                                                                     <  1:1024600/93‑1 (MQ=255)
       cTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAAc                                                                                     <  2:5365702/93‑1 (MQ=255)
        tGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACg                                                                                    >  1:3112931/1‑93 (MQ=255)
          ggATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCt                                                                                  >  2:4896749/1‑93 (MQ=255)
             tAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCaa                                                                               >  2:4604091/1‑93 (MQ=255)
              aaTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAat                                                                              <  2:2345918/93‑1 (MQ=255)
              aaTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAat                                                                              <  2:1237595/93‑1 (MQ=255)
               aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAatc                                                                             >  1:3339938/1‑93 (MQ=255)
               aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAatc                                                                             >  2:1559924/1‑93 (MQ=255)
                tCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAatca                                                                            <  2:2889549/93‑1 (MQ=255)
                tCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAatca                                                                            >  1:5467950/1‑93 (MQ=255)
                tCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAatca                                                                            <  1:3213585/93‑1 (MQ=255)
                  aTGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAatcatc                                                                          >  2:2026639/1‑93 (MQ=255)
                    ggCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCAc                                                                        <  2:5669778/93‑1 (MQ=255)
                     gcgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACt                                                                       >  1:5864129/1‑93 (MQ=255)
                      cgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACtt                                                                      <  1:6885410/93‑1 (MQ=255)
                       gcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACttt                                                                     >  2:6910687/1‑93 (MQ=255)
                        ccAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTg                                                                    >  2:748064/1‑93 (MQ=255)
                        ccAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTg                                                                    >  1:6569949/1‑93 (MQ=255)
                         cAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGa                                                                   >  1:3805958/1‑93 (MQ=255)
                         cAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGa                                                                   <  1:1290325/93‑1 (MQ=255)
                         cAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGa                                                                   >  2:4125135/1‑93 (MQ=255)
                         cAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGa                                                                   <  1:4553787/93‑1 (MQ=255)
                          aaTGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGAt                                                                  >  1:2924874/1‑93 (MQ=255)
                           aTGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATa                                                                 <  1:3281962/93‑1 (MQ=255)
                            tGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAg                                                                >  2:3174816/1‑93 (MQ=255)
                            tGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAg                                                                >  2:2679855/1‑93 (MQ=255)
                            tGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAg                                                                >  2:4801918/1‑93 (MQ=255)
                             ggCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAgg                                                               <  1:1584525/93‑1 (MQ=255)
                             ggCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAgg                                                               <  1:2083153/93‑1 (MQ=255)
                              gCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAggg                                                              <  1:4042641/93‑1 (MQ=255)
                               cTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAgggg                                                             <  1:5914972/93‑1 (MQ=255)
                               cTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAgggg                                                             >  2:7117273/1‑93 (MQ=255)
                                 gTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGa                                                           <  1:6251992/93‑1 (MQ=255)
                                     cATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGa                                                       >  1:4794890/1‑93 (MQ=255)
                                     cATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGa                                                       <  2:5759826/93‑1 (MQ=255)
                                     cATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGa                                                       <  1:926330/93‑1 (MQ=255)
                                       tAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAgg                                                     >  2:3842062/1‑93 (MQ=255)
                                       tAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAgg                                                     >  1:4310052/1‑93 (MQ=255)
                                        aaCCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGc                                                    >  2:4182600/1‑93 (MQ=255)
                                         aCCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCt                                                   >  2:5460075/1‑93 (MQ=255)
                                          ccAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTg                                                  >  1:2573105/1‑93 (MQ=255)
                                          ccAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTg                                                  <  1:7267081/93‑1 (MQ=255)
                                           cAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGt                                                 >  1:4230467/1‑93 (MQ=255)
                                            aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                <  2:3058145/93‑1 (MQ=255)
                                            aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                <  1:7125035/93‑1 (MQ=255)
                                            aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                >  2:5260753/1‑93 (MQ=255)
                                            aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                <  1:2220613/93‑1 (MQ=255)
                                            aaCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTc                                                >  2:894171/1‑93 (MQ=255)
                                             aCGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCg                                               >  1:3550573/1‑93 (MQ=255)
                                               gCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGAt                                             <  1:6776830/93‑1 (MQ=255)
                                                      gCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATg                                      >  1:2438724/1‑93 (MQ=255)
                                                       cAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGa                                     <  1:6671057/93‑1 (MQ=255)
                                                        aGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGAc                                    >  2:2411634/1‑93 (MQ=255)
                                                         gggCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGGCTGAGGCTGTCGATCGCCATGACg                                   >  2:139423/1‑93 (MQ=255)
                                                          ggCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGa                                  <  2:5176431/93‑1 (MQ=255)
                                                          ggCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGa                                  <  2:4470463/93‑1 (MQ=255)
                                                          ggCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGa                                  >  2:7031974/1‑93 (MQ=255)
                                                           gCCTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGAt                                 <  1:5317710/93‑1 (MQ=255)
                                                           gCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGAt                                 <  2:2208141/93‑1 (MQ=255)
                                                             cTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCa                               <  1:4948686/93‑1 (MQ=255)
                                                             cTGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCa                               >  2:6748882/1‑93 (MQ=255)
                                                              tGCGCCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAc                              >  2:7365740/1‑93 (MQ=255)
                                                               gcgcCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCAcc                             <  2:5073283/93‑1 (MQ=255)
                                                                cgcCACCGTAAAAGACAGGCGGGGCGGGGCCAAACCGCTAAATCATCCTTTTGAGGGGGGGACGGGGGCTGTCGACCGCCATGACGATCACCa                            <  1:5110568/93‑1 (MQ=255)
                                                                cgcCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCa                            <  2:600627/93‑1 (MQ=255)
                                                                cgcCAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCa                            <  1:2296434/93‑1 (MQ=255)
                                                                 gcCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAg                           <  1:6391910/93‑1 (MQ=255)
                                                                  ccAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGt                          >  1:1938007/1‑93 (MQ=255)
                                                                  ccAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGt                          >  2:2562784/1‑93 (MQ=255)
                                                                  ccAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGt                          >  1:3995572/1‑93 (MQ=255)
                                                                  ccAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGt                          <  1:1232251/93‑1 (MQ=255)
                                                                   cAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTc                         >  2:6276736/1‑93 (MQ=255)
                                                                   cAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTc                         >  2:1753937/1‑93 (MQ=255)
                                                                   cAACGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTc                         <  1:2642920/93‑1 (MQ=255)
                                                                   cAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTc                         >  1:4848257/1‑93 (MQ=255)
                                                                   cAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTc                         <  2:4071670/93‑1 (MQ=255)
                                                                    aaCGTCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTcc                        >  2:3370135/1‑93 (MQ=255)
                                                                        tCAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATa                    <  1:4096603/93‑1 (MQ=255)
                                                                        tCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATa                    <  1:6075873/93‑1 (MQ=255)
                                                                        tCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATa                    <  2:7079651/93‑1 (MQ=255)
                                                                           aaGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCg                 >  2:3003228/1‑93 (MQ=255)
                                                                           aaGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCg                 <  1:1535917/93‑1 (MQ=255)
                                                                           aaGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCg                 >  1:6472081/1‑93 (MQ=255)
                                                                             gACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGAc               >  2:1896172/1‑93 (MQ=255)
                                                                             gACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGAc               <  1:3633102/93‑1 (MQ=255)
                                                                               cATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACaa             >  2:1651866/1‑93 (MQ=255)
                                                                                 tGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACaaaa           <  2:2786900/93‑1 (MQ=255)
                                                                                  gCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAAc          <  1:5048958/93‑1 (MQ=255)
                                                                                  gCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAAc          <  2:7167417/93‑1 (MQ=255)
                                                                                   cAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAAcc         <  2:944273/93‑1 (MQ=255)
                                                                                     gggCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCaa       <  1:352307/93‑1 (MQ=255)
                                                                                     gggCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCaa       >  1:113473/1‑93 (MQ=255)
                                                                                       gCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATc     <  2:2073188/93‑1 (MQ=255)
                                                                                        cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATcc    >  1:1366656/1‑93 (MQ=255)
                                                                                          |                                                                                            
AAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCA  >  NC_000913/4527513‑4527695

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