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,9020GA99.6% 956.2 / ‑5.7 258intergenic (‑235/‑4)prpR/prpBregulator for prp operon/putative phosphonomutase 2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (146/111);  total (147/111)
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.78e-01

Only 100 of 258 total aligned reads displayed.

AAATTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGGGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTAA  >  NC_000913/347812‑347991
                                                                                          |                                                                                         
aaaTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGc                                                                                         >  2:2067020/1‑93 (MQ=255)
 aaTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCt                                                                                        <  1:4338338/93‑1 (MQ=255)
 aaTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCt                                                                                        <  2:233704/93‑1 (MQ=255)
  aTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTa                                                                                       <  2:3735867/93‑1 (MQ=255)
  aTTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTa                                                                                       >  1:635321/1‑93 (MQ=255)
   ttAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTAt                                                                                      >  2:244988/1‑93 (MQ=255)
   ttAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTAt                                                                                      >  2:6534652/1‑93 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGGGCTATg                                                                                     >  1:630633/1‑93 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATg                                                                                     <  2:4950395/93‑1 (MQ=255)
    tAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATg                                                                                     >  2:512896/1‑93 (MQ=255)
     aaaCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGt                                                                                    >  1:3921653/1‑93 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  <  1:1108084/93‑1 (MQ=255)
       aCATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtct                                                                                  >  1:4285700/1‑93 (MQ=255)
        cATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtctc                                                                                 <  1:2730696/93‑1 (MQ=255)
        cATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGtctc                                                                                 >  1:6663293/1‑93 (MQ=255)
          tttAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTa                                                                               >  2:6262036/1‑93 (MQ=255)
          tttAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTa                                                                               >  1:4380696/1‑93 (MQ=255)
             aaTTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTacac                                                                            >  2:3286901/1‑93 (MQ=255)
                 ttattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACtctc                                                                        <  2:6596172/93‑1 (MQ=255)
                  tattaAGGCAATTTTGGCACACCCATTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       <  1:1061619/93‑1 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       >  2:2165425/1‑93 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       >  2:5964135/1‑93 (MQ=255)
                  tattaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTcc                                                                       >  1:1729304/1‑93 (MQ=255)
                   attaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCa                                                                      <  2:4000567/93‑1 (MQ=255)
                   attaAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCa                                                                      <  2:678110/93‑1 (MQ=255)
                     taAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAgg                                                                    <  1:4043922/93‑1 (MQ=255)
                      aaGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGt                                                                   <  1:1963163/93‑1 (MQ=255)
                      aaGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGt                                                                   <  2:3479208/93‑1 (MQ=255)
                       aGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTa                                                                  >  2:767506/1‑93 (MQ=255)
                       aGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTa                                                                  >  1:3711919/1‑93 (MQ=255)
                        ggCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaa                                                                 >  1:5159396/1‑93 (MQ=255)
                         gCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaaa                                                                >  2:1627435/1‑93 (MQ=255)
                         gCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTaaa                                                                >  1:2217353/1‑93 (MQ=255)
                          cAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAg                                                               <  2:4721170/93‑1 (MQ=255)
                          cAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAg                                                               >  2:3545443/1‑93 (MQ=255)
                          cAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAg                                                               <  1:4146100/93‑1 (MQ=255)
                            aTTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCg                                                             >  2:1016101/1‑93 (MQ=255)
                            aTTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCg                                                             >  1:885602/1‑93 (MQ=255)
                              tgtgGCTCTCCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGtt                                                           <  1:4804644/93‑1 (MQ=255)
                              tgtgGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGtt                                                           <  1:3374977/93‑1 (MQ=255)
                                tgGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTc                                                         <  1:2940850/93‑1 (MQ=255)
                                   cacaCCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTcgcg                                                      <  1:1651774/93‑1 (MQ=255)
                                     cacCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGCTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCt                                                    <  2:4703510/93‑1 (MQ=255)
                                     cacCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTCAAGCGTATCGCGCt                                                    >  1:719414/1‑93 (MQ=255)
                                         ccTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCAc                                                <  2:4227958/93‑1 (MQ=255)
                                         ccTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCAc                                                >  1:3463551/1‑93 (MQ=255)
                                          cTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACt                                               >  2:6416785/1‑93 (MQ=255)
                                              cTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACt                                           >  2:629794/1‑93 (MQ=255)
                                               tttGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTa                                          >  1:1939993/1‑93 (MQ=255)
                                                ttGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTaa                                         <  2:4256738/93‑1 (MQ=255)
                                                ttGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTaa                                         >  2:1152109/1‑93 (MQ=255)
                                                ttGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTaa                                         >  2:4373395/1‑93 (MQ=255)
                                                  gTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAg                                       >  1:5844594/1‑93 (MQ=255)
                                                  gTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAg                                       >  1:3672040/1‑93 (MQ=255)
                                                   tCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGa                                      >  1:540320/1‑93 (MQ=255)
                                                    cTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaa                                     >  1:2197860/1‑93 (MQ=255)
                                                    cTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaa                                     >  1:5739348/1‑93 (MQ=255)
                                                     tttATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGaaa                                    >  2:6018576/1‑93 (MQ=255)
                                                       tATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAAt                                  >  2:5879532/1‑93 (MQ=255)
                                                       tATCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAAt                                  <  1:4965336/93‑1 (MQ=255)
                                                        aTCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATc                                 <  1:1184882/93‑1 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATcc                                >  1:1814464/1‑93 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATcc                                >  1:3206497/1‑93 (MQ=255)
                                                         tCAACGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATcc                                >  1:2619549/1‑93 (MQ=255)
                                                           aaCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAt                              >  2:2735074/1‑93 (MQ=255)
                                                            aCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAtt                             <  2:1334512/93‑1 (MQ=255)
                                                            aCGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCAtt                             <  2:5250235/93‑1 (MQ=255)
                                                             cGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTg                            <  2:3806608/93‑1 (MQ=255)
                                                             cGCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTg                            <  1:5627555/93‑1 (MQ=255)
                                                              gCAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGc                           >  2:460385/1‑93 (MQ=255)
                                                               cAAATAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCa                          <  1:4860712/93‑1 (MQ=255)
                                                                 aaTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGa                        <  1:2902329/93‑1 (MQ=255)
                                                                 aaTAACAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGa                        <  1:5985744/93‑1 (MQ=255)
                                                                    aaCAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttg                     >  2:5364613/1‑93 (MQ=255)
                                                                     aCAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgt                    >  2:301894/1‑93 (MQ=255)
                                                                     aCAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgt                    >  2:6374678/1‑93 (MQ=255)
                                                                     aCAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgt                    <  2:505399/93‑1 (MQ=255)
                                                                      cAAGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgtt                   >  2:5660606/1‑93 (MQ=255)
                                                                       aaGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgttg                  <  2:6171127/93‑1 (MQ=255)
                                                                       aaGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgttg                  >  2:855164/1‑93 (MQ=255)
                                                                       aaGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgttg                  >  2:1913400/1‑93 (MQ=255)
                                                                       aaGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGAttgttg                  <  1:5304140/93‑1 (MQ=255)
                                                                        aGTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTgg                 >  1:11984/1‑93 (MQ=255)
                                                                         gTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGc                <  2:4055207/93‑1 (MQ=255)
                                                                         gTTGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGc                >  1:3235337/1‑93 (MQ=255)
                                                                           tGATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCAc              >  1:3785580/1‑93 (MQ=255)
                                                                            gATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCAcc             >  2:4493379/1‑93 (MQ=255)
                                                                            gATAACAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCAcc             <  2:2087063/93‑1 (MQ=255)
                                                                                 cAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCaa        >  2:3179758/1‑93 (MQ=255)
                                                                                 cAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCaa        <  1:2415055/93‑1 (MQ=255)
                                                                                 cAAAGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCaa        >  2:4857884/1‑93 (MQ=255)
                                                                                  aaaGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAAc       <  2:6037412/93‑1 (MQ=255)
                                                                                    aGGATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGc     >  1:4440002/1‑93 (MQ=255)
                                                                                     ggATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCt    >  2:4171502/1‑93 (MQ=255)
                                                                                     ggATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCt    >  2:559522/1‑93 (MQ=255)
                                                                                     ggATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCt    >  2:2689975/1‑93 (MQ=255)
                                                                                      gATGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTa   >  1:6443856/1‑93 (MQ=255)
                                                                                       aTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTaa  <  1:2996085/93‑1 (MQ=255)
                                                                                       aTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTaa  <  1:3317777/93‑1 (MQ=255)
                                                                                       aTGAGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTaa  <  1:858977/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAATTAAACATTTAATTTTATTAAGGCAATTGTGGCACACCCCTTGCTTTGTCTTTATCAACGCAAATAACAAGTTGATAACAAAGGATGGGCTATGTCTCTACACTCTCCAGGTAAAGCGTTTCGCGCTGCACTGACTAAAGAAAATCCATTGCAGATTGTTGGCACCATCAACGCTAA  >  NC_000913/347812‑347991

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: