Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 367,441 A→G D106D (GAT→GAC mhpR ← transcriptional regulator for mhp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913367,4410AG98.0% 481.3 / ‑0.7 152D106D (GAT→GACmhpRtranscriptional regulator for mhp operon
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (79/69);  total (81/71)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 164 total aligned reads displayed.

GCCCAGCAGTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTATCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAGTT  >  NC_000913/367351‑367529
                                                                                          |                                                                                        
 cccAGCAGTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtca                                                                                       >  1:4648790/1‑93 (MQ=255)
  ccAGCAGTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtcat                                                                                      <  2:4007283/93‑1 (MQ=255)
   cagcagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtcatc                                                                                     >  2:1358621/1‑93 (MQ=255)
   cagcagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtcatc                                                                                     <  1:341707/93‑1 (MQ=255)
   cagcagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtcatc                                                                                     <  1:3061578/93‑1 (MQ=255)
   cagcagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTgtcatc                                                                                     <  2:3499648/93‑1 (MQ=255)
     gcagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCgg                                                                                   <  2:1507106/93‑1 (MQ=255)
      cagTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGa                                                                                  <  2:629222/93‑1 (MQ=255)
       agTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAg                                                                                 >  1:925059/1‑93 (MQ=255)
        gTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAgg                                                                                <  2:2213508/93‑1 (MQ=255)
          ggggCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAgggg                                                                              >  1:2761073/1‑93 (MQ=255)
          ggggCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAgggg                                                                              >  2:3583004/1‑93 (MQ=255)
           gggCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGc                                                                             >  2:3663631/1‑93 (MQ=255)
           gggCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGc                                                                             <  2:564343/93‑1 (MQ=255)
             gccgccAGTGGAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTa                                                                           <  1:2016049/93‑1 (MQ=255)
             gccgccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCCCTTTGATGGTCAGTCGAAAACTGTAATCGGAGGGGCTa                                                                           <  2:1591891/93‑1 (MQ=255)
              ccgccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAc                                                                          <  1:109991/93‑1 (MQ=255)
               cgccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcg                                                                         <  2:2267299/93‑1 (MQ=255)
               cgccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATAGGAGGGGCTAcg                                                                         <  1:198029/93‑1 (MQ=255)
               cgccAGTCCAGAAACCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcg                                                                         <  2:1731521/93‑1 (MQ=255)
                gccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcgg                                                                        <  1:2114421/93‑1 (MQ=255)
                gccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcgg                                                                        >  2:2202736/1‑93 (MQ=255)
                 ccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcggc                                                                       >  2:434260/1‑93 (MQ=255)
                 ccAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcggc                                                                       >  1:3153931/1‑93 (MQ=255)
                  cAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcggcg                                                                      >  1:1794454/1‑93 (MQ=255)
                   aGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTAcggcgg                                                                     >  1:2199035/1‑93 (MQ=255)
                     tGCAGAAATCCACTGTTCGTCACGAAGTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGAc                                                                   <  1:4330486/93‑1 (MQ=255)
                        aGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACata                                                                <  2:2735743/93‑1 (MQ=255)
                        aGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACata                                                                <  2:1400010/93‑1 (MQ=255)
                         gAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACatat                                                               >  2:27762/1‑93 (MQ=255)
                         gAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACatat                                                               >  1:750322/1‑93 (MQ=255)
                         gAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACatat                                                               <  1:506377/93‑1 (MQ=255)
                         gAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACatac                                                               >  1:4165316/1‑92 (MQ=255)
                            aTCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATccc                                                            <  1:1167958/93‑1 (MQ=255)
                             tCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCct                                                           >  2:4005386/1‑93 (MQ=255)
                               cACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCctct                                                         <  2:741313/93‑1 (MQ=255)
                               cACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCctct                                                         <  1:1032273/93‑1 (MQ=255)
                                aCTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCtt                                                        >  2:4152629/1‑93 (MQ=255)
                                  tGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTcc                                                      >  1:660754/1‑93 (MQ=255)
                                   gTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCAGAGGGGCTACGGCGGACATATCCCTCTTCCt                                                     >  2:1079265/1‑93 (MQ=255)
                                      tgTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCa                                                  <  2:2567315/92‑1 (MQ=255)
                                      cGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCa                                                  <  1:3714891/93‑1 (MQ=255)
                                      cGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCa                                                  >  1:2713825/1‑93 (MQ=255)
                                        tCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGc                                                >  1:1652874/1‑93 (MQ=255)
                                         cACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCg                                               >  2:4519228/1‑93 (MQ=255)
                                           cGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtc                                             <  1:4556342/93‑1 (MQ=255)
                                           cGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtc                                             <  1:3008513/93‑1 (MQ=255)
                                           cGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtc                                             >  2:1457504/1‑93 (MQ=255)
                                           cGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATAGCCCTCTTCCTGCAGCGtc                                             <  1:4556341/93‑1 (MQ=255)
                                            gAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtct                                            <  2:1987072/93‑1 (MQ=255)
                                             aaaTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtctc                                           >  1:716969/1‑93 (MQ=255)
                                             aaaTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGtctc                                           <  2:2332313/93‑1 (MQ=255)
                                              aaTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTcc                                          <  2:3749824/93‑1 (MQ=255)
                                              aaTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGAGTCTcc                                          >  2:3270915/1‑93 (MQ=255)
                                               aTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAGAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTTTCca                                         <  1:1380518/93‑1 (MQ=255)
                                               aTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCca                                         <  1:3865781/93‑1 (MQ=255)
                                               aTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCca                                         >  1:2395592/1‑93 (MQ=255)
                                               aTCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCca                                         >  2:3302446/1‑93 (MQ=255)
                                                 ccTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCcagc                                       <  2:218093/93‑1 (MQ=255)
                                                  cTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCcagca                                      >  1:970121/1‑93 (MQ=255)
                                                  cTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCcagca                                      >  1:4227257/1‑93 (MQ=255)
                                                   ttCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCcagcag                                     <  1:1236795/93‑1 (MQ=255)
                                                     cGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTc                                   <  2:4428210/93‑1 (MQ=255)
                                                     cGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTc                                   <  1:2722595/93‑1 (MQ=255)
                                                      gCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGCGGGGCTACGGCGGACATATCCCTCTTCCTGCAGAGTCTCCAGGAGTcg                                  >  1:385404/1‑93 (MQ=255)
                                                       cTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgc                                 >  2:2017900/1‑93 (MQ=255)
                                                       cTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgc                                 >  2:1867171/1‑93 (MQ=255)
                                                       cTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgc                                 <  1:3696490/93‑1 (MQ=255)
                                                        ttAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGATGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgcc                                >  1:491594/1‑93 (MQ=255)
                                                        ttAATTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgcc                                <  2:2998349/93‑1 (MQ=255)
                                                        ttAATTGCCGCACTTTGATGGTCAGTCGAAAACTATCATCGGAGGGGCTACGGCGGACATATTCCTATTCCTGCAGCGGCTCCAGCACTcgcc                                >  2:4372998/1‑93 (MQ=255)
                                                          aaTTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgccgc                              <  2:2843871/93‑1 (MQ=255)
                                                           aTTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTcgccgcc                             >  1:890457/1‑92 (MQ=255)
                                                           aTTGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTGGCCGca                             >  2:2097210/1‑93 (MQ=255)
                                                             tGCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGcaca                           <  2:3395072/93‑1 (MQ=255)
                                                              gCCGCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAg                          >  2:3892126/1‑93 (MQ=255)
                                                                 gCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAgtgg                       >  1:914550/1‑93 (MQ=255)
                                                                 gCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAgtgg                       >  2:1800134/1‑93 (MQ=255)
                                                                 gCACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGAGTCTCCAGCCGTCGCCGCACAgtgg                       >  1:3653539/1‑93 (MQ=255)
                                                                  cACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAgtggt                      <  1:4664648/93‑1 (MQ=255)
                                                                  cACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAgtggg                      >  1:2618529/1‑92 (MQ=255)
                                                                  cACTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCATGCAGCGTCTCCAGCAGTCGCCGCACAggggg                      >  2:1785477/1‑92 (MQ=255)
                                                                    cTTTGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGc                    >  1:3677385/1‑93 (MQ=255)
                                                                     tttGATGGTCAGTCGAAAGCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCTGCACAGTGGTGCg                   <  1:2451889/93‑1 (MQ=255)
                                                                       tGATGGTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGAt                 <  2:3252173/93‑1 (MQ=255)
                                                                           ggTCAGTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAg             <  2:61264/93‑1 (MQ=255)
                                                                                gTCGAAAACTGTCATCGGAGGGGCTAGGGCGGACATATCCCTCTTCCTGCAGGGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGc        >  1:2231636/1‑93 (MQ=255)
                                                                                gTCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGc        >  1:3569983/1‑93 (MQ=255)
                                                                                agcGAAAACCGTCATCGCAGGGGCTACTGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGc        <  1:1691878/91‑1 (MQ=255)
                                                                                 tCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGcg       >  2:3253667/1‑92 (MQ=255)
                                                                                 tCGAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCt       >  2:4143038/1‑93 (MQ=255)
                                                                                   gAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTgg     >  2:1964633/1‑92 (MQ=255)
                                                                                   gAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGa     >  2:3864645/1‑93 (MQ=255)
                                                                                   gAAAACTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGa     >  2:1091346/1‑93 (MQ=255)
                                                                                    aaaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAg    >  1:3308620/1‑93 (MQ=255)
                                                                                    aaaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAg    >  1:3308621/1‑93 (MQ=255)
                                                                                    aaaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAg    <  2:4487641/93‑1 (MQ=255)
                                                                                    aaaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAg    >  2:4244800/1‑93 (MQ=255)
                                                                                     aaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAGt   >  1:3920838/1‑93 (MQ=255)
                                                                                     aaaCTGTCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAGt   <  1:510150/93‑1 (MQ=255)
                                                                                          |                                                                                        
GCCCAGCAGTGGGGCCGCCAGTGCAGAAATCCACTGTTCGTCACGAAATCCTTCGCTTAATTGCCGCACTTTGATGGTCAGTCGAAAACTATCATCGGAGGGGCTACGGCGGACATATCCCTCTTCCTGCAGCGTCTCCAGCAGTCGCCGCACAGTGGTGCGATGCAGGCCGCTGAGTT  >  NC_000913/367351‑367529

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