Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,071,928 T→C intergenic (‑217/+68) epd ← / ← yggC D‑erythrose 4‑phosphate dehydrogenase/putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,071,9280TC96.4% 687.1 / 10.1 224intergenic (‑217/+68)epd/yggCD‑erythrose 4‑phosphate dehydrogenase/putative kinase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base C (115/101);  total (119/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 231 total aligned reads displayed.

AAACGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATAC  >  NC_000913/3071837‑3072014
                                                                                            |                                                                                      
aaaCGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATc                                                                                       >  1:3798821/1‑93 (MQ=255)
 aaCGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCCCACTTTTAccc                                                                                      <  1:1618527/93‑1 (MQ=255)
 aaCGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTAccc                                                                                      <  2:311535/93‑1 (MQ=255)
 aaCGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTAccc                                                                                      >  2:1653248/1‑93 (MQ=255)
   cGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCTAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACcccc                                                                                    <  2:226983/93‑1 (MQ=255)
   cGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACcccc                                                                                    <  1:1247517/93‑1 (MQ=255)
    gCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCt                                                                                   >  1:1495222/1‑93 (MQ=255)
      ttCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCttt                                                                                 >  2:4325690/1‑93 (MQ=255)
      ttCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCttt                                                                                 >  2:3741007/1‑93 (MQ=255)
      ttCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCttt                                                                                 <  2:1642064/93‑1 (MQ=255)
         aGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTtgt                                                                              <  1:1078844/93‑1 (MQ=255)
          gCTAGGATAAGGGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAGGACATTTATCTAACTCACATACACCTTTTACCCCCTTTtgtg                                                                             >  2:4642211/1‑93 (MQ=255)
             aGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGa                                                                          >  1:358579/1‑93 (MQ=255)
              ggATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGaa                                                                         >  1:892653/1‑93 (MQ=255)
                aTAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGc                                                                       <  1:2871579/93‑1 (MQ=255)
                 tAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCt                                                                      <  1:2822243/93‑1 (MQ=255)
                 tAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCt                                                                      >  1:4344150/1‑93 (MQ=255)
                  aaGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCtt                                                                     <  1:3384533/93‑1 (MQ=255)
                  aaGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCtt                                                                     >  2:2916962/1‑93 (MQ=255)
                  aaGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGAATCACATCACACTTTTATCCCCTTTTGTGGGAAGCtt                                                                     >  2:2407865/1‑93 (MQ=255)
                   aGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCttt                                                                    >  1:3961375/1‑93 (MQ=255)
                   aGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCttt                                                                    >  1:3655341/1‑93 (MQ=255)
                     cGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTAt                                                                  >  1:3926889/1‑93 (MQ=255)
                     cGAAACGTGGAATAAAAGGAATGTTTGACCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTAt                                                                  <  2:4362785/93‑1 (MQ=255)
                      gAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTGTTGTGGGAAGCTTTAtt                                                                 >  1:4225349/1‑93 (MQ=255)
                       aaaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTc                                                                <  1:1931523/93‑1 (MQ=255)
                       aaaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTc                                                                >  1:453045/1‑93 (MQ=255)
                        aaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTcc                                                               >  1:3254902/1‑93 (MQ=255)
                         aCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCa                                                              >  2:4328407/1‑93 (MQ=255)
                         aCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCa                                                              <  2:2732060/93‑1 (MQ=255)
                          cGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAg                                                             >  2:4043844/1‑93 (MQ=255)
                           gTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAgg                                                            >  1:1948165/1‑93 (MQ=255)
                           gTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAgg                                                            <  2:603199/93‑1 (MQ=255)
                            tGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGc                                                           <  2:343958/93‑1 (MQ=255)
                             ggAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCt                                                          >  1:4607476/1‑93 (MQ=255)
                             ggAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGATTTATTCCAGGCt                                                          >  1:4405710/1‑93 (MQ=255)
                              gAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTg                                                         >  1:785/1‑93 (MQ=255)
                              gAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTg                                                         >  1:2910617/1‑93 (MQ=255)
                                 tAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCg                                                      >  2:902142/1‑93 (MQ=255)
                                   aaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGTTGGCGta                                                    >  1:988532/1‑93 (MQ=255)
                                   aaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGta                                                    >  1:2613895/1‑93 (MQ=255)
                                    aaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaa                                                   <  2:3384051/93‑1 (MQ=255)
                                     aGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaat                                                  >  1:1419950/1‑93 (MQ=255)
                                     aGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaat                                                  <  1:63033/93‑1 (MQ=255)
                                      ggAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaata                                                 <  2:2432924/93‑1 (MQ=255)
                                      ggAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaata                                                 >  1:2314506/1‑93 (MQ=255)
                                       gAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataa                                                >  1:2773508/1‑93 (MQ=255)
                                       gAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataa                                                <  2:1555190/93‑1 (MQ=255)
                                        aaTGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataat                                               <  2:2126514/93‑1 (MQ=255)
                                        aaTGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataat                                               <  2:2108787/93‑1 (MQ=255)
                                        aaTGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataat                                               >  1:569786/1‑93 (MQ=255)
                                         aTGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataata                                              <  1:3304716/93‑1 (MQ=255)
                                          tGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGtaataataa                                             >  2:2607200/1‑93 (MQ=255)
                                            tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAAcc                                           >  2:1179208/1‑93 (MQ=255)
                                            tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAAcc                                           >  1:248800/1‑93 (MQ=255)
                                            tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAAcc                                           >  2:4402283/1‑93 (MQ=255)
                                            tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTACAGGCTGGCGTAATACTAAcc                                           >  1:4096214/1‑93 (MQ=255)
                                             ttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATTATAAccc                                          <  2:4194219/93‑1 (MQ=255)
                                                  cAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACaa                                     <  2:2511778/93‑1 (MQ=255)
                                                  cAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACaa                                     <  1:2500498/93‑1 (MQ=255)
                                                   aGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAAt                                    >  2:3276772/1‑93 (MQ=255)
                                                   aGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGATTTATTCCAGGCTGGCGTAATAATAACCCTACAAt                                    >  1:2237679/1‑93 (MQ=255)
                                                      cGAAGAAGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAAc                                 >  2:2734765/1‑93 (MQ=255)
                                                       gaagaaGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACt                                <  2:1312461/93‑1 (MQ=255)
                                                       gaagaaGACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACt                                >  2:4474983/1‑93 (MQ=255)
                                                        aagaagACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTg                               <  2:1246587/93‑1 (MQ=255)
                                                        aagaagACATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTg                               <  1:1629001/93‑1 (MQ=255)
                                                         agaagaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTgg                              >  1:318425/1‑93 (MQ=255)
                                                         agaagaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTgg                              <  1:2390856/93‑1 (MQ=255)
                                                           aagaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGaa                            <  1:4695846/93‑1 (MQ=255)
                                                           aagaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGaa                            <  2:785/93‑1 (MQ=255)
                                                            agtcATTTATCTGACTCACATCACATTTTTACCTCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACTCTACAATAACTGGAAt                           >  1:4206710/4‑93 (MQ=255)
                                                            agaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAAt                           >  2:500206/1‑93 (MQ=255)
                                                             gaCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATa                          >  1:413391/1‑93 (MQ=255)
                                                              aCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaa                         >  1:3525479/1‑93 (MQ=255)
                                                              aCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaa                         <  1:495062/93‑1 (MQ=255)
                                                              aCATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaa                         >  1:703260/1‑93 (MQ=255)
                                                               cATTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaaa                        >  1:2369596/1‑93 (MQ=255)
                                                                aTTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAAt                       <  1:4506306/93‑1 (MQ=255)
                                                                aTTTATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAAt                       >  2:368463/1‑93 (MQ=255)
                                                                  ttATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTg                     >  1:3675513/1‑93 (MQ=255)
                                                                   tATCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGt                    >  2:3211856/1‑93 (MQ=255)
                                                                    aTCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTc                   <  1:3283880/93‑1 (MQ=255)
                                                                     tCTGACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCa                  >  1:188823/1‑93 (MQ=255)
                                                                        gACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCt               <  2:3786245/93‑1 (MQ=255)
                                                                        gACTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCt               <  2:1851129/93‑1 (MQ=255)
                                                                         gcgcacaacacaCTTTTGCCCCCTTTTGTGAGAAGCTTTATTACAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCtt              <  2:2814544/90‑1 (MQ=255)
                                                                          cTCACATCACACTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtg             <  1:2169792/93‑1 (MQ=255)
                                                                           tcacatcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgt            >  1:4167588/1‑93 (MQ=255)
                                                                           tcacatcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgt            >  2:1809515/1‑93 (MQ=255)
                                                                            cacatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgtg           <  1:3089858/93‑1 (MQ=255)
                                                                              catcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGaa         <  2:1397355/93‑1 (MQ=255)
                                                                              catcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGaa         >  1:3728788/1‑93 (MQ=255)
                                                                               atcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAg        <  2:4522532/93‑1 (MQ=255)
                                                                                tcacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGt       >  2:4600050/1‑93 (MQ=255)
                                                                                 cacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTg      <  2:4112064/93‑1 (MQ=255)
                                                                                 cacaCTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTg      >  2:4254473/1‑93 (MQ=255)
                                                                                  acacTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGa     <  1:2913965/93‑1 (MQ=255)
                                                                                   cacTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGAt    <  2:1922427/93‑1 (MQ=255)
                                                                                    acTTTTACCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATa   <  1:1349133/93‑1 (MQ=255)
                                                                                            |                                                                                      
AAACGCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATAC  >  NC_000913/3071837‑3072014

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