Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 200,585 A→G Q35R (CAG→CGG)  hlpA → histone‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913200,5850AG98.8% 529.3 / ‑3.0 166Q35R (CAG→CGG) hlpAhistone‑like protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (66/98);  total (67/99)
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.12e-01

Only 100 of 180 total aligned reads displayed.

GTTATTAGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGC  >  NC_000913/200493‑200674
                                                                                            |                                                                                         
gTTATTAGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCACAATGGCCAGCCTGTTCcagca                                                                                           >  1:2472313/1‑93 (MQ=255)
  tattaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCggg                                                                                         >  1:302447/1‑93 (MQ=255)
  tattaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCggg                                                                                         >  1:425498/1‑93 (MQ=255)
   attaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCgggc                                                                                        >  2:4292093/1‑92 (MQ=255)
   attaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGt                                                                                        <  1:3074157/93‑1 (MQ=255)
   attaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGt                                                                                        >  1:964073/1‑93 (MQ=255)
    ttaGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTa                                                                                       <  2:1161182/93‑1 (MQ=255)
          gCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAg                                                                                 >  2:3305523/1‑93 (MQ=255)
            aGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGaa                                                                               >  1:1993811/1‑93 (MQ=255)
             ggTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGaaa                                                                              >  2:2218509/1‑93 (MQ=255)
             ggTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGaaa                                                                              >  2:4518754/1‑93 (MQ=255)
             ggTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGaaa                                                                              <  1:3763340/93‑1 (MQ=255)
               tctcGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAAc                                                                            <  2:432895/93‑1 (MQ=255)
                ctcGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGACAAcc                                                                           >  1:1003388/1‑93 (MQ=255)
                ctcGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAAcc                                                                           <  2:2180661/93‑1 (MQ=255)
                  cGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCgg                                                                         >  2:4438239/1‑93 (MQ=255)
                    gTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGgtg                                                                       >  1:1108637/1‑93 (MQ=255)
                      ttAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCATCGGGTAGCGCAGAAAACCGGTGtt                                                                     <  2:4640741/93‑1 (MQ=255)
                         gCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAATCGGTGTTTCt                                                                  >  1:4195398/1‑93 (MQ=255)
                         gCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCt                                                                  >  2:3824186/1‑93 (MQ=255)
                         gCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCt                                                                  <  1:1781018/93‑1 (MQ=255)
                         gCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCt                                                                  <  1:3281797/93‑1 (MQ=255)
                           aCTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTaa                                                                <  2:933193/93‑1 (MQ=255)
                            cTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAac                                                               >  2:2724909/1‑93 (MQ=255)
                            cTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAac                                                               >  1:4029771/1‑93 (MQ=255)
                             tGGCAACTTGTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACAGGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAaca                                                              <  2:3898188/93‑1 (MQ=255)
                             tGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAaca                                                              <  1:3762685/93‑1 (MQ=255)
                                cAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGc                                                           <  2:3662844/93‑1 (MQ=255)
                                cAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGc                                                           <  1:4531773/93‑1 (MQ=255)
                                cAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGAAGCGCAGAAACCCGGTGTTTCTAGCACGc                                                           >  1:1814831/1‑93 (MQ=255)
                                  aCTTCTGCTCAGTCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTg                                                         >  1:3108060/1‑93 (MQ=255)
                                  aCTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTg                                                         <  1:1748247/93‑1 (MQ=255)
                                  aCTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTg                                                         >  1:3624337/1‑93 (MQ=255)
                                  aCTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTg                                                         <  1:3429358/93‑1 (MQ=255)
                                   cTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTgg                                                        <  2:1370360/93‑1 (MQ=255)
                                    ttCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGa                                                       >  1:2814117/1‑93 (MQ=255)
                                     tCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGaa                                                      >  2:3276673/1‑93 (MQ=255)
                                      cTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGaaa                                                     >  2:4408313/1‑93 (MQ=255)
                                        gCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAAt                                                   >  2:3069064/1‑93 (MQ=255)
                                            aGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGt                                               <  2:542885/93‑1 (MQ=255)
                                            aGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGt                                               <  2:4494782/93‑1 (MQ=255)
                                            aGGCGGCTGACAAAATTGAAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGt                                               <  1:1965541/93‑1 (MQ=255)
                                             ggcggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGtt                                              >  2:1830377/1‑93 (MQ=255)
                                              gcggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTc                                             <  1:2399330/93‑1 (MQ=255)
                                              gcggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTc                                             <  1:820227/93‑1 (MQ=255)
                                               cggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCa                                            <  1:2351738/93‑1 (MQ=255)
                                                 gcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCaaa                                          <  2:3258828/93‑1 (MQ=255)
                                                 gcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCaaa                                          >  1:2488522/1‑93 (MQ=255)
                                                 gcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCaaa                                          <  2:3654406/93‑1 (MQ=255)
                                                  cTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAg                                         <  2:125489/93‑1 (MQ=255)
                                                  cTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAg                                         <  1:1243186/93‑1 (MQ=255)
                                                  cTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAg                                         <  2:521977/93‑1 (MQ=255)
                                                   tGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAgg                                        >  2:1919237/1‑93 (MQ=255)
                                                    gACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGc                                       <  2:695872/93‑1 (MQ=255)
                                                      cAAAATTGCAATCGTCAACATGGGCCGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGACg                                     <  2:2491915/93‑1 (MQ=255)
                                                      cAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCg                                     >  2:2808287/1‑93 (MQ=255)
                                                      aaaaaTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCg                                     <  1:2979569/92‑1 (MQ=255)
                                                       aaaaTTTCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGt                                    <  1:1931902/93‑1 (MQ=255)
                                                       aaaaTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGt                                    <  1:4697026/93‑1 (MQ=255)
                                                       aaaaTTGCAATCCGCAACGGGTGCAGCCCGTTCTAGGGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTGCAAAGGCCGt                                    <  1:4697025/93‑1 (MQ=255)
                                                        aaaTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTg                                   <  2:1328460/93‑1 (MQ=255)
                                                          aTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGcc                                 >  1:344248/1‑93 (MQ=255)
                                                            tGCAGTCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAg                               <  2:2375021/93‑1 (MQ=255)
                                                                 tCGTCACCATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAGAGGCCGTGCCAGCGAAc                          <  2:4222175/93‑1 (MQ=255)
                                                                 tCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAAc                          >  2:1253550/1‑93 (MQ=255)
                                                                 tCGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAAc                          <  1:2171048/93‑1 (MQ=255)
                                                                  cGTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACt                         <  2:3139690/93‑1 (MQ=255)
                                                                   gTCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTg                        <  1:89592/93‑1 (MQ=255)
                                                                    tCAAGATGGGCCGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACCGc                       <  2:1781823/93‑1 (MQ=255)
                                                                    tCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGc                       >  2:2601811/1‑93 (MQ=255)
                                                                    tCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGc                       <  1:849555/93‑1 (MQ=255)
                                                                    tCAACATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGc                       <  1:1115807/93‑1 (MQ=255)
                                                                      aaCATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCTGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAg                     >  1:3613814/1‑93 (MQ=255)
                                                                      aaCATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAg                     >  1:2263017/1‑93 (MQ=255)
                                                                      aaCATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAg                     <  2:3726834/93‑1 (MQ=255)
                                                                       aCATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGc                    >  1:639391/1‑93 (MQ=255)
                                                                        cATGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCg                   <  2:4402196/93‑1 (MQ=255)
                                                                         aTGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGt                  >  2:763562/1‑93 (MQ=255)
                                                                         aTGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGt                  >  1:2653827/1‑93 (MQ=255)
                                                                         aTGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGt                  >  1:4145545/1‑93 (MQ=255)
                                                                          tGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTa                 >  1:3387084/1‑93 (MQ=255)
                                                                          tGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTa                 >  1:770884/1‑93 (MQ=255)
                                                                          tGGGCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTa                 >  1:4180458/1‑93 (MQ=255)
                                                                            ggCAGCCTGTTCCAGCGGCTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATg               <  2:3147863/93‑1 (MQ=255)
                                                                               aGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGaa            >  2:1622612/1‑93 (MQ=255)
                                                                                gCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGCAAATGAGTTCAAAGGCCGTGCCAGAGAACTGCAGCGTATGGaaa           >  1:4044875/1‑93 (MQ=255)
                                                                                gCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGaaa           >  2:2750546/1‑93 (MQ=255)
                                                                                gCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGaaa           <  2:1272465/93‑1 (MQ=255)
                                                                                 ccTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAAc          >  1:2416512/1‑93 (MQ=255)
                                                                                 ccTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAAc          <  2:4464463/93‑1 (MQ=255)
                                                                                  cTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGCAAcc         <  1:4697844/93‑1 (MQ=255)
                                                                                  cTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAAcc         <  1:4613871/93‑1 (MQ=255)
                                                                                  cTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAAcc         >  2:290579/1‑93 (MQ=255)
                                                                                    gTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGa       >  1:1882130/1‑93 (MQ=255)
                                                                                      tCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATc     <  1:13620/93‑1 (MQ=255)
                                                                                       ccAGCGGGTAGCGCAGAAAGCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAGACCGATCt    <  1:2735088/93‑1 (MQ=255)
                                                                                       ccAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCt    <  1:2735087/93‑1 (MQ=255)
                                                                                        cagcggGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTg   >  1:3879991/1‑93 (MQ=255)
                                                                                         agcggGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGc  >  1:2849260/1‑93 (MQ=255)
                                                                                         agcggGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGc  <  2:2998423/93‑1 (MQ=255)
                                                                                            |                                                                                         
GTTATTAGCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGC  >  NC_000913/200493‑200674

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