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,0150GA95.5% 521.6 / 6.2 156G46S (GGT→AGT) b3051putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base A (80/69);  total (86/70)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 160 total aligned reads displayed.

ATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCT  >  NC_000913/3190925‑3191103
                                                                                          |                                                                                        
attattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGt                                                                                        <  2:2052463/93‑1 (MQ=255)
attattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGt                                                                                        >  2:2770264/1‑93 (MQ=255)
 ttattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGtt                                                                                       <  1:2742997/93‑1 (MQ=255)
  tattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGttt                                                                                      >  1:4476902/1‑93 (MQ=255)
   attGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTa                                                                                     >  1:438744/1‑93 (MQ=255)
   attGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTa                                                                                     <  1:3195996/93‑1 (MQ=255)
   attGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTa                                                                                     <  1:1221058/93‑1 (MQ=255)
    ttGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAg                                                                                    >  1:3827567/1‑93 (MQ=255)
    ttGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAg                                                                                    <  1:597131/93‑1 (MQ=255)
      gCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGt                                                                                  <  2:3310543/93‑1 (MQ=255)
       ccGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTg                                                                                 <  2:490197/93‑1 (MQ=255)
       ccGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTg                                                                                 <  2:2413324/93‑1 (MQ=255)
        cGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTgg                                                                                <  2:296189/93‑1 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               >  2:589967/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               <  1:2457627/93‑1 (MQ=255)
         gTATGCATTCTGTTTATTACTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               <  1:1628006/93‑1 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGc                                                                              >  2:1011070/1‑93 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGc                                                                              >  1:3370126/1‑93 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGc                                                                              >  1:397381/1‑93 (MQ=255)
           aTGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCa                                                                             <  2:2576423/93‑1 (MQ=255)
            tGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCaa                                                                            >  1:3069722/1‑93 (MQ=255)
            tGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaa                                                                            >  2:1306481/1‑93 (MQ=255)
             gCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCATTTGTTCGTACTAGTTTAGGTGGGCaaa                                                                           >  1:364973/1‑93 (MQ=255)
              cATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaaaa                                                                          <  2:1097287/93‑1 (MQ=255)
                 tCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGt                                                                       <  2:1743341/93‑1 (MQ=255)
                  cTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTg                                                                      >  1:3514991/1‑93 (MQ=255)
                    gTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGt                                                                    >  1:3622927/1‑93 (MQ=255)
                     tttATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTa                                                                   >  1:335561/1‑93 (MQ=255)
                     tttATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTa                                                                   <  2:829392/93‑1 (MQ=255)
                     tttATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTa                                                                   >  1:1395452/1‑93 (MQ=255)
                        attattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATg                                                                >  1:4392645/1‑93 (MQ=255)
                          tattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAg                                                              <  1:4703151/93‑1 (MQ=255)
                           attGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGc                                                             <  1:2990467/93‑1 (MQ=255)
                            ttGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGCGGGCAAAAAGTGGTAATGAGCg                                                            <  2:2993065/93‑1 (MQ=255)
                               gAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTg                                                         <  1:2115714/93‑1 (MQ=255)
                               gAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTg                                                         >  2:4117154/1‑93 (MQ=255)
                                aaTTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTgg                                                        >  1:1902361/1‑93 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgc                                                       <  2:1615948/93‑1 (MQ=255)
                                   tattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcgc                                                     <  2:3884441/93‑1 (MQ=255)
                                   tattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcgc                                                     >  1:1074364/1‑93 (MQ=255)
                                    attTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGATCGGTGgcgga                                                    >  2:245914/1‑91 (MQ=255)
                                       ttCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATc                                                 >  1:794524/1‑93 (MQ=255)
                                         cGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGt                                               <  2:3019514/93‑1 (MQ=255)
                                         cGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGt                                               <  2:3019513/93‑1 (MQ=255)
                                          gCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTg                                              >  1:3838797/1‑93 (MQ=255)
                                          gCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGAGGTAATGAGCGGTGGCGCAATCGTg                                              >  1:3838707/1‑93 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                             <  1:3773827/93‑1 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                             <  2:4271470/93‑1 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                             >  1:3480833/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            >  1:4167807/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            >  2:3917211/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            >  1:643625/1‑93 (MQ=255)
                                             aGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATg                                           >  1:1082592/1‑93 (MQ=255)
                                             aGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATg                                           >  1:1916917/1‑93 (MQ=255)
                                              ggCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAGAAGTGGTAATGAGCGGTGGCGCAATCGTGATGc                                          >  1:4568697/1‑93 (MQ=255)
                                              ggCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGc                                          >  1:1047751/1‑93 (MQ=255)
                                                ctctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCg                                        >  1:2551193/1‑93 (MQ=255)
                                                 tctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGa                                       >  1:4677286/1‑93 (MQ=255)
                                                 tctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGa                                       >  2:1185519/1‑93 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGAtt                                     >  1:244126/1‑92 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATc                                     >  1:3737162/1‑93 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    <  1:297431/93‑1 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    >  1:1315746/1‑93 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    >  1:2793549/1‑93 (MQ=255)
                                                       gtcgCGATTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTTGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTc                                 >  2:4651223/1‑93 (MQ=255)
                                                        tcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                <  2:4425977/93‑1 (MQ=255)
                                                         cgcgCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCAt                               <  2:4324240/93‑1 (MQ=255)
                                                          gcgcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATg                              >  2:3358491/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                           >  1:2569441/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                           >  2:2813728/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                           <  1:4179810/93‑1 (MQ=255)
                                                              ttCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaat                          >  1:1169445/1‑93 (MQ=255)
                                                              ttCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaat                          <  2:2850082/93‑1 (MQ=255)
                                                              ttCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaat                          >  1:993402/1‑93 (MQ=255)
                                                              ttCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaat                          >  2:2639828/1‑93 (MQ=255)
                                                               tCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaata                         <  2:1592091/93‑1 (MQ=255)
                                                               tCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaata                         <  2:1109618/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGGaataat                       <  1:895758/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataat                       <  1:4106025/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataat                       <  2:2127237/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataat                       <  1:1798796/93‑1 (MQ=255)
                                                                 ggCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataat                       >  1:3601018/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      >  1:2662953/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      >  2:3789533/1‑93 (MQ=255)
                                                                    agagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATccc                    <  1:807554/93‑1 (MQ=255)
                                                                    agagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGGGATTCCGATCTTTCATGAAATAATctc                    >  1:485404/1‑91 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   <  1:4661750/93‑1 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   >  1:3894010/1‑93 (MQ=255)
                                                                      agCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCa                  <  2:547884/93‑1 (MQ=255)
                                                                       gCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCAt                 >  1:2223382/1‑93 (MQ=255)
                                                                         aaGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCa               >  1:3305252/1‑93 (MQ=255)
                                                                          aGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCaa              >  2:1399877/1‑93 (MQ=255)
                                                                          aGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCaa              <  1:2851385/93‑1 (MQ=255)
                                                                            cTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGTCGATCTTTCATGAAATAATCCCCATCAata            <  1:1776864/93‑1 (MQ=255)
                                                                              tttGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATg          <  1:2393291/93‑1 (MQ=255)
                                                                                 gTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAAt       >  1:1554216/1‑93 (MQ=255)
                                                                                 gTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAAt       <  2:4143346/93‑1 (MQ=255)
                                                                                   tCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATAc     <  1:1851312/93‑1 (MQ=255)
                                                                                   tCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATAc     <  2:564444/93‑1 (MQ=255)
                                                                                      tACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATActct  <  2:2324814/93‑1 (MQ=255)
                                                                                          |                                                                                        
ATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCT  >  NC_000913/3190925‑3191103

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