Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,689,553 A→G intergenic (+169/+57) ydgA → / ← uidC hypothetical protein/membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,5530AG99.2% 402.5 / ‑5.4 128intergenic (+169/+57)ydgA/uidChypothetical protein/membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (51/76);  total (51/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 134 total aligned reads displayed.

GCTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAT  >  NC_000913/1689463‑1689642
                                                                                          |                                                                                         
 cTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAgcgc                                                                                        <  1:68676/93‑1 (MQ=255)
  ttATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGcgcg                                                                                       >  1:2438926/1‑93 (MQ=255)
  ttATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGcgcg                                                                                       >  1:2443852/1‑93 (MQ=255)
   tATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGGCGCATCCGGCATAAAAAgcgcgc                                                                                      >  1:250441/1‑93 (MQ=255)
      cAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACt                                                                                   <  2:1248133/93‑1 (MQ=255)
      cAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACt                                                                                   <  1:3437709/93‑1 (MQ=255)
       aGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACtt                                                                                  <  2:417402/93‑1 (MQ=255)
       aGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTACGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACtt                                                                                  <  1:4435135/93‑1 (MQ=255)
         gCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTg                                                                                <  1:3107975/93‑1 (MQ=255)
         gCCTACCAATACCTGCAACACATTGAATTTATGAGTATTTGTAGGCAGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTg                                                                                <  2:4343774/93‑1 (MQ=18)
           cTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTc                                                                              >  1:4562059/1‑93 (MQ=255)
           cTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTc                                                                              >  1:655081/1‑93 (MQ=255)
           cTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATACGGCATAAAAGGCGCGCACTTTGTc                                                                              >  1:1985359/1‑93 (MQ=255)
             aCCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaa                                                                            >  2:4622630/1‑93 (MQ=255)
             aCCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaa                                                                            >  2:4622631/1‑93 (MQ=255)
              ccTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaac                                                                           <  1:1078130/93‑1 (MQ=255)
               cTTTTCCTGCAACACTTTGAATTTATGAGGTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaaca                                                                          >  1:1899748/1‑93 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         <  1:956096/93‑1 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         >  2:3848539/1‑93 (MQ=255)
                 tttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAAt                                                                        <  2:2131093/93‑1 (MQ=255)
                  ttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATc                                                                       <  1:2226944/93‑1 (MQ=255)
                  ttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATc                                                                       <  1:2257973/93‑1 (MQ=255)
                   tCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCt                                                                      <  2:275769/93‑1 (MQ=255)
                   tCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCt                                                                      <  1:4387610/93‑1 (MQ=255)
                    ccTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTg                                                                     <  2:2959769/93‑1 (MQ=255)
                    ccTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTg                                                                     <  1:387893/93‑1 (MQ=255)
                     cTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGa                                                                    >  2:1746066/1‑93 (MQ=255)
                     cTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGa                                                                    <  2:784727/93‑1 (MQ=255)
                      tGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaa                                                                   <  2:1233878/93‑1 (MQ=255)
                      tGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTAAGACGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaa                                                                   >  2:821380/1‑93 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  1:91233/93‑1 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  >  1:2419612/1‑93 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  1:303178/93‑1 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  2:3091097/93‑1 (MQ=255)
                        cAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAAc                                                                 >  1:3183493/1‑93 (MQ=255)
                        cAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAAc                                                                 <  1:1172010/93‑1 (MQ=255)
                          acacTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGc                                                               >  1:1053419/1‑93 (MQ=255)
                          acacTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGc                                                               >  1:4479928/1‑93 (MQ=255)
                           cacTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGcc                                                              <  2:4136129/93‑1 (MQ=255)
                            acTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCg                                                             >  2:3494334/1‑93 (MQ=255)
                            acTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCg                                                             >  1:125018/1‑93 (MQ=255)
                            acTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCg                                                             <  2:246817/93‑1 (MQ=255)
                             cTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCgg                                                            <  2:1639646/93‑1 (MQ=255)
                              tttGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGga                                                           <  2:4638721/93‑1 (MQ=255)
                               ttGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGgag                                                          >  2:726754/1‑93 (MQ=255)
                                 gAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAt                                                        <  1:1387181/93‑1 (MQ=255)
                                  aaTTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGGCGGAGAtt                                                       >  1:2141741/1‑93 (MQ=255)
                                   aTTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAttt                                                      <  2:1334309/93‑1 (MQ=255)
                                      tATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtc                                                   <  1:752532/93‑1 (MQ=255)
                                      tATGAGTTTTTGTAGGCTGGATAAGGCGTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtct                                                  <  2:740409/93‑1 (MQ=255)
                                       aTGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtct                                                  >  2:1891430/1‑93 (MQ=255)
                                          aGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtctctc                                               >  1:1074281/1‑93 (MQ=255)
                                            tttttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCg                                             <  2:1947615/93‑1 (MQ=255)
                                                tGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGt                                         >  2:3966110/1‑93 (MQ=255)
                                                 gTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTCGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGtt                                        <  2:2065671/93‑1 (MQ=255)
                                                  tAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTa                                       <  1:1332465/93‑1 (MQ=255)
                                                  tAGGCTGGATAAGGCGTTGACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTa                                       <  2:3615660/93‑1 (MQ=255)
                                                   aGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAt                                      <  1:786565/93‑1 (MQ=255)
                                                    ggCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtt                                     >  1:953722/1‑93 (MQ=255)
                                                     gCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAttt                                    >  2:1966251/1‑93 (MQ=255)
                                                     gCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAttt                                    >  1:1537309/1‑93 (MQ=255)
                                                     gCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAttt                                    <  2:4531295/93‑1 (MQ=255)
                                                     gCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATCTTTCTCTCCGGCGTTAttt                                    >  1:4304865/1‑93 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  1:4556535/93‑1 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  2:2162452/93‑1 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  1:285331/93‑1 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  1:285332/93‑1 (MQ=255)
                                                          aTAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTAc                               >  2:2856114/1‑93 (MQ=255)
                                                            aaGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACtt                             >  2:3189886/1‑93 (MQ=255)
                                                             aGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTc                            >  1:612253/1‑93 (MQ=255)
                                                             aGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTc                            <  2:3326017/93‑1 (MQ=255)
                                                              ggCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGATATTTTTCTCTCCGGCGTTATTTTTTACTTCa                           <  1:2474711/93‑1 (MQ=255)
                                                                 gTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTATTTCAGct                        >  1:3389941/1‑92 (MQ=255)
                                                                 gTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCa                        <  1:1894043/93‑1 (MQ=255)
                                                                  tttACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCAt                       <  2:3469327/93‑1 (MQ=255)
                                                                     aCGCGGCATCCGGCATAAAAAGCGCGCACTGTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATaaa                    <  1:233524/93‑1 (MQ=255)
                                                                     aCGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATaaa                    <  2:1670033/93‑1 (MQ=255)
                                                                       gccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGt                  <  2:2821633/93‑1 (MQ=255)
                                                                       gccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGt                  <  1:888083/93‑1 (MQ=255)
                                                                       cccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGt                  <  2:786795/92‑1 (MQ=255)
                                                                        ccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTc                 <  1:632654/93‑1 (MQ=255)
                                                                         cgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCa                <  1:1759912/93‑1 (MQ=255)
                                                                         cgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCa                <  1:4064197/93‑1 (MQ=255)
                                                                         cgcATCCGGCATAAAAAGCGCGCACTTGGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCa                <  1:4064339/93‑1 (MQ=255)
                                                                             tCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACt            <  2:3192150/93‑1 (MQ=255)
                                                                             tCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACt            >  2:2250259/1‑93 (MQ=255)
                                                                             tCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACt            <  1:3557407/93‑1 (MQ=255)
                                                                              ccGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTg           <  1:1165854/93‑1 (MQ=255)
                                                                               cGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGa          <  2:1022886/93‑1 (MQ=255)
                                                                               cGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGa          >  1:2837706/1‑93 (MQ=255)
                                                                                ggCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGaa         <  1:4314445/93‑1 (MQ=255)
                                                                                 gCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAAt        <  1:1400794/93‑1 (MQ=255)
                                                                                 gCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAAt        <  1:1400793/93‑1 (MQ=255)
                                                                                  cATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATa       <  2:3728286/93‑1 (MQ=255)
                                                                                   aTAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATAc      >  1:2632865/1‑93 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  1:117541/1‑93 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  2:4281164/1‑93 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  1:1661849/1‑93 (MQ=255)
                                                                                    tAAAAAACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     <  1:4260329/93‑1 (MQ=255)
                                                                                       aaagCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAt  <  1:4439799/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAT  >  NC_000913/1689463‑1689642

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