Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,511,365 A→G L173P (CTT→CCT)  fecB ← FecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,3650AG99.5% 617.7 / NA 198L173P (CTT→CCT) fecBFecB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (110/87);  total (110/88)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

GCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATA  >  NC_000913/4511275‑4511454
                                                                                          |                                                                                         
gCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgc                                                                                         <  2:646826/93‑1 (MQ=255)
 ccAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgct                                                                                        <  2:1719062/93‑1 (MQ=255)
     cACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgatggct                                                                                    >  1:3216210/1‑93 (MQ=255)
      aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                   >  1:3875053/1‑93 (MQ=255)
      aCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGgctggctg                                                                                   >  1:1266257/1‑93 (MQ=255)
          tGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGccc                                                                               >  2:4296807/1‑93 (MQ=255)
           gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCa                                                                              >  1:435600/1‑93 (MQ=255)
             cGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACt                                                                            >  1:2994016/1‑93 (MQ=255)
             cGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACt                                                                            <  2:196468/93‑1 (MQ=255)
               gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                          >  2:4496032/1‑93 (MQ=255)
               gTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgc                                                                          <  1:393072/93‑1 (MQ=255)
                tCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgcg                                                                         <  1:1023473/93‑1 (MQ=255)
                tCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgcg                                                                         >  2:3622348/1‑93 (MQ=255)
                 ccAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTgcgc                                                                        >  1:4598493/1‑93 (MQ=255)
                    ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCAt                                                                     >  1:407989/1‑93 (MQ=255)
                      tctcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATcc                                                                   <  1:4558925/93‑1 (MQ=255)
                      tctcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATcc                                                                   <  1:2790311/93‑1 (MQ=255)
                       ctcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCct                                                                  <  2:1479921/93‑1 (MQ=255)
                        tcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctc                                                                 <  2:2295795/93‑1 (MQ=255)
                        tcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctc                                                                 >  1:3431054/1‑93 (MQ=255)
                        tcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctc                                                                 <  2:1507701/93‑1 (MQ=255)
                        tcCCGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATGctc                                                                 >  2:3656453/1‑93 (MQ=255)
                          cTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctc                                                               <  1:2678007/93‑1 (MQ=255)
                          cTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctc                                                               <  2:1822198/93‑1 (MQ=255)
                           tGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctct                                                              >  1:1862176/1‑93 (MQ=255)
                           tGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCctctct                                                              >  1:4257832/1‑93 (MQ=255)
                            gAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCCCTGCGCCATCCTCTCtt                                                             >  2:3988695/1‑93 (MQ=255)
                             aGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCttt                                                            <  1:3949121/93‑1 (MQ=255)
                              gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCCCGCGTGTCCCTTTGGGGGGCTGGCTGGCCCCCTGCGCCATCCCCTCtttt                                                           >  1:3312699/1‑92 (MQ=255)
                              gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCCTCCTCTCTTTa                                                           >  1:2471169/1‑93 (MQ=255)
                              gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCtttt                                                           >  2:4179923/1‑92 (MQ=255)
                               tATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAt                                                          >  2:844570/1‑93 (MQ=255)
                               tATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAt                                                          >  2:2581152/1‑93 (MQ=255)
                                aTGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtg                                                         >  2:454254/1‑93 (MQ=255)
                                aTGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtg                                                         <  1:295093/93‑1 (MQ=255)
                                 tGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtgt                                                        <  2:4660854/93‑1 (MQ=255)
                                 tGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtgt                                                        >  2:1245176/1‑93 (MQ=255)
                                  gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtgtt                                                       >  1:2516747/1‑93 (MQ=255)
                                  gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtgtt                                                       <  1:1589533/93‑1 (MQ=255)
                                    aGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTAtgttgt                                                     >  2:282998/1‑93 (MQ=255)
                                     ggTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATgttgtt                                                    >  2:2948716/1‑93 (MQ=255)
                                       ttGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGTGCCATCCTCTCTTTATGTTGTTcc                                                  >  2:1605632/1‑93 (MQ=255)
                                       ttGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTcc                                                  <  2:3199981/93‑1 (MQ=255)
                                       ttGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTcc                                                  <  1:1651594/93‑1 (MQ=255)
                                        tGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCa                                                 <  1:614560/93‑1 (MQ=255)
                                         gAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAg                                                >  2:1912736/1‑93 (MQ=255)
                                         gAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAg                                                >  2:3884445/1‑93 (MQ=255)
                                          aaTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGa                                               <  2:2082674/93‑1 (MQ=255)
                                           aTTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGAc                                              >  1:1179705/1‑93 (MQ=255)
                                            ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACg                                             >  1:3605739/1‑93 (MQ=255)
                                            ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACg                                             >  1:4486186/1‑93 (MQ=255)
                                              gCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTg                                           >  2:555050/1‑93 (MQ=255)
                                                tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGcc                                         >  1:720087/1‑93 (MQ=255)
                                                tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGcc                                         >  1:1910576/1‑93 (MQ=255)
                                                tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGcc                                         >  2:2646430/1‑93 (MQ=255)
                                                 gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCt                                        >  1:3131212/1‑93 (MQ=255)
                                                 gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCt                                        <  2:2434682/93‑1 (MQ=255)
                                                 gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCt                                        >  1:1992166/1‑93 (MQ=255)
                                                  ttCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTg                                       >  2:3407894/1‑93 (MQ=255)
                                                  ttCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTg                                       >  2:4222590/1‑93 (MQ=255)
                                                  ttCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTg                                       >  1:3651225/1‑93 (MQ=255)
                                                    cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCa                                     <  1:4348506/93‑1 (MQ=255)
                                                    cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCa                                     <  1:311709/93‑1 (MQ=255)
                                                     gcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAt                                    >  2:2803124/1‑93 (MQ=255)
                                                      cgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtc                                   >  1:1975968/1‑93 (MQ=255)
                                                      cgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtc                                   >  1:775399/1‑93 (MQ=255)
                                                       gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtct                                  >  1:4469823/1‑93 (MQ=255)
                                                        tGATGTGCCAAAGGCCACGCGTGTGCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtctc                                 >  2:4517993/1‑93 (MQ=255)
                                                        tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtctc                                 >  2:1558209/1‑93 (MQ=255)
                                                          aTGTGCTAAAGGCCACGCGTGTCCCTTTGGGAGGCTTGCTGGCCCCCTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtctctc                               <  2:2615760/93‑1 (MQ=255)
                                                          aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCAtctctc                               >  1:1006278/1‑93 (MQ=255)
                                                            gtgCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGc                             <  1:3908809/93‑1 (MQ=255)
                                                            gtgCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGc                             >  1:4173687/1‑93 (MQ=255)
                                                            gtgCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGc                             <  1:766325/93‑1 (MQ=255)
                                                            gtgCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGc                             <  2:2893884/93‑1 (MQ=255)
                                                              gCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCtt                           >  1:1077867/1‑93 (MQ=255)
                                                                cAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCTCtttt                         >  2:1188803/1‑93 (MQ=255)
                                                                cAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCtttt                         >  1:959565/1‑93 (MQ=255)
                                                                 aaaGGCCGCGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCttttt                        <  1:229114/93‑1 (MQ=255)
                                                                 aaaGGCCACGCGTGTCCCTTTGGGAGGTTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCttttt                        <  1:827830/93‑1 (MQ=255)
                                                                 aaaGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCttttt                        >  2:34807/1‑93 (MQ=255)
                                                                  aaGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTa                       <  1:353103/93‑1 (MQ=255)
                                                                   aGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAc                      >  2:3820524/1‑93 (MQ=255)
                                                                   aGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAc                      >  1:1348035/1‑93 (MQ=255)
                                                                    ggCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTAcc                     >  2:1934273/1‑93 (MQ=255)
                                                                     gCCACCCGTTTCCCTTTGGGAGGCTGGCTGGGCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGCTGTGCCTGCATCTCTCGCttttttact                    >  2:3943963/1‑89 (MQ=255)
                                                                      ccACGCGTGTCCCTTTGGGAGTCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACcca                   <  2:756014/93‑1 (MQ=255)
                                                                      ccACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACcca                   <  2:28192/93‑1 (MQ=255)
                                                                        aCGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacc                 >  1:947596/1‑93 (MQ=255)
                                                                         cgcgTGTCCCTTTGGGAGGCTGGCGGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACccacca                <  1:1215572/93‑1 (MQ=255)
                                                                          gcgTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAt               <  2:2456043/93‑1 (MQ=255)
                                                                           cgTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCAtt              <  2:375379/93‑1 (MQ=255)
                                                                               tCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGc          <  1:2211265/93‑1 (MQ=255)
                                                                                cccTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc         <  2:250736/93‑1 (MQ=255)
                                                                                cccTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc         >  1:3481276/1‑93 (MQ=255)
                                                                                  cTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCga       <  2:3277730/93‑1 (MQ=255)
                                                                                     tGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatga    <  1:214532/93‑1 (MQ=255)
                                                                                      gggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat   <  1:1322394/93‑1 (MQ=255)
                                                                                      gggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat   >  2:2135891/1‑93 (MQ=255)
                                                                                       ggAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATa  <  2:1009729/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATA  >  NC_000913/4511275‑4511454

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