Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,372,245 A→G S44S (TCA→TCG ycjP → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,372,2450AG87.8% 689.1 / 77.2 272S44S (TCA→TCGycjPputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (20/13);  new base G (108/130);  total (128/144)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 272 total aligned reads displayed.

CTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCC  >  NC_000913/1372155‑1372336
                                                                                          |                                                                                           
cTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCt                                                                                           <  2:7184352/93‑1 (MQ=255)
 tGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTg                                                                                          >  2:5240803/1‑93 (MQ=255)
 tGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTg                                                                                          >  1:2106990/1‑93 (MQ=255)
  gCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGc                                                                                         <  1:2429338/93‑1 (MQ=255)
  gCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGc                                                                                         >  2:2412663/1‑93 (MQ=255)
   cGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCa                                                                                        <  2:2588830/93‑1 (MQ=255)
    gggCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCAt                                                                                       <  1:6982751/93‑1 (MQ=255)
     ggCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATc                                                                                      <  1:1716137/93‑1 (MQ=255)
      gCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATcc                                                                                     <  2:373647/93‑1 (MQ=255)
       cTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCt                                                                                    <  1:423424/93‑1 (MQ=255)
        tGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTa                                                                                   >  1:7347488/1‑93 (MQ=255)
           cgcTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgc                                                                                >  1:4146292/1‑93 (MQ=255)
            gcTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgct                                                                               <  2:6704328/93‑1 (MQ=255)
              tGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgctgc                                                                             >  2:2012409/1‑93 (MQ=255)
              tGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgctgc                                                                             >  1:5742887/1‑93 (MQ=255)
              tGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgctgc                                                                             >  1:5460641/1‑93 (MQ=255)
              tGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgctgc                                                                             <  1:7117687/93‑1 (MQ=255)
              tGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACgctgc                                                                             <  2:6127850/93‑1 (MQ=255)
                tttCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGctgctg                                                                           <  1:550317/93‑1 (MQ=255)
                 ttCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCtgctgc                                                                          >  2:459744/1‑93 (MQ=255)
                 ttCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCtgctgc                                                                          <  2:1236787/93‑1 (MQ=255)
                    tcatcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgc                                                                       >  1:3420080/1‑93 (MQ=255)
                     catcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgca                                                                      <  2:2181996/93‑1 (MQ=255)
                     catcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgca                                                                      >  2:5599087/1‑93 (MQ=255)
                      atcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgcag                                                                     <  1:2065054/93‑1 (MQ=255)
                       tcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgcagc                                                                    >  1:2328028/1‑93 (MQ=255)
                       tcatcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgcagc                                                                    <  2:1660043/93‑1 (MQ=255)
                        catcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgcagca                                                                   <  2:3651337/93‑1 (MQ=255)
                        catcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCgcagca                                                                   <  1:6663029/93‑1 (MQ=255)
                         atcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGcagcag                                                                  <  1:6501902/93‑1 (MQ=255)
                         atcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGcagcag                                                                  <  2:3112602/93‑1 (MQ=255)
                          tcaCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGt                                                                 >  2:4081423/1‑93 (MQ=255)
                            aCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTgg                                                               <  2:5439567/93‑1 (MQ=255)
                             cGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGa                                                              <  1:995089/93‑1 (MQ=255)
                             cGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGa                                                              <  2:357070/93‑1 (MQ=255)
                                tGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGc                                                           <  2:1276638/93‑1 (MQ=255)
                                tGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGc                                                           >  2:862320/1‑93 (MQ=255)
                                tGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGc                                                           >  1:2567123/1‑93 (MQ=255)
                                   ttCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTgg                                                        <  1:4918072/93‑1 (MQ=255)
                                    tCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGa                                                       <  1:146724/93‑1 (MQ=255)
                                     ccATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAg                                                      >  2:2807697/1‑93 (MQ=255)
                                     ccATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAg                                                      <  2:6820541/93‑1 (MQ=255)
                                       aTTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCa                                                    <  1:4840079/93‑1 (MQ=255)
                                       aTTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCa                                                    >  1:6014254/1‑93 (MQ=255)
                                        ttttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGTAGTGGACGCTGGAGCAt                                                   <  2:4904373/93‑1 (MQ=255)
                                        ttttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAt                                                   <  2:3765775/93‑1 (MQ=255)
                                        ttttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAt                                                   <  2:5130783/93‑1 (MQ=255)
                                        ttttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAt                                                   <  2:1015849/93‑1 (MQ=255)
                                        ttttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAt                                                   <  2:442919/93‑1 (MQ=255)
                                         tttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAtt                                                  <  2:572349/93‑1 (MQ=255)
                                          ttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTa                                                 <  2:3536257/93‑1 (MQ=255)
                                          ttttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTa                                                 <  1:7477203/93‑1 (MQ=255)
                                           tttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTAc                                                <  2:3862408/93‑1 (MQ=255)
                                            ttGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACg                                               >  1:4014230/1‑93 (MQ=255)
                                             tgtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGt                                              <  1:1182844/93‑1 (MQ=255)
                                             tgtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGt                                              >  2:7245201/1‑93 (MQ=255)
                                              gtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTc                                             <  1:2210260/93‑1 (MQ=255)
                                              gtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTc                                             >  2:4675773/1‑93 (MQ=255)
                                              gtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTc                                             >  2:6205050/1‑93 (MQ=255)
                                              gtgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTc                                             >  2:402648/1‑93 (MQ=255)
                                               tgATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCg                                            <  1:5058518/93‑1 (MQ=255)
                                                   gCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAt                                        >  1:6286625/1‑93 (MQ=255)
                                                    cTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAtt                                       >  1:1541107/1‑93 (MQ=255)
                                                     tgatgaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAttt                                      <  2:4230568/93‑1 (MQ=255)
                                                     tgatgaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAttt                                      >  1:3055345/1‑93 (MQ=255)
                                                     tgatgaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAttt                                      <  2:4161408/93‑1 (MQ=255)
                                                      gatgaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAtttt                                     >  2:398565/1‑93 (MQ=255)
                                                      gatgaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAtttt                                     >  1:1097156/1‑93 (MQ=255)
                                                         gaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTaa                                  <  1:1323041/93‑1 (MQ=255)
                                                         gaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTaa                                  >  2:6487707/1‑93 (MQ=255)
                                                         gaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTaa                                  >  1:1582400/1‑93 (MQ=255)
                                                         gaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTaa                                  <  1:682935/93‑1 (MQ=255)
                                                         gaCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTaa                                  >  1:2956190/1‑93 (MQ=255)
                                                              cGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCga                             <  1:3624908/93‑1 (MQ=255)
                                                              cGTTCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCga                             <  2:5756923/93‑1 (MQ=255)
                                                              cGTTCAAGGGCGCGAAAGAGGCGATCTCCCTGCAACCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCga                             <  1:896691/93‑1 (MQ=255)
                                                                ttCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatg                           >  2:3422539/1‑93 (MQ=255)
                                                                ttCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatg                           >  1:2637947/1‑93 (MQ=255)
                                                                ttCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatg                           >  1:7232360/1‑93 (MQ=255)
                                                                 tCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatga                          <  1:6896107/93‑1 (MQ=255)
                                                                 tCAAGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatga                          >  2:1793302/1‑93 (MQ=255)
                                                                 tCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCgatga                          >  2:7381091/1‑93 (MQ=255)
                                                                   aaGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAtt                        >  1:4752422/1‑93 (MQ=255)
                                                                   aaGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAtt                        >  1:6082732/1‑93 (MQ=255)
                                                                    aGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAttt                       <  2:1478127/93‑1 (MQ=255)
                                                                    aGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAttt                       >  2:7426019/1‑93 (MQ=255)
                                                                    aGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAttt                       >  2:2822578/1‑93 (MQ=255)
                                                                    aGGGCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGAttt                       >  2:3956078/1‑93 (MQ=255)
                                                                      ggCGCGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTc                     <  1:5200288/93‑1 (MQ=255)
                                                                       gcgcGAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTcc                    <  2:5837039/93‑1 (MQ=255)
                                                                         gcgAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCg                  <  2:4713218/93‑1 (MQ=255)
                                                                         gcgAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCg                  >  2:802902/1‑93 (MQ=255)
                                                                           gAAAGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGtt                >  1:5943506/1‑93 (MQ=255)
                                                                             aaGAGGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTg              <  2:6823684/93‑1 (MQ=255)
                                                                             aaGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTg              >  1:6436495/1‑93 (MQ=255)
                                                                              agagGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGt             >  1:3271764/1‑93 (MQ=255)
                                                                                agGCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCg           >  2:2363413/1‑93 (MQ=255)
                                                                                 ggCGATCTCGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGa          <  1:6605203/93‑1 (MQ=255)
                                                                                     aTCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGactac      >  1:4350142/1‑93 (MQ=255)
                                                                                      tctcGCTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGactact     >  1:849892/1‑93 (MQ=255)
                                                                                          |                                                                                           
CTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCC  >  NC_000913/1372155‑1372336

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