Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 347,902 G→A intergenic (‑235/‑4) prpR ← / → prpB regulator for prp operon/putative phosphonomutase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913347,9020GA97.6% 563.2 / ‑2.2 164intergenic (‑235/‑4)prpR/prpBregulator for prp operon/putative phosphonomutase 2
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (84/76);  total (86/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 165 total aligned reads displayed.

AAATTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGGGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTAA  >  NC_000913/347812‑347991
                                                                                          |                                                                                         
aaaTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGc                                                                                         <  2:1546143/93‑1 (MQ=255)
 aaTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCt                                                                                        <  2:3529311/93‑1 (MQ=255)
   ttAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTAt                                                                                      >  1:2171503/1‑93 (MQ=255)
   ttAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTAt                                                                                      >  2:3831918/1‑93 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAACTAACAAGTTGATAACAAAGGATGAGCTATg                                                                                     >  1:4489682/1‑93 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATg                                                                                     <  2:501690/93‑1 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATg                                                                                     >  2:4491817/1‑93 (MQ=255)
     aaaCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGt                                                                                    >  1:533447/1‑93 (MQ=255)
      aaCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtc                                                                                   >  2:4534626/1‑93 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  >  2:924541/1‑93 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  <  2:3034576/93‑1 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  <  2:4349869/93‑1 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  >  2:580342/1‑93 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  >  2:736051/1‑93 (MQ=255)
           ttAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTac                                                                              >  1:1295350/1‑93 (MQ=255)
           ttAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTac                                                                              >  2:87222/1‑93 (MQ=255)
              aTTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACAct                                                                           >  1:358987/1‑93 (MQ=255)
                tttATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACActct                                                                         <  2:584260/93‑1 (MQ=255)
                 ttattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACtctc                                                                        <  1:770873/93‑1 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       >  2:306978/1‑93 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       <  2:302452/93‑1 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGCTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       >  2:2279854/1‑93 (MQ=255)
                   attaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCa                                                                      >  2:4090183/1‑93 (MQ=255)
                   attaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCa                                                                      <  2:1512999/93‑1 (MQ=255)
                   attaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCa                                                                      >  1:1032408/1‑93 (MQ=255)
                    ttaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAg                                                                     <  1:1644331/93‑1 (MQ=255)
                    ttaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAg                                                                     >  2:3702528/1‑93 (MQ=255)
                     taAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAgg                                                                    <  2:4502671/93‑1 (MQ=255)
                      aaGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGt                                                                   <  1:1764754/93‑1 (MQ=255)
                      aaGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGt                                                                   >  1:1930793/1‑93 (MQ=255)
                       aGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTa                                                                  >  2:1254963/1‑93 (MQ=255)
                        ggCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaa                                                                 <  1:2353253/93‑1 (MQ=255)
                         gCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaaa                                                                >  2:3094280/1‑93 (MQ=255)
                         gCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaaa                                                                <  1:728880/93‑1 (MQ=255)
                          cAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAg                                                               <  1:3031125/93‑1 (MQ=255)
                           aaTTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGc                                                              >  2:3581855/1‑93 (MQ=255)
                             ttGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGt                                                            >  2:285518/1‑93 (MQ=255)
                              tgtgGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGtt                                                           <  1:2169928/93‑1 (MQ=255)
                                tgGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACGCTCTCCAGGTAAAGCGTTTc                                                         <  1:3314739/93‑1 (MQ=255)
                                tgGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTc                                                         <  1:3974499/93‑1 (MQ=255)
                                 ggCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTcg                                                        >  2:2047953/1‑93 (MQ=255)
                                  gCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTGTCCAGGTAAAGCGTTTcgc                                                       <  1:4598074/93‑1 (MQ=255)
                                  gCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTcgc                                                       <  2:3448161/93‑1 (MQ=255)
                                  gCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTcgc                                                       <  1:1940497/93‑1 (MQ=255)
                                    acacCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCgcgc                                                     <  2:2698977/93‑1 (MQ=255)
                                     cacCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCt                                                    <  2:2384269/93‑1 (MQ=255)
                                      acCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTg                                                   <  2:3935773/93‑1 (MQ=255)
                                      acCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTg                                                   >  1:3754893/1‑93 (MQ=255)
                                         ccTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGccc                                                >  2:2830804/1‑91 (MQ=255)
                                         ccTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCAc                                                >  2:1851699/1‑93 (MQ=255)
                                          cTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACt                                               >  1:4676702/1‑93 (MQ=255)
                                           ttGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTg                                              <  1:4419415/93‑1 (MQ=255)
                                            tGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGa                                             <  2:4683605/93‑1 (MQ=255)
                                            tGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGa                                             <  2:1293087/93‑1 (MQ=255)
                                              cTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACt                                           <  1:1862116/93‑1 (MQ=255)
                                              cTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACt                                           >  2:332036/1‑93 (MQ=255)
                                                ttGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTaa                                         <  1:1465697/93‑1 (MQ=255)
                                                  gTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAg                                       <  1:4052985/93‑1 (MQ=255)
                                                  gTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAg                                       <  2:468116/93‑1 (MQ=255)
                                                  gTCTCTATCAACGCAAATAACAAGTGGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAg                                       <  1:4052984/93‑1 (MQ=255)
                                                   tCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGa                                      >  1:1179894/1‑93 (MQ=255)
                                                   tCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGa                                      >  1:2393089/1‑93 (MQ=255)
                                                    cTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaa                                     <  1:578137/93‑1 (MQ=255)
                                                     tttATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaaa                                    <  1:1785434/93‑1 (MQ=255)
                                                     tttATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaaa                                    >  1:3026338/1‑93 (MQ=255)
                                                       tATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAAt                                  >  2:1824852/1‑93 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACATTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATca                                >  2:3822690/1‑92 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATcc                                >  2:3822693/1‑93 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATcc                                <  2:3783375/93‑1 (MQ=255)
                                                          cAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCa                               <  2:3355405/93‑1 (MQ=255)
                                                           aaCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAt                              <  2:3169357/93‑1 (MQ=255)
                                                            aCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAtt                             >  2:2645528/1‑93 (MQ=255)
                                                            aCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAtt                             <  1:1513504/93‑1 (MQ=255)
                                                            aCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAtt                             <  1:3908889/93‑1 (MQ=255)
                                                             cGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTg                            >  1:3715587/1‑93 (MQ=255)
                                                               cAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCa                          <  2:3983130/93‑1 (MQ=255)
                                                                aaaTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGccg                         >  2:717042/1‑91 (MQ=255)
                                                                aaaTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAg                         <  1:913022/93‑1 (MQ=255)
                                                                  aTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAt                       <  1:1228712/93‑1 (MQ=255)
                                                                  aTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAt                       >  1:100000/1‑93 (MQ=255)
                                                                   tAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCATAtt                      >  2:1167612/1‑93 (MQ=255)
                                                                   tAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAtt                      <  2:4301925/93‑1 (MQ=255)
                                                                   tAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAtt                      <  2:1988535/93‑1 (MQ=255)
                                                                      cAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgtt                   <  1:449000/93‑1 (MQ=255)
                                                                      cAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgtt                   <  2:1549539/93‑1 (MQ=255)
                                                                         gTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGc                <  1:2428727/93‑1 (MQ=255)
                                                                          ttGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCa               >  1:1997149/1‑93 (MQ=255)
                                                                          ttGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCa               <  1:4004450/93‑1 (MQ=255)
                                                                           tGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCAc              >  1:2716729/1‑93 (MQ=255)
                                                                              tAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCAt           <  1:3117363/93‑1 (MQ=255)
                                                                               aaCAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATc          <  2:3852678/93‑1 (MQ=255)
                                                                                aCAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCa         >  2:3874825/1‑93 (MQ=255)
                                                                                 cAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCaa        >  2:2242765/1‑93 (MQ=255)
                                                                                  aaaGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAAc       <  1:2820708/93‑1 (MQ=255)
                                                                                  aaaGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAAc       >  1:3322894/1‑93 (MQ=255)
                                                                                  aaaGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAAc       >  2:4307169/1‑93 (MQ=255)
                                                                                      gATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTa   <  1:368179/93‑1 (MQ=255)
                                                                                      gATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTa   >  1:4699160/1‑93 (MQ=255)
                                                                                       aTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTaa  <  2:3501836/93‑1 (MQ=255)
                                                                                       aTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTaa  >  1:2557393/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAATTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGGGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTAA  >  NC_000913/347812‑347991

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