Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,013,176 C→G G182R (GGA→CGA)  ysgA ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,013,1760CG99.7% 2052.7 / ‑5.7 596G182R (GGA→CGA) ysgAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (314/280);  minor base T (2/0);  total (316/280)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 620 total aligned reads displayed.

CATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCGGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013090‑4013260
                                                                                      |                                                                                     
cATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACACCCATATACGCCGAGAATCGGGCggt                                                                                >  1:5917128/1‑91 (MQ=255)
  ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTa                                                                               >  1:122517/1‑93 (MQ=255)
  ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTa                                                                               <  1:6065514/93‑1 (MQ=255)
   tAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGCGAATCGGCGCGTTaa                                                                              <  1:8587977/93‑1 (MQ=255)
   tAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                              >  1:5765094/1‑93 (MQ=255)
   tAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                              <  2:1825358/93‑1 (MQ=255)
   tAGCAGCCCGCAGCGCCTGGCCCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                              <  2:4877587/93‑1 (MQ=255)
    agcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAg                                                                             <  1:4855537/93‑1 (MQ=255)
      cagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGAt                                                                           <  2:363060/93‑1 (MQ=255)
       agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          >  1:2033530/1‑93 (MQ=255)
       agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          >  1:1596741/1‑93 (MQ=255)
       agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          <  2:5906738/93‑1 (MQ=255)
       agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          >  2:6049470/1‑93 (MQ=255)
       agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          <  2:6226154/93‑1 (MQ=255)
        gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  2:919328/1‑93 (MQ=255)
        gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         <  1:8560980/93‑1 (MQ=255)
        gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  1:6412568/1‑93 (MQ=255)
        gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  1:5600004/1‑93 (MQ=255)
        gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  1:2933904/1‑93 (MQ=255)
            gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTg                                                                     >  2:1021408/1‑93 (MQ=255)
            gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTg                                                                     <  2:3389293/93‑1 (MQ=255)
             cAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGc                                                                    >  1:6720811/1‑93 (MQ=255)
               gcgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGa                                                                  >  2:3897996/1‑93 (MQ=255)
                cgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGATCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCgtt                                                                 >  2:2935603/1‑91 (MQ=255)
                cgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGat                                                                 <  2:124247/93‑1 (MQ=255)
                    tGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCa                                                             >  2:5552041/1‑93 (MQ=255)
                      gcggATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAGATCGGTTGCGATATCAACa                                                          <  2:5130479/93‑1 (MQ=255)
                      gcgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAAc                                                           <  1:8364780/93‑1 (MQ=255)
                       cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACa                                                          <  1:276411/93‑1 (MQ=255)
                       cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACa                                                          <  2:8756715/93‑1 (MQ=255)
                       cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACa                                                          <  2:4636299/93‑1 (MQ=255)
                       cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACa                                                          <  1:7553656/93‑1 (MQ=255)
                        gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                         >  2:2363813/1‑93 (MQ=255)
                         cATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAgg                                                        >  2:7161060/1‑93 (MQ=255)
                         cATGGTTTCAACGCTCTCCTGCGGAAGTCTGCTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAgg                                                        <  1:5240573/93‑1 (MQ=255)
                          aTGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgt                                                       >  1:1351816/1‑93 (MQ=255)
                           tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      <  1:156967/93‑1 (MQ=255)
                           tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  1:5714768/1‑93 (MQ=255)
                           tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  1:3169718/1‑93 (MQ=255)
                           tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  1:7950000/1‑93 (MQ=255)
                           tGGTTTCAACGCTCTCCTGCGGAATGCTGTTAGCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      <  1:3065021/93‑1 (MQ=255)
                            ggTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgttg                                                     >  1:7103759/1‑93 (MQ=255)
                                   aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              <  2:2866664/93‑1 (MQ=255)
                                    cGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTg                                             <  1:2405089/93‑1 (MQ=255)
                                     gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                                            >  1:4592294/1‑93 (MQ=255)
                                     gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                                            >  1:5454430/1‑93 (MQ=255)
                                      ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                           >  2:8119186/1‑93 (MQ=255)
                                      ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                           <  2:3724429/93‑1 (MQ=255)
                                        ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  2:694782/93‑1 (MQ=255)
                                        ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  1:4751973/93‑1 (MQ=255)
                                        ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  1:369387/93‑1 (MQ=255)
                                         tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        <  2:4097459/93‑1 (MQ=255)
                                          ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       >  2:4494190/1‑93 (MQ=255)
                                           cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCac                                      >  1:7443223/1‑92 (MQ=255)
                                           cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                      <  2:8679051/93‑1 (MQ=255)
                                           cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                      >  2:2446877/1‑93 (MQ=255)
                                            tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                     >  1:2303248/1‑93 (MQ=255)
                                            tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                     <  2:1663120/93‑1 (MQ=255)
                                             gCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                    >  1:3396950/1‑93 (MQ=255)
                                              cGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGa                                   <  1:1140018/93‑1 (MQ=255)
                                                 aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                >  1:103925/1‑93 (MQ=255)
                                                 aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                <  1:5083414/93‑1 (MQ=255)
                                                 aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                <  1:772401/93‑1 (MQ=255)
                                                 aaTGCTGCTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                <  1:8442569/93‑1 (MQ=255)
                                                   tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              >  1:8161082/1‑93 (MQ=255)
                                                    gCTGATATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             <  2:2267304/93‑1 (MQ=255)
                                                     cTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                            <  2:1085773/93‑1 (MQ=255)
                                                     cTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                            >  1:7226262/1‑93 (MQ=255)
                                                       ttgaTCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGcc                          <  1:8646370/90‑1 (MQ=255)
                                                       gTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGCCTTGTCGcc                          >  1:4233731/1‑93 (MQ=255)
                                                       gTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGcc                          >  1:1991420/1‑93 (MQ=255)
                                                       gTTATCCTGACCACCATATAAGACGAGAATCGGCGCGCTAAGATCGGATGCGATATCAACAGGTTGTTTCGGTGAATTTAGAGACTTGTCggc                          >  2:3281926/1‑91 (MQ=255)
                                                         tATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg                        >  2:3118607/1‑93 (MQ=255)
                                                          aTCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGt                       >  1:1752983/1‑93 (MQ=255)
                                                            ccTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCa                     <  1:4372214/93‑1 (MQ=255)
                                                            ccTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGGTGTTTCGGTGAATTCAGCGACTTGTCGCCCGtaa                     >  2:2675096/1‑91 (MQ=255)
                                                               gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  >  2:7696105/1‑93 (MQ=255)
                                                               gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  >  2:2180153/1‑93 (MQ=255)
                                                                accaccATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGttt                 <  2:737529/93‑1 (MQ=255)
                                                                  caccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGCATTCAGCGACTTGTCGCCCGTCAGTTTGc               >  2:7944763/1‑93 (MQ=255)
                                                                    ccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCg             >  2:5847036/1‑93 (MQ=255)
                                                                    ccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCg             >  1:6568300/1‑93 (MQ=255)
                                                                     caTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGt            >  2:259129/1‑93 (MQ=255)
                                                                      atatAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa           >  1:2512858/1‑93 (MQ=255)
                                                                       tataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAc          >  2:8575445/1‑93 (MQ=255)
                                                                        ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         <  1:3592304/93‑1 (MQ=255)
                                                                        ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         >  2:6843905/1‑93 (MQ=255)
                                                                         taAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        >  1:4063232/1‑93 (MQ=255)
                                                                         taAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        >  1:8809585/1‑93 (MQ=255)
                                                                          aaGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCAc       <  2:7008184/93‑1 (MQ=255)
                                                                          aaGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCAc       <  2:6601386/93‑1 (MQ=255)
                                                                           aGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACg      >  2:7549845/1‑93 (MQ=255)
                                                                             ccGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGcc    >  1:6932681/1‑93 (MQ=255)
                                                                             ccGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGcc    <  1:4922095/93‑1 (MQ=255)
                                                                              cGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCa   <  2:1392120/93‑1 (MQ=255)
                                                                              cGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCa   >  1:2699843/1‑93 (MQ=255)
                                                                              cGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCa   <  2:5299494/93‑1 (MQ=255)
                                                                               gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  <  2:4267596/93‑1 (MQ=255)
                                                                               gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  <  1:1434459/93‑1 (MQ=255)
                                                                               gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  >  2:5665255/1‑93 (MQ=255)
                                                                                      |                                                                                     
CATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCGGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013090‑4013260

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