Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,051 G→A A329T (GCA→ACA)  b1372 → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,0510GA99.6% 851.2 / ‑5.2 230A329T (GCA→ACA) b1372putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (121/108);  total (122/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 245 total aligned reads displayed.

GGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTC  >  NC_000913/1427961‑1428139
                                                                                          |                                                                                        
ggCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACa                                                                                        <  1:4532770/93‑1 (MQ=255)
 gCTCAAAAACAGCGTCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAg                                                                                       <  1:3889677/93‑1 (MQ=255)
 gCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAg                                                                                       >  2:2429757/1‑93 (MQ=255)
 gCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAg                                                                                       >  2:5930343/1‑93 (MQ=255)
 gCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAg                                                                                       >  2:109573/1‑93 (MQ=255)
  cTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGa                                                                                      >  1:3234468/1‑93 (MQ=255)
   tCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGaa                                                                                     <  1:3782223/93‑1 (MQ=255)
    cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGaaa                                                                                    <  2:5637463/93‑1 (MQ=255)
    cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGaaa                                                                                    <  2:454588/93‑1 (MQ=255)
     aaaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAg                                                                                   <  1:6380660/93‑1 (MQ=255)
      aaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGGAAATCGACAGAAAgc                                                                                  >  1:2228322/1‑93 (MQ=255)
       aaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAgcg                                                                                 <  2:6091654/93‑1 (MQ=255)
       aaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAgcg                                                                                 <  1:4486306/93‑1 (MQ=255)
       aaaCAGCGGCTGCGTCGCCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGAAAGAAAgcg                                                                                 <  1:6277361/93‑1 (MQ=255)
         aCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGcgcc                                                                               >  2:5141593/1‑93 (MQ=255)
         aCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGcgcc                                                                               <  1:6005404/93‑1 (MQ=255)
           aGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGcgccgc                                                                             >  2:4648570/1‑93 (MQ=255)
           aGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGcgccgc                                                                             >  1:5556862/1‑93 (MQ=255)
             cGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCAt                                                                           <  1:2653880/93‑1 (MQ=255)
              ggCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATc                                                                          >  1:2953437/1‑93 (MQ=255)
              ggCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATc                                                                          >  1:4768685/1‑93 (MQ=255)
                cTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGt                                                                        <  1:5889757/93‑1 (MQ=255)
                 tGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTc                                                                       <  1:181474/93‑1 (MQ=255)
                 ggCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTc                                                                       <  2:2938910/92‑1 (MQ=255)
                  gCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCt                                                                      >  2:3005131/1‑93 (MQ=255)
                  gCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCt                                                                      >  1:3731727/1‑93 (MQ=255)
                  gCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCt                                                                      <  1:5142086/93‑1 (MQ=255)
                    gtcgtcTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGc                                                                    >  1:1079521/1‑93 (MQ=255)
                       gtcTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTc                                                                 >  1:6040711/1‑93 (MQ=255)
                        tcTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCa                                                                >  2:4231625/1‑93 (MQ=255)
                          tGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAAc                                                              <  2:1934427/93‑1 (MQ=255)
                           cccAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACa                                                             <  1:5527693/92‑1 (MQ=255)
                              aGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGcc                                                          >  2:5295972/1‑93 (MQ=255)
                               gTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCca                                                         >  1:4704956/1‑93 (MQ=255)
                                tGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGTCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                        <  2:4170332/93‑1 (MQ=255)
                                   aGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAAc                                                     <  2:2597628/93‑1 (MQ=255)
                                    gCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACg                                                    >  1:37265/1‑93 (MQ=255)
                                    gCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACATCCACAACg                                                    >  2:5512348/1‑93 (MQ=255)
                                    gCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACg                                                    >  2:1539290/1‑93 (MQ=255)
                                      gTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGaa                                                  <  2:2286983/93‑1 (MQ=255)
                                       tCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAg                                                 >  1:4499271/1‑93 (MQ=255)
                                        caccaaGTGCCGGGCAGCCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAgg                                                <  1:307114/90‑1 (MQ=255)
                                        caacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAgg                                                <  1:1890826/93‑1 (MQ=255)
                                         aacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                               >  1:1232994/1‑93 (MQ=255)
                                         aacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                               >  1:4621497/1‑93 (MQ=255)
                                         aacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                               >  1:3552498/1‑93 (MQ=255)
                                         aacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                               <  2:3845715/93‑1 (MQ=255)
                                         aacaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                               >  1:2386712/1‑93 (MQ=255)
                                          acaaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCt                                              <  1:851110/93‑1 (MQ=255)
                                           caaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTg                                             <  2:6184858/93‑1 (MQ=255)
                                           caaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTg                                             <  2:5192488/93‑1 (MQ=255)
                                           caaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTg                                             <  2:275076/93‑1 (MQ=255)
                                            aaGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTgg                                            >  2:4846110/1‑93 (MQ=255)
                                            aaGGGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTgg                                            <  2:2535366/93‑1 (MQ=255)
                                             aGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGc                                           >  2:2033112/1‑93 (MQ=255)
                                               tGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGa                                         >  2:4981448/1‑93 (MQ=255)
                                                  cGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGc                                      <  1:1965928/93‑1 (MQ=255)
                                                  cGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGc                                      <  2:6422646/93‑1 (MQ=255)
                                                   gggCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGcc                                     <  1:5619149/93‑1 (MQ=255)
                                                    ggCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGccc                                    >  1:1163156/1‑92 (MQ=255)
                                                    ggCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCa                                    >  1:425937/1‑93 (MQ=255)
                                                    ggCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCa                                    <  2:3519177/93‑1 (MQ=255)
                                                     gCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCATCAGCCACAACGAAGGCTGGCGAAGCCAc                                   <  1:5791926/93‑1 (MQ=255)
                                                        ggCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGa                                >  2:5332968/1‑93 (MQ=255)
                                                          ccTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAAc                              >  1:2605539/1‑93 (MQ=255)
                                                           cTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACa                             >  2:1274270/1‑93 (MQ=255)
                                                           cTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACa                             <  2:477181/93‑1 (MQ=255)
                                                           cTCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACa                             <  1:129916/93‑1 (MQ=255)
                                                            tCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAg                            >  2:1798023/1‑93 (MQ=255)
                                                            tCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAg                            <  1:745351/93‑1 (MQ=255)
                                                            tCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAg                            >  1:2054027/1‑93 (MQ=255)
                                                            tCAGCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAg                            >  2:4206361/1‑93 (MQ=255)
                                                            ccAGCCAGTGCCCCCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAg                            <  2:1011901/92‑1 (MQ=255)
                                                               gCCAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGcc                         <  2:3991597/93‑1 (MQ=255)
                                                                ccAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCa                        <  2:3501036/93‑1 (MQ=255)
                                                                 cAGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAg                       >  1:336976/1‑93 (MQ=255)
                                                                  aGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgc                      >  1:1437560/1‑93 (MQ=255)
                                                                  aGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgc                      >  2:701224/1‑93 (MQ=255)
                                                                  aGTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgc                      >  2:5296594/1‑93 (MQ=255)
                                                                   gTGCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgcg                     >  2:2170599/1‑93 (MQ=255)
                                                                     gCCACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgca                   >  2:377935/1‑93 (MQ=255)
                                                                       cACCGCCGCCGGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcagc                 <  1:3329503/93‑1 (MQ=255)
                                                                       cACCGCCGCCGAAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcagc                 <  1:826858/93‑1 (MQ=255)
                                                                          cgccgccgGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCagcagc              <  2:4299068/93‑1 (MQ=255)
                                                                           gccgccgGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCg             >  1:4183268/1‑93 (MQ=255)
                                                                            ccgccgAAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGa            <  2:3324696/93‑1 (MQ=255)
                                                                              gccgGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAgg          >  1:2305857/1‑93 (MQ=255)
                                                                               ccgGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGt         <  1:5917735/93‑1 (MQ=255)
                                                                                cgGAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTc        <  1:6456742/93‑1 (MQ=255)
                                                                                 ggAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       >  2:4512653/1‑93 (MQ=255)
                                                                                 ggAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       >  1:5088649/1‑93 (MQ=255)
                                                                                 ggAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       >  1:5305815/1‑93 (MQ=255)
                                                                                 ggAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       >  2:2310338/1‑93 (MQ=255)
                                                                                 ggAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       >  2:2324250/1‑93 (MQ=255)
                                                                                 gaaaaaTCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt       <  2:5699152/91‑1 (MQ=255)
                                                                                  gAAAATCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTg      >  1:3986556/1‑93 (MQ=255)
                                                                                   aaaaGCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGc     <  2:3697573/93‑1 (MQ=255)
                                                                                    aaaTCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCt    >  1:48657/1‑93 (MQ=255)
                                                                                     aaTCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCtt   >  2:5252433/1‑93 (MQ=255)
                                                                                     aaTCGACAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGACCCTGCACAGGGCAGCGCAGCAGCGAGGTCTGCtt   >  2:3170830/1‑93 (MQ=255)
                                                                                          |                                                                                        
GGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTC  >  NC_000913/1427961‑1428139

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑