Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,362,606 A→G D722G (GAC→GGC)  yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,6060AG99.6% 856.7 / ‑4.6 248D722G (GAC→GGC) yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (119/128);  total (120/128)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 258 total aligned reads displayed.

AAATGATCAAACAAGGCCGTGGTGCAGTAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTC  >  NC_000913/3362516‑3362694
                                                                                          |                                                                                        
aaaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCg                                                                                        <  1:5573842/93‑1 (MQ=255)
 aaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCcggcgg                                                                                       >  2:3412785/1‑93 (MQ=255)
 aaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCcggcgg                                                                                       >  1:5759094/1‑93 (MQ=255)
 aaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCcggcgg                                                                                       >  1:6037194/1‑93 (MQ=255)
  aTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggcggc                                                                                      >  2:244787/1‑93 (MQ=255)
    gATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCaa                                                                                    >  1:1573261/1‑93 (MQ=255)
     aTCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCaaa                                                                                   >  2:4978273/1‑93 (MQ=255)
      tCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCaaaa                                                                                  >  1:5275391/1‑93 (MQ=255)
         aaCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAAc                                                                               >  1:4578755/1‑93 (MQ=255)
         aaCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAAc                                                                               >  2:6344946/1‑93 (MQ=255)
          aCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACt                                                                              <  1:4124438/93‑1 (MQ=255)
             aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc                                                                           >  2:1214544/1‑93 (MQ=255)
             aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc                                                                           <  1:3281766/93‑1 (MQ=255)
             aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc                                                                           >  2:1270203/1‑93 (MQ=255)
             aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc                                                                           <  2:4694485/93‑1 (MQ=255)
              ggCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGccc                                                                          >  1:2579528/1‑93 (MQ=255)
              ggCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGccc                                                                          <  1:5672557/93‑1 (MQ=255)
              ggCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGccc                                                                          >  2:808571/1‑93 (MQ=255)
               gCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCa                                                                         <  2:1978718/93‑1 (MQ=255)
               gCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCa                                                                         >  1:1171831/1‑93 (MQ=255)
                  gtggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCagag                                                                      <  1:6176786/93‑1 (MQ=255)
                  gtggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCagag                                                                      <  1:520999/93‑1 (MQ=255)
                   tggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAgg                                                                     <  1:810659/93‑1 (MQ=255)
                    ggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGc                                                                    <  2:237313/93‑1 (MQ=255)
                      tgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGt                                                                  >  2:5961488/1‑93 (MQ=255)
                       gcagcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGtt                                                                 <  1:2145707/93‑1 (MQ=255)
                        cagcagGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTg                                                                >  2:2953683/1‑93 (MQ=255)
                        cagcagGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTg                                                                <  1:4366577/93‑1 (MQ=255)
                         agcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGc                                                               <  1:4832425/90‑1 (MQ=255)
                         agcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGc                                                               <  1:2779562/90‑1 (MQ=255)
                         agcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGc                                                               <  2:5663023/90‑1 (MQ=255)
                          gcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCg                                                              <  1:5460183/91‑1 (MQ=255)
                          gcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCg                                                              >  2:5986871/3‑93 (MQ=255)
                             gggAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGAtt                                                           >  1:4942662/1‑93 (MQ=255)
                              ggAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTg                                                          <  1:19882/93‑1 (MQ=255)
                              ggAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTg                                                          <  1:4812116/93‑1 (MQ=255)
                                aaaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGaa                                                        <  1:4978273/93‑1 (MQ=255)
                                aaaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGaa                                                        >  2:2192219/1‑93 (MQ=255)
                                 aaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAg                                                       >  2:3597466/1‑93 (MQ=255)
                                 aaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAg                                                       <  2:3984225/93‑1 (MQ=255)
                                 aaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAg                                                       <  2:2350215/93‑1 (MQ=255)
                                 aaGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAg                                                       >  1:1036709/1‑93 (MQ=255)
                                   gTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGAt                                                     <  2:195361/93‑1 (MQ=255)
                                    tCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATa                                                    <  2:6238716/93‑1 (MQ=255)
                                    tCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATa                                                    <  2:6439018/93‑1 (MQ=255)
                                       gATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTg                                                 >  1:6348250/1‑93 (MQ=255)
                                       gATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTg                                                 >  2:55377/1‑93 (MQ=255)
                                       gATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCACAGGCGTTGCGATTGAAGATAGTg                                                 >  1:3856283/1‑93 (MQ=255)
                                        aTTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGa                                                >  2:151953/1‑93 (MQ=255)
                                         ttCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGaa                                               >  1:2451944/1‑93 (MQ=255)
                                         ttCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGaa                                               <  2:1773102/93‑1 (MQ=255)
                                          tCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAg                                              <  1:4605528/93‑1 (MQ=255)
                                          tCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAg                                              <  2:266780/93‑1 (MQ=255)
                                           cAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAgg                                             <  1:5298717/93‑1 (MQ=255)
                                            aGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGc                                            <  2:3794851/93‑1 (MQ=255)
                                            aGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGc                                            <  2:3879266/93‑1 (MQ=255)
                                             gTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCa                                           >  1:2088074/1‑93 (MQ=255)
                                                 gATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTa                                       <  2:5838012/93‑1 (MQ=255)
                                                  aTAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTAt                                      >  1:650156/1‑93 (MQ=255)
                                                    aaCACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCt                                    <  1:3408460/93‑1 (MQ=255)
                                                    aaCACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCt                                    >  2:534152/1‑93 (MQ=255)
                                                      cacaACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGa                                  <  2:5441659/93‑1 (MQ=255)
                                                       acaACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGac                                 >  1:4084485/1‑93 (MQ=255)
                                                       acaACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGac                                 <  2:3099679/93‑1 (MQ=255)
                                                        caACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGaca                                >  1:1519916/1‑93 (MQ=255)
                                                        caACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGaca                                <  1:3528147/93‑1 (MQ=255)
                                                         aaCGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGacaa                               <  2:4857326/93‑1 (MQ=255)
                                                         aaCGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGacaa                               >  2:4488958/1‑93 (MQ=255)
                                                         aaCGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGacaa                               >  2:2125149/1‑93 (MQ=255)
                                                          aCGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGacaac                              >  2:4146490/1‑93 (MQ=255)
                                                          aCGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGAGAGTGAAGGCAACTATCTGacaac                              >  1:1499550/1‑93 (MQ=255)
                                                             ccGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATc                           >  2:4943701/1‑93 (MQ=255)
                                                             ccGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATc                           <  1:3216658/93‑1 (MQ=255)
                                                               gTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAg                         <  1:4638031/93‑1 (MQ=255)
                                                               gTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAg                         <  1:1948581/93‑1 (MQ=255)
                                                               gTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAg                         >  2:5203679/1‑93 (MQ=255)
                                                               gTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAg                         >  2:2048600/1‑93 (MQ=255)
                                                               gTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAg                         <  2:1019924/93‑1 (MQ=255)
                                                                tGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGt                        <  1:37728/93‑1 (MQ=255)
                                                                tGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGt                        <  1:2027136/93‑1 (MQ=255)
                                                                tGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGt                        <  1:4261400/93‑1 (MQ=255)
                                                                 gCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTg                       <  2:5728514/93‑1 (MQ=255)
                                                                     tACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATg                   <  1:3009519/93‑1 (MQ=255)
                                                                      aCTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGa                  <  1:1696280/93‑1 (MQ=255)
                                                                       cTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGAc                 <  1:5849325/93‑1 (MQ=255)
                                                                       cTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGAc                 <  1:4319713/93‑1 (MQ=255)
                                                                        tCAATATCACACTTGCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACg                <  1:3918124/93‑1 (MQ=255)
                                                                          aaTATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgt              <  2:2495328/93‑1 (MQ=255)
                                                                          aaTATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgt              <  1:24764/93‑1 (MQ=255)
                                                                          aaTATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgt              >  2:6041606/1‑93 (MQ=255)
                                                                           atatCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgtg             <  2:2829172/93‑1 (MQ=255)
                                                                            tatCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgtgt            >  2:3241402/1‑93 (MQ=255)
                                                                            tatCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGgtgt            >  1:2846424/1‑93 (MQ=255)
                                                                             atCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGtgtt           <  2:1755689/93‑1 (MQ=255)
                                                                             atCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGtgtt           <  2:1497911/93‑1 (MQ=255)
                                                                              tCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGtgttg          >  1:3981900/1‑93 (MQ=255)
                                                                                  acTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGGGTTGTAtt      >  2:941574/1‑93 (MQ=255)
                                                                                   cTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTc     <  2:1693006/93‑1 (MQ=255)
                                                                                     ttCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCt   <  2:3309823/93‑1 (MQ=255)
                                                                                     ttCCGGCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCt   <  1:5933822/93‑1 (MQ=255)
                                                                                          |                                                                                        
AAATGATCAAACAAGGCCGTGGTGCAGTAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTC  >  NC_000913/3362516‑3362694

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