Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,510,339 T→C V221A (GTG→GCG)  b1440 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,510,3390TC99.5% 572.1 / ‑5.9 183V221A (GTG→GCG) b1440putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (88/93);  total (90/93)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 189 total aligned reads displayed.

TTGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGT  >  NC_000913/1510249‑1510427
                                                                                          |                                                                                        
ttGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGgcgc                                                                                        >  1:5657460/1‑93 (MQ=255)
   cGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGa                                                                                     <  2:6320393/93‑1 (MQ=255)
   cGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGa                                                                                     <  2:1211431/93‑1 (MQ=255)
   cGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGa                                                                                     >  1:3846024/1‑93 (MQ=255)
    ggACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGaa                                                                                    >  2:3450560/1‑93 (MQ=255)
    ggACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaa                                                                                    <  1:7021458/93‑1 (MQ=255)
     gACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGaaa                                                                                   >  2:5947834/1‑93 (MQ=255)
      aCGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAg                                                                                  <  1:4998585/93‑1 (MQ=255)
      aCGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAg                                                                                  <  1:4591762/93‑1 (MQ=255)
       cgcgGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGt                                                                                 <  2:1081396/93‑1 (MQ=255)
       cGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGt                                                                                 >  1:2691562/1‑93 (MQ=255)
       cGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAAAAGTACCAGGCGGCGCTGAAAGt                                                                                 <  1:1748817/93‑1 (MQ=255)
         cTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgc                                                                               <  1:2530444/93‑1 (MQ=255)
          tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                              >  1:608431/1‑93 (MQ=255)
          tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                              <  1:4488090/93‑1 (MQ=255)
          tGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgct                                                                              >  2:1582051/1‑93 (MQ=255)
           gCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctg                                                                             >  2:9903939/1‑93 (MQ=255)
           gCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGtgctg                                                                             >  1:4916580/1‑93 (MQ=255)
              ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                          >  1:9072253/1‑93 (MQ=255)
              ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                          <  1:1107421/93‑1 (MQ=255)
              ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                          <  2:8089048/93‑1 (MQ=255)
              ttgttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgc                                                                          >  2:7458099/1‑93 (MQ=255)
                gttCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTgcgcgc                                                                        >  2:2245591/1‑93 (MQ=255)
                     ttAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAac                                                                   >  2:6400844/1‑93 (MQ=255)
                      tAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAaca                                                                  <  1:9137221/93‑1 (MQ=255)
                       aaaGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAacac                                                                 <  2:6242123/93‑1 (MQ=255)
                        aaGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAAcaca                                                                <  1:940492/93‑1 (MQ=255)
                          gCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGt                                                              <  2:6616846/93‑1 (MQ=255)
                           ccgcTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtt                                                             <  2:7253448/90‑1 (MQ=255)
                           ccACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGtt                                                             <  1:7665150/93‑1 (MQ=255)
                            cACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGttt                                                            >  1:9912660/1‑93 (MQ=255)
                            cACTCACCGCAGTTGGGCAACAGGCATCCGTATCAACTCCCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTg                                                           <  1:6820479/93‑1 (MQ=255)
                              cTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGa                                                          >  2:4064724/1‑93 (MQ=255)
                              cTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGa                                                          <  1:5775412/93‑1 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         <  2:6881435/93‑1 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         <  2:4152991/93‑1 (MQ=255)
                               tCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGAt                                                         >  2:8721857/1‑93 (MQ=255)
                                cAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATc                                                        >  1:3586972/1‑93 (MQ=255)
                                cAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATc                                                        <  2:9526213/93‑1 (MQ=255)
                                 aGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATcc                                                       <  2:5770080/93‑1 (MQ=255)
                                  gCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCa                                                      <  2:2993436/93‑1 (MQ=255)
                                   ccGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCAt                                                     <  1:5140537/93‑1 (MQ=255)
                                    cGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATc                                                    >  1:9244730/1‑93 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   >  2:2812443/1‑93 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   <  2:1940765/93‑1 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   <  1:7408854/93‑1 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   >  2:8554693/1‑93 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   >  1:3083101/1‑93 (MQ=255)
                                     gCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCg                                                   >  1:5530559/1‑93 (MQ=255)
                                         ttGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                               >  1:4071263/1‑93 (MQ=255)
                                         ttGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTAc                                               >  2:1724415/1‑93 (MQ=255)
                                          tGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACt                                              <  1:4161111/93‑1 (MQ=255)
                                           gggCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTg                                             >  2:4185373/1‑93 (MQ=255)
                                           gggCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTg                                             >  1:9550578/1‑93 (MQ=255)
                                              cATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCa                                          <  1:9493446/93‑1 (MQ=255)
                                               aTCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCAt                                         <  1:5293686/93‑1 (MQ=255)
                                                tCAGCGTTCCGTTTCAACTTTTCGTAGTTCAGTTTCTGGCGGCGCTGTAAGTGCCGCGCGCTCAACACAGTTTGATCTATCGCTATAGGCagg                                        >  1:8992808/1‑91 (MQ=255)
                                                tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                        >  1:7260363/1‑93 (MQ=255)
                                                tCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATg                                        >  2:6798305/1‑93 (MQ=255)
                                                 cAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGa                                       >  1:9226807/1‑93 (MQ=255)
                                                  aGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGac                                      <  2:9917570/93‑1 (MQ=255)
                                                  aGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGac                                      >  1:5446779/1‑93 (MQ=255)
                                                   gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGaca                                     >  1:8749571/1‑93 (MQ=255)
                                                   gCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGaca                                     >  1:2134694/1‑93 (MQ=255)
                                                      aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTTATCCATCCCTACTGGCATGACACtc                                  >  1:5127233/1‑92 (MQ=255)
                                                      aTCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTa                                  >  2:4909003/1‑93 (MQ=255)
                                                       tCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAc                                 <  1:2329200/93‑1 (MQ=255)
                                                       tCCGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAc                                 <  1:9386380/93‑1 (MQ=255)
                                                        ccGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAcc                                <  2:2128012/93‑1 (MQ=255)
                                                        ccGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTAcc                                <  1:1647004/93‑1 (MQ=255)
                                                         cGTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCg                               <  2:5387449/93‑1 (MQ=255)
                                                          gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCgc                              >  2:7283197/1‑92 (MQ=255)
                                                          gTATCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGt                              >  2:5090647/1‑93 (MQ=255)
                                                            aTCAACTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGc                            <  1:4925499/93‑1 (MQ=255)
                                                               aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                         <  2:3085294/93‑1 (MQ=255)
                                                               aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                         <  1:7170777/93‑1 (MQ=255)
                                                               aaCTCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCaaa                         <  1:7018748/93‑1 (MQ=255)
                                                                aCTCACCTAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAAt                        <  2:7632988/93‑1 (MQ=255)
                                                                  tCACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGa                      >  1:5978645/1‑93 (MQ=255)
                                                                   cACCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAg                     >  2:7400214/1‑93 (MQ=255)
                                                                    aCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGc                    <  1:5183811/93‑1 (MQ=255)
                                                                    aCCGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGc                    <  1:689924/93‑1 (MQ=255)
                                                                      cGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGTGa                  >  2:2449995/1‑93 (MQ=255)
                                                                      cGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGa                  >  2:7105026/1‑93 (MQ=255)
                                                                      cGAAGAACAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGa                  <  1:9690535/93‑1 (MQ=255)
                                                                        aagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGGGCGAtt                <  2:9759921/93‑1 (MQ=255)
                                                                        aagaaCAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGAtt                <  1:1075767/93‑1 (MQ=255)
                                                                             cAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaa           >  1:817761/1‑93 (MQ=255)
                                                                             cAGTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaa           <  1:9196479/93‑1 (MQ=255)
                                                                               gTACCAGGCGGCGCTGAAAGTGCTGCGGGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa         <  1:5726534/93‑1 (MQ=255)
                                                                               gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa         >  1:2272242/1‑93 (MQ=255)
                                                                               gTACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCaaaaa         >  1:2071844/1‑93 (MQ=255)
                                                                                tACCAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAAc        <  2:3905808/93‑1 (MQ=255)
                                                                                   cAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAg     <  1:6293818/93‑1 (MQ=255)
                                                                                   cAGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAg     >  1:4811279/1‑93 (MQ=255)
                                                                                   cAGGCGGCGCTGAAAGGTCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAg     <  2:2205488/93‑1 (MQ=255)
                                                                                    aGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACTCTACCGTGCAAATGAGCGATTTCAAAAACGAgg    >  2:2329183/1‑93 (MQ=255)
                                                                                    aGGCGGCGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAgg    >  2:1636835/1‑93 (MQ=255)
                                                                                     ggcggcGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAggg   <  1:6600364/93‑1 (MQ=255)
                                                                                      gcggcgCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGgt  <  2:3291573/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTGCGGACGCTGCGTTGTTCGTTAAAGCCACTCAGCCGCAGTTGGGCATCAGCGATCCGTATCAACTCACCGAAGAACAGTACCAGGCGGTGCTGAAAGTGCTGCGCGCTCAACACAGTTTGATCCATCGCTACTGGCATGACACTACCGTGCAAATGAGCGATTTCAAAAACGAGGGT  >  NC_000913/1510249‑1510427

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