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,0150GA99.4% 541.5 / NA 153G46S (GGT→AGT) b3051putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (97/55);  total (97/56)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ATTATTGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCT  >  NC_000913/3190925‑3191103
                                                                                          |                                                                                        
 ttattGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGtt                                                                                       <  1:348591/93‑1 (MQ=255)
    ttGCCGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAg                                                                                    <  2:2892254/93‑1 (MQ=255)
        cGTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTgg                                                                                >  2:4045708/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               >  1:839218/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               >  1:927989/1‑93 (MQ=255)
         gTATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTggg                                                                               <  1:1463594/93‑1 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTgggg                                                                              >  2:3187558/1‑92 (MQ=255)
          tATGCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGc                                                                              <  1:4030494/93‑1 (MQ=255)
             gCATTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaaa                                                                           >  2:4628411/1‑93 (MQ=255)
               aTTCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCaaaaa                                                                         <  2:1627135/93‑1 (MQ=255)
                ttCTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAg                                                                        >  1:3146047/1‑93 (MQ=255)
                  cTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTg                                                                      >  1:455996/1‑93 (MQ=255)
                  cTGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTg                                                                      <  1:2641954/93‑1 (MQ=255)
                   tGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTgg                                                                     >  1:3479040/1‑93 (MQ=255)
                   tGTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTgg                                                                     >  1:3127451/1‑93 (MQ=255)
                    gTTTATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGt                                                                    >  1:4053557/1‑93 (MQ=255)
                     tttATTATTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTa                                                                   >  1:2488582/1‑93 (MQ=255)
                      ttattaTTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTaa                                                                  >  2:3202135/1‑93 (MQ=255)
                       tattatTGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAAt                                                                 <  2:4303022/93‑1 (MQ=255)
                         ttattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGa                                                               >  2:4263946/1‑93 (MQ=255)
                         ttattGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGa                                                               >  2:1310941/1‑93 (MQ=255)
                           attGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGc                                                             >  2:2651832/1‑93 (MQ=255)
                           attGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGc                                                             >  1:1217856/1‑93 (MQ=255)
                           attGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGc                                                             >  2:4138649/1‑93 (MQ=255)
                            ttGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCg                                                            >  1:3635736/1‑93 (MQ=255)
                             tGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCgg                                                           >  2:2245777/1‑93 (MQ=255)
                             tGGAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCgg                                                           >  2:3074811/1‑93 (MQ=255)
                              ggAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCggg                                                          >  2:81494/1‑92 (MQ=255)
                               gAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTg                                                         >  1:1464310/1‑93 (MQ=255)
                               gAATTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTg                                                         <  1:2549130/93‑1 (MQ=255)
                                aaTTATTTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTgg                                                        >  1:4270261/1‑93 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTATGTGGGCAAAAAGGGGTAATGAGCGGTgga                                                       >  2:1404260/1‑92 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgc                                                       >  1:3279972/1‑93 (MQ=255)
                                 attattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgc                                                       >  2:2140323/1‑93 (MQ=255)
                                  ttattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcg                                                      >  2:2104027/1‑93 (MQ=255)
                                   tattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcgc                                                     >  1:82680/1‑93 (MQ=255)
                                   tattTTCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGgcgc                                                     >  1:3069045/1‑93 (MQ=255)
                                      tttCGCCAGGCTCTATCGTGGCGCTTCGGCAGAGCCAGCTTTTGTTCGCAGTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAAt                                                  >  1:1701058/1‑93 (MQ=255)
                                      tttCGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAAt                                                  <  2:3482489/93‑1 (MQ=255)
                                         cGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGt                                               <  1:3943197/93‑1 (MQ=255)
                                         cGCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGt                                               <  2:1033074/93‑1 (MQ=255)
                                          gCCAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGGGCAATCGTg                                              >  2:362344/1‑93 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGGAAAAAGTGGTAATGAGAGGTGGCGCAATCGTGa                                             >  2:2045652/1‑93 (MQ=255)
                                           ccAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGa                                             >  2:85864/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            <  2:2721146/93‑1 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            <  2:1267756/93‑1 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            >  1:4068879/1‑93 (MQ=255)
                                            cAGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGAt                                            <  1:1800233/93‑1 (MQ=255)
                                             aGGCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATg                                           >  2:3999576/1‑93 (MQ=255)
                                              ggCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGc                                          >  1:1285934/1‑93 (MQ=255)
                                               gCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGcc                                         >  2:684218/1‑93 (MQ=255)
                                               gCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGcc                                         <  2:784410/93‑1 (MQ=255)
                                               gCTCTATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGcc                                         >  1:2311370/1‑93 (MQ=255)
                                                 tctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGGGGTAATGAGCGGTGGCGCAATCGTGATGCCGa                                       >  1:364035/1‑93 (MQ=255)
                                                  ctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGAt                                      >  1:1107114/1‑93 (MQ=255)
                                                  ctATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGAt                                      <  2:2952672/93‑1 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATc                                     >  2:2854848/1‑93 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATc                                     >  2:2583419/1‑93 (MQ=255)
                                                   tATCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATc                                     >  1:3982291/1‑93 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    <  1:671628/93‑1 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    >  1:2129187/1‑93 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    <  2:1425434/93‑1 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCt                                    <  2:4084171/93‑1 (MQ=255)
                                                    aTCGTCGCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGACGATCt                                    >  2:2812770/1‑93 (MQ=255)
                                                     tcgtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCtt                                   >  1:2379872/1‑93 (MQ=255)
                                                     tcgtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCtt                                   >  1:1762490/1‑93 (MQ=255)
                                                       gtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTc                                 <  2:2841101/93‑1 (MQ=255)
                                                       gtcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTc                                 >  2:2713274/1‑93 (MQ=255)
                                                        tcgCGCTTCGGCAGATCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                <  2:3176102/93‑1 (MQ=255)
                                                        tcgCGCTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTGTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCa                                <  2:1767457/93‑1 (MQ=255)
                                                            gcTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaa                            >  2:2353211/1‑93 (MQ=255)
                                                             cTTCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGaaa                           <  2:3699931/93‑1 (MQ=255)
                                                              ttCGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAaat                          >  1:2900438/1‑93 (MQ=255)
                                                                cGGCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGAGGTGGCGCAATCGTGATGCCGATCTTTCATGAaataa                        <  2:2337138/93‑1 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      <  2:4300806/93‑1 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      <  1:2326491/93‑1 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      >  2:4630223/1‑93 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      <  1:2263829/93‑1 (MQ=255)
                                                                  gCAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATc                      <  2:3328394/93‑1 (MQ=255)
                                                                   cAGAGCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcc                     >  1:1197323/1‑93 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   >  1:3664484/1‑93 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   >  2:604869/1‑93 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   >  2:3084438/1‑93 (MQ=255)
                                                                     gagCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATcccc                   >  2:2964448/1‑93 (MQ=255)
                                                                      agCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCa                  <  1:3066936/93‑1 (MQ=255)
                                                                       gCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCAt                 >  1:438770/1‑93 (MQ=255)
                                                                       gCAAGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCAt                 >  2:798081/1‑93 (MQ=255)
                                                                        cAAGCTTTTGTTCGTACTGGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATc                <  1:1222516/93‑1 (MQ=255)
                                                                          aGCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCaa              >  1:4122855/1‑93 (MQ=255)
                                                                           gCTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAat             >  1:1014306/1‑93 (MQ=255)
                                                                            cTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAata            <  1:1018664/93‑1 (MQ=255)
                                                                            cTTTTGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAata            >  2:1142148/1‑93 (MQ=255)
                                                                                tGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGaa        <  1:1786465/93‑1 (MQ=255)
                                                                                tGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGaa        <  1:508554/93‑1 (MQ=255)
                                                                                tGTTCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGaa        <  1:529484/93‑1 (MQ=255)
                                                                                   tCGTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATAc     <  2:3824629/93‑1 (MQ=255)
                                                                                     gTACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATActc   <  2:4477619/93‑1 (MQ=255)
                                                                                      tACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATActct  <  2:4449664/93‑1 (MQ=255)
                                                                                      tACTAGTTTAGGTGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATActct  >  1:2904166/1‑93 (MQ=255)
                                                                                      tACTAGTTTAGGGGGGCAAAAAGTGGTAATGAGCGGTGGCGCAATCGTGATGCCGATCTTTCATGAAATAATCACCATCAATATGAATActct  >  2:2273911/1‑93 (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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: