Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,458,282 A→G G651G (GGT→GGC)  nrdD ← anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,458,2820AG99.6% 802.7 / ‑4.8 251G651G (GGT→GGC) nrdDanaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (114/136);  total (115/136)
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.68e-01

Only 100 of 262 total aligned reads displayed.

CCGCACACGCGGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCTTC  >  NC_000913/4458193‑4458371
                                                                                         |                                                                                         
 cGCACACGCGGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCgg                                                                                       <  1:4360705/93‑1 (MQ=255)
  gCACACGCGGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGt                                                                                      <  1:3943515/93‑1 (MQ=255)
   cacaCGCGGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTa                                                                                     >  1:4912287/1‑93 (MQ=255)
       cgcgGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGc                                                                                 <  2:344745/93‑1 (MQ=255)
       cgcgGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGc                                                                                 <  2:3046083/93‑1 (MQ=255)
        gcggcgAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCAGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGcc                                                                                >  1:1256002/1‑93 (MQ=255)
         cggcgAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCca                                                                               >  2:3653122/1‑93 (MQ=255)
           gcgAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCTGGTAAAGCcaca                                                                             <  1:5926077/93‑1 (MQ=255)
            cgAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCacac                                                                            <  1:5256526/93‑1 (MQ=255)
               gTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCg                                                                         <  1:4349698/93‑1 (MQ=255)
                 tACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTa                                                                       >  2:6227538/1‑93 (MQ=255)
                 tACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGATGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTa                                                                       <  1:1725229/93‑1 (MQ=255)
                  aCCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAg                                                                      <  2:5249015/93‑1 (MQ=255)
                  aCCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAg                                                                      <  2:4886178/93‑1 (MQ=255)
                  aCCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAg                                                                      <  2:2878110/93‑1 (MQ=255)
                   ccGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGc                                                                     <  2:2475276/93‑1 (MQ=255)
                    cGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCa                                                                    >  2:3185205/1‑93 (MQ=255)
                     gACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCAc                                                                   >  2:4786492/1‑93 (MQ=255)
                      acacacGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACt                                                                  >  2:3347001/1‑93 (MQ=255)
                       cacacGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACtc                                                                 >  2:4145960/1‑93 (MQ=255)
                        acacGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACtca                                                                <  1:989699/93‑1 (MQ=255)
                         cacGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACtcat                                                               >  1:5447595/1‑93 (MQ=255)
                         cacGGGAGGCGTCATGGTTACCACATTTCGGGCAAGCGAAGCCTTTGCTCGTGCACTCGAACTCGCCGGTAAAGCCACACTCCTAGCACtcat                                                               >  1:109148/1‑93 (MQ=255)
                          acGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACtcatc                                                              <  1:5750131/93‑1 (MQ=255)
                           cGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACtcatca                                                             <  1:1350563/93‑1 (MQ=255)
                            gggAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCaa                                                            <  2:64489/93‑1 (MQ=255)
                             ggAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAAt                                                           >  1:317280/1‑93 (MQ=255)
                                ggCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCgg                                                        <  1:1596157/93‑1 (MQ=255)
                                 gCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGgt                                                       <  2:3239711/93‑1 (MQ=255)
                                   gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGgtgt                                                     >  2:2217160/1‑93 (MQ=255)
                                   gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGgtgt                                                     >  1:1123407/1‑93 (MQ=255)
                                    tCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTa                                                    >  2:1608648/1‑93 (MQ=255)
                                    tCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTa                                                    >  1:5594515/1‑93 (MQ=255)
                                    tCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTa                                                    >  1:500885/1‑93 (MQ=255)
                                     cATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTAt                                                   <  1:890632/93‑1 (MQ=255)
                                      aTGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTAtt                                                  >  1:4133864/1‑93 (MQ=255)
                                       tGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTg                                                 <  1:6123696/93‑1 (MQ=255)
                                         gTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGt                                               <  1:821300/93‑1 (MQ=255)
                                         gTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGt                                               <  1:1203867/93‑1 (MQ=255)
                                           tACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGc                                             <  2:5074950/93‑1 (MQ=255)
                                           tACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGc                                             <  2:2107287/93‑1 (MQ=255)
                                            aCCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGcc                                            <  2:858939/93‑1 (MQ=255)
                                            aCCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGcc                                            <  2:1870022/93‑1 (MQ=255)
                                              cacaTTTTGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGt                                          >  1:3442308/1‑93 (MQ=255)
                                               ccaTTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTa                                         <  2:496611/92‑1 (MQ=255)
                                                caTTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTaa                                        <  1:2481680/93‑1 (MQ=255)
                                                 aTTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAAt                                       >  1:4604491/1‑93 (MQ=255)
                                                 aTTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAAt                                       <  1:2311150/93‑1 (MQ=255)
                                                  tttCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATa                                      >  1:5188177/1‑93 (MQ=255)
                                                  tttCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATa                                      <  1:5316687/93‑1 (MQ=255)
                                                    tCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACg                                    >  1:2181674/1‑93 (MQ=255)
                                                     cGGGCAAGGGAAGCCTTTGCTGGTGCCCTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACgg                                   <  2:1223291/93‑1 (MQ=255)
                                                      gggCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGt                                  >  1:712643/1‑93 (MQ=255)
                                                      gggCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGt                                  <  2:3388140/93‑1 (MQ=255)
                                                         cAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACa                               <  2:2665439/93‑1 (MQ=255)
                                                         cAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACa                               <  1:3047410/93‑1 (MQ=255)
                                                         cAAGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACa                               <  2:5348372/93‑1 (MQ=255)
                                                          aaGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACAt                              <  2:5857472/93‑1 (MQ=255)
                                                          aaGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACAt                              >  1:3016365/1‑93 (MQ=255)
                                                           aGTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATg                             >  2:1534025/1‑93 (MQ=255)
                                                            gTGAAGCCTTTGCTGGTGCCCTCGACCTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGc                            <  2:5676663/93‑1 (MQ=255)
                                                            gTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGc                            >  1:2546956/1‑93 (MQ=255)
                                                            gTGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACCCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGc                            <  1:4051521/93‑1 (MQ=255)
                                                             tGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCt                           <  1:3319661/93‑1 (MQ=255)
                                                             tGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCt                           <  2:1061566/93‑1 (MQ=255)
                                                             tGAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCt                           <  2:713662/93‑1 (MQ=255)
                                                              gAAGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTg                          <  1:3554446/93‑1 (MQ=255)
                                                                aGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGAt                        <  2:4211213/93‑1 (MQ=255)
                                                                aGCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGAt                        >  1:1950951/1‑93 (MQ=255)
                                                                 gCCTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATa                       <  1:4007585/93‑1 (MQ=255)
                                                                  ccTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAg                      >  2:3960840/1‑93 (MQ=255)
                                                                  ccTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAg                      >  1:2630943/1‑93 (MQ=255)
                                                                  ccTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAg                      <  2:2355946/93‑1 (MQ=255)
                                                                   cTTTGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGc                     >  1:3887086/1‑93 (MQ=255)
                                                                    tttGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCt                    >  1:4762660/1‑93 (MQ=255)
                                                                    tttGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCt                    <  2:1419939/93‑1 (MQ=255)
                                                                    tttGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCt                    <  2:5989564/93‑1 (MQ=255)
                                                                     ttGCTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTg                   <  1:2825770/93‑1 (MQ=255)
                                                                        cTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTaa                <  2:2424884/93‑1 (MQ=255)
                                                                        cTGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTaa                <  2:2125177/93‑1 (MQ=255)
                                                                         tGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAAt               >  1:2929035/1‑93 (MQ=255)
                                                                         tGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAAt               >  2:5494223/1‑93 (MQ=255)
                                                                         tGGTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTAGTGCCGTAATACGGTACATGCTGATAGCTGTAAt               >  1:5771381/1‑93 (MQ=255)
                                                                          ggTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATc              <  2:2547161/93‑1 (MQ=255)
                                                                          ggTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATc              >  2:1708435/1‑93 (MQ=255)
                                                                           gTGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATcc             >  1:4699982/1‑93 (MQ=255)
                                                                            tGCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATccc            <  1:3225582/93‑1 (MQ=255)
                                                                             gCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCa           <  2:370683/93‑1 (MQ=255)
                                                                             gCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCa           <  1:1628722/93‑1 (MQ=255)
                                                                             gCACTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCa           <  2:6076331/93‑1 (MQ=255)
                                                                                cTCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGAc        >  2:3325784/1‑93 (MQ=255)
                                                                                cTCGAACTCACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGAc        >  1:466251/1‑93 (MQ=255)
                                                                                 tCGAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACa       <  2:6098415/93‑1 (MQ=255)
                                                                                   gAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATc     <  1:1700204/93‑1 (MQ=255)
                                                                                   gAACTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATc     <  2:4363072/93‑1 (MQ=255)
                                                                                    aaCTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCt    <  2:3611461/93‑1 (MQ=255)
                                                                                     aCTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCtt   <  1:334693/93‑1 (MQ=255)
                                                                                     aCTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCtt   >  2:2920637/1‑93 (MQ=255)
                                                                                      cTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCTTc  >  1:5965054/1‑93 (MQ=255)
                                                                                      cTCGCCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCTTc  >  1:2692396/1‑93 (MQ=255)
                                                                                         |                                                                                         
CCGCACACGCGGCGAGTTACCGACACACGGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAGCTGTAATCCCAGACATCTTC  >  NC_000913/4458193‑4458371

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: ‑