Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,191,015 G→A G46S (GGT→AGT)  b3051 → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,191,0150GA98.2% 999.0 / ‑1.0 283G46S (GGT→AGT) b3051putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (161/117);  total (165/118)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 290 total aligned reads displayed.

ATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTG  >  NC_000913/3190925‑3191104
                                                                                          |                                                                                         
 ttattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGtt                                                                                        <  2:6443256/93‑1 (MQ=255)
  tattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGttt                                                                                       >  1:1352328/1‑93 (MQ=255)
   attGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTa                                                                                      <  1:7102018/93‑1 (MQ=255)
    ttGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAg                                                                                     >  1:6462950/1‑93 (MQ=255)
     tGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAgg                                                                                    <  2:6871996/93‑1 (MQ=255)
     tGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAgg                                                                                    <  1:3258568/93‑1 (MQ=255)
     tGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAgg                                                                                    <  2:285009/93‑1 (MQ=255)
      gCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGt                                                                                   >  2:2389316/1‑93 (MQ=255)
        cGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTgg                                                                                 <  1:5366090/93‑1 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                                >  1:5841509/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                                >  1:3509505/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGCACTAGTTTAGGTggg                                                                                >  1:1792838/1‑93 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTgggt                                                                               >  1:6713639/1‑92 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGc                                                                               >  2:3913455/1‑93 (MQ=255)
            tGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaa                                                                             >  2:3576985/1‑93 (MQ=255)
            tGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaa                                                                             >  1:4174118/1‑93 (MQ=255)
            tGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaa                                                                             >  1:7510965/1‑93 (MQ=255)
               aTTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaaaaa                                                                          >  2:6014339/1‑93 (MQ=255)
                ttCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAg                                                                         >  1:495760/1‑93 (MQ=255)
                ttCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAg                                                                         >  2:5012453/1‑93 (MQ=255)
                ttCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAg                                                                         >  1:5676724/1‑93 (MQ=255)
                 tCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAgg                                                                        >  2:265173/1‑92 (MQ=255)
                  cTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTg                                                                       >  2:1631068/1‑93 (MQ=255)
                    gTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGt                                                                     >  1:820468/1‑93 (MQ=255)
                      ttattaTTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGCACTAGTTTAGGTGGGCAAAAAGTGGTaa                                                                   >  1:2003723/1‑93 (MQ=255)
                       tattatTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAAt                                                                  >  2:2719046/1‑93 (MQ=255)
                       tattatTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAAt                                                                  >  2:4583749/1‑93 (MQ=255)
                        attattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATg                                                                 >  1:5481193/1‑93 (MQ=255)
                        atagttGGAAGTATTTTCGCCAGGCGCTGTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATg                                                                 <  2:6316812/89‑1 (MQ=255)
                         ttattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGa                                                                >  1:1190157/1‑93 (MQ=255)
                         ttattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGa                                                                >  2:1990872/1‑93 (MQ=255)
                          tattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAg                                                               >  2:808965/1‑93 (MQ=255)
                          tattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAg                                                               >  1:108900/1‑93 (MQ=255)
                          tattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAg                                                               <  1:2119124/93‑1 (MQ=255)
                             tGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCgg                                                            >  2:1886451/1‑93 (MQ=255)
                              ggAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGt                                                           <  2:6799566/93‑1 (MQ=255)
                               gAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTg                                                          <  2:7544623/93‑1 (MQ=255)
                                aaTTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTgg                                                         <  1:4774941/93‑1 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTggg                                                        >  2:7382868/1‑92 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgc                                                        >  2:2446909/1‑93 (MQ=255)
                                  ttattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcg                                                       >  1:7339310/1‑93 (MQ=255)
                                  ttattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcg                                                       <  1:2911780/93‑1 (MQ=255)
                                     ttttCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCaa                                                    >  2:2953804/1‑93 (MQ=255)
                                      tttCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAAt                                                   >  2:735055/1‑93 (MQ=255)
                                       ttCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATc                                                  <  2:5351613/93‑1 (MQ=255)
                                         cGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGt                                                <  1:1936222/93‑1 (MQ=255)
                                          gCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTg                                               >  1:7686708/1‑93 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                              <  2:6245405/93‑1 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                              <  2:3325457/93‑1 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                             >  1:4065071/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                             <  1:5290013/93‑1 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                             >  1:4489785/1‑93 (MQ=255)
                                             aGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATg                                            >  2:594121/1‑93 (MQ=255)
                                              ggCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGc                                           >  1:2706156/1‑93 (MQ=255)
                                                 tctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGa                                        >  1:1676653/1‑93 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATc                                      >  1:2475477/1‑93 (MQ=255)
                                                     tcgtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCtt                                    >  1:517834/1‑93 (MQ=255)
                                                       gtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTc                                  >  2:5826158/1‑93 (MQ=255)
                                                        tcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                 <  1:3108027/93‑1 (MQ=255)
                                                        tcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                 >  1:5957028/1‑93 (MQ=255)
                                                        tcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                 >  1:6804938/1‑93 (MQ=255)
                                                         cgcgCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCAt                                >  1:4662398/1‑93 (MQ=255)
                                                          gcgcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATg                               <  1:3012862/93‑1 (MQ=255)
                                                           cgcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGa                              <  1:5086943/93‑1 (MQ=255)
                                                            gcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaa                             >  1:6431693/1‑93 (MQ=255)
                                                            gcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaa                             <  2:4194893/93‑1 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                            >  1:239645/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                            >  1:2039071/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                            <  1:4906758/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataat                        >  2:3681006/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                       >  1:6932985/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                       >  1:3333884/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                       >  1:7219367/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                       <  1:1032938/93‑1 (MQ=255)
                                                                   cAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcc                      <  2:4266054/93‑1 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                    <  1:4425951/93‑1 (MQ=255)
                                                                      agCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCa                   >  2:1273515/1‑93 (MQ=255)
                                                                      agCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCa                   <  2:5165052/93‑1 (MQ=255)
                                                                       gCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCAt                  >  2:5599765/1‑93 (MQ=255)
                                                                        cAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATc                 <  2:2612032/93‑1 (MQ=255)
                                                                        cAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATc                 <  1:1575254/93‑1 (MQ=255)
                                                                         aaGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCa                <  2:4062417/93‑1 (MQ=255)
                                                                         aaGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCa                >  2:2115727/1‑93 (MQ=255)
                                                                            cTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAata             <  1:3575745/93‑1 (MQ=255)
                                                                            cTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAata             <  2:6470738/93‑1 (MQ=255)
                                                                             ttttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAatat            <  2:471812/93‑1 (MQ=255)
                                                                             ttttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAatat            >  1:3376080/1‑93 (MQ=255)
                                                                             ttttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAatat            >  1:3770918/1‑93 (MQ=255)
                                                                              tttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATg           >  2:5459366/1‑93 (MQ=255)
                                                                              tttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATg           >  1:286601/1‑93 (MQ=255)
                                                                               ttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGa          <  2:5523677/93‑1 (MQ=255)
                                                                                tGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGaa         >  1:3469512/1‑93 (MQ=255)
                                                                                 gTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAAt        >  1:275259/1‑93 (MQ=255)
                                                                                  ttCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATa       <  1:7278550/93‑1 (MQ=255)
                                                                                  ttCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATa       <  1:2214681/93‑1 (MQ=255)
                                                                                  ttCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATa       >  1:6963477/1‑93 (MQ=255)
                                                                                  ttCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATa       <  2:1122825/93‑1 (MQ=255)
                                                                                   tCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATAc      <  2:4528991/93‑1 (MQ=255)
                                                                                   tCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATAc      <  2:7419008/93‑1 (MQ=255)
                                                                                     gTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATActc    <  1:1463863/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTG  >  NC_000913/3190925‑3191104

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