Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,428,098 G→A G126G (GGC→GGT aroE ← dehydroshikimate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,428,0980GA99.4% 1202.4 / ‑4.4 341G126G (GGC→GGTaroEdehydroshikimate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (143/195);  total (144/197)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 347 total aligned reads displayed.

TTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCGCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCCAG  >  NC_000913/3428008‑3428187
                                                                                          |                                                                                         
ttCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCAcc                                                                                         <  2:1886838/93‑1 (MQ=255)
 tCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCg                                                                                        <  2:449556/93‑1 (MQ=255)
 tCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCg                                                                                        >  1:7775422/1‑93 (MQ=255)
  ccGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGa                                                                                       <  1:1799598/93‑1 (MQ=255)
     cgcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATaa                                                                                    <  1:1515298/93‑1 (MQ=255)
      gcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAg                                                                                   <  2:3517674/93‑1 (MQ=255)
       cgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGc                                                                                  >  2:7819022/1‑93 (MQ=255)
       cgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGc                                                                                  <  2:4229712/93‑1 (MQ=255)
        gggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCa                                                                                 <  1:977791/93‑1 (MQ=255)
         ggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAg                                                                                >  2:4848156/1‑93 (MQ=255)
           aTACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGaa                                                                              <  1:5527959/93‑1 (MQ=255)
           aTACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGaa                                                                              <  2:5680899/93‑1 (MQ=255)
             aCCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATa                                                                            <  2:5256964/93‑1 (MQ=255)
              ccGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATAc                                                                           >  1:5903490/1‑93 (MQ=255)
              ccGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATAc                                                                           >  2:5743632/1‑93 (MQ=255)
              ccGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATAc                                                                           >  2:4081364/1‑93 (MQ=255)
                gTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGt                                                                         <  1:2398310/93‑1 (MQ=255)
                gTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGt                                                                         >  2:2688099/1‑93 (MQ=255)
                 tCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTa                                                                        >  2:5351250/1‑93 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                     >  2:1447254/1‑93 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                     >  1:6912284/1‑93 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                     >  1:7281863/1‑93 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                     <  2:2951817/93‑1 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                     <  1:447329/93‑1 (MQ=255)
                     aTTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAcc                                                                    >  1:2029777/1‑93 (MQ=255)
                      ttAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCCCCAGCACCGATAAGCAGAATACGTAAACCa                                                                   >  2:6578441/1‑93 (MQ=255)
                      ttAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCa                                                                   <  1:1208040/93‑1 (MQ=255)
                      ttAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCa                                                                   >  1:7012316/1‑93 (MQ=255)
                        aGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAgg                                                                 <  1:2589002/93‑1 (MQ=255)
                          ttATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGc                                                               <  1:7836974/93‑1 (MQ=255)
                          ttATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGc                                                               <  2:5062753/93‑1 (MQ=255)
                          ttATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGc                                                               >  2:3458421/1‑93 (MQ=255)
                              tGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGAt                                                           >  2:6861654/1‑93 (MQ=255)
                               gTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATa                                                          >  2:6154587/1‑93 (MQ=255)
                                tCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaa                                                         <  2:6880859/93‑1 (MQ=255)
                                tCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaa                                                         >  1:5287715/1‑93 (MQ=255)
                                tCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaa                                                         <  2:5460681/93‑1 (MQ=255)
                                 cACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaa                                                        <  2:6983195/93‑1 (MQ=255)
                                  aCCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaaa                                                       <  1:5096541/93‑1 (MQ=255)
                                  aCCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaaa                                                       >  1:4112503/1‑93 (MQ=255)
                                   ccGCACACTCCAGGGAAAGGAGTGGCAGAAGTACGCCGCGAGCCGCTCCACCAGCACCGATAAGCAGAATACAGCAACTAGGGCGGACaaaac                                                      >  2:2508366/1‑92 (MQ=255)
                                    cGCACAGTCCAGGGAAAGGATTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAaga                                                     <  1:5269358/93‑1 (MQ=255)
                                     gCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagac                                                    <  1:544510/93‑1 (MQ=255)
                                     gCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagac                                                    <  1:3447677/93‑1 (MQ=255)
                                      cacaGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagaca                                                   <  2:7381676/93‑1 (MQ=255)
                                       acaGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacag                                                  >  1:7417525/1‑93 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                 <  1:126359/93‑1 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                 <  2:5616829/93‑1 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                 <  2:2283393/93‑1 (MQ=255)
                                         aGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCGCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacagac                                                <  2:2820209/93‑1 (MQ=255)
                                             cAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTc                                            >  2:6470863/1‑93 (MQ=255)
                                              aGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTcc                                           <  2:5892026/93‑1 (MQ=255)
                                              aGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTcc                                           <  1:1456135/93‑1 (MQ=255)
                                               gggAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCa                                          >  1:3618459/1‑93 (MQ=255)
                                               gggAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCa                                          <  1:5801865/93‑1 (MQ=255)
                                                ggAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAg                                         >  1:2246253/1‑93 (MQ=255)
                                                 gAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGa                                        <  1:384831/93‑1 (MQ=255)
                                                  aaaGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGAt                                       <  1:3738191/93‑1 (MQ=255)
                                                   aaGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATc                                      <  2:2011243/93‑1 (MQ=255)
                                                   aaGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATc                                      >  2:2331828/1‑93 (MQ=255)
                                                    aGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCg                                     >  2:4359279/1‑93 (MQ=255)
                                                     ggAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGc                                    >  2:2185062/1‑93 (MQ=255)
                                                       aGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCtt                                  <  1:1007697/93‑1 (MQ=255)
                                                       aGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCtt                                  <  1:881159/93‑1 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTa                                 >  2:892252/1‑93 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTa                                 >  2:2617639/1‑93 (MQ=255)
                                                          ggCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAAc                               <  2:7618682/93‑1 (MQ=255)
                                                          ggCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAAc                               <  2:3690034/93‑1 (MQ=255)
                                                          ggCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAAc                               <  1:802148/93‑1 (MQ=255)
                                                            cAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACaa                             <  1:1679364/93‑1 (MQ=255)
                                                             agtagtACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAg                            <  2:2229221/93‑1 (MQ=255)
                                                             agtagtACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAg                            >  2:4960936/1‑93 (MQ=255)
                                                              gtagtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGc                           >  1:2953093/1‑93 (MQ=255)
                                                               tagtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGcc                          <  2:6811774/93‑1 (MQ=255)
                                                               tagtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGcc                          >  1:3501917/1‑93 (MQ=255)
                                                               tagtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGcc                          <  1:6675108/93‑1 (MQ=255)
                                                               tagtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGcc                          <  1:5665752/93‑1 (MQ=255)
                                                                agtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCt                         >  2:1706828/1‑93 (MQ=255)
                                                                agtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCt                         <  2:5126897/93‑1 (MQ=255)
                                                                     gccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACAcc                    >  2:5489557/1‑93 (MQ=255)
                                                                     gccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACAcc                    <  1:7094781/93‑1 (MQ=255)
                                                                     cccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACAcc                    <  1:3685491/92‑1 (MQ=255)
                                                                      ccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCa                   <  1:4087855/93‑1 (MQ=255)
                                                                       cgcgAGATGCTCCACCAGCACCTATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCAt                  >  1:1387021/1‑93 (MQ=255)
                                                                       cgcgAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCAt                  <  1:1136489/93‑1 (MQ=255)
                                                                        gcgAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGAATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATc                 <  2:3250922/93‑1 (MQ=255)
                                                                          gagaTGCTCCACCAGCCCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCgg               >  2:3850895/1‑93 (MQ=255)
                                                                           agaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGt              >  1:3296012/1‑93 (MQ=255)
                                                                             aTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTAt            <  2:7479405/93‑1 (MQ=255)
                                                                              tGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTAtt           >  1:3066401/1‑93 (MQ=255)
                                                                                 tCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc        >  1:6470530/1‑93 (MQ=255)
                                                                                 tCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc        <  1:4832291/93‑1 (MQ=255)
                                                                                 tCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc        >  1:6373840/1‑93 (MQ=255)
                                                                                 cccaccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc        <  1:731578/92‑1 (MQ=255)
                                                                                 cccaccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc        <  1:4376346/92‑1 (MQ=255)
                                                                                   caccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCAc      <  2:1214836/93‑1 (MQ=255)
                                                                                   caccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCAc      >  1:5211061/1‑93 (MQ=255)
                                                                                    accaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCAcc     >  1:4732072/1‑93 (MQ=255)
                                                                                     ccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCAccc    >  2:362135/1‑93 (MQ=255)
                                                                                       aGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCcag  >  2:728036/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCGCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCCAG  >  NC_000913/3428008‑3428187

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