Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,761,1480GA59.2% 344.9 / inf 412T34T (ACC→ACT) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (82/86);  new base A (97/147);  total (179/233)
Fisher's exact test for biased strand distribution p-value = 6.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 415 total aligned reads displayed.

TCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTC  >  NC_000913/3761050‑3761247
                                                                                                  |                                                                                                   
tcttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAgt                                                                                                   <  1:7126642/101‑1 (MQ=255)
 cttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGgtg                                                                                                  >  1:7127207/1‑101 (MQ=255)
 cttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGgtg                                                                                                  >  2:7127206/1‑101 (MQ=255)
 cttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGgtg                                                                                                  <  2:7126479/101‑1 (MQ=255)
 cttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAgtg                                                                                                  <  1:7126051/101‑1 (MQ=255)
   tGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGtgtg                                                                                                <  2:7126129/101‑1 (MQ=255)
    gcgcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGa                                                                                               <  2:7126840/101‑1 (MQ=255)
    gcgcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGa                                                                                               <  2:7126581/101‑1 (MQ=255)
      gcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGAcc                                                                                             >  1:7127215/1‑101 (MQ=255)
        ccAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCAt                                                                                           >  2:7127217/1‑101 (MQ=255)
         cAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATa                                                                                          <  1:7126572/101‑1 (MQ=255)
         cAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATa                                                                                          <  2:7126498/101‑1 (MQ=255)
          aGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATaa                                                                                         <  2:7126821/101‑1 (MQ=255)
          aGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATaa                                                                                         >  1:7127221/1‑101 (MQ=255)
            ttttCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGt                                                                                       <  1:7126484/101‑1 (MQ=255)
                 acacCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCAc                                                                                  >  1:7127226/1‑101 (MQ=255)
                    ccacGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACg                                                                                 <  2:7126862/99‑1 (MQ=255)
                    ccacGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCa                                                                               <  2:7126234/101‑1 (MQ=255)
                      acGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGc                                                                                >  2:7127232/1‑98 (MQ=255)
                       cGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCaaag                                                                            <  2:7126633/101‑1 (MQ=255)
                        gcaggcagCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAg                                                                                        <  1:7127100/88‑1 (MQ=255)
                         caggcagCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCaaagaa                                                                          >  1:7127238/1‑101 (MQ=255)
                         caggcagCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCaaagaa                                                                          <  2:7126792/101‑1 (MQ=255)
                            gcagCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAaagg                                                                       >  2:7127242/1‑101 (MQ=255)
                              agCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAaaggaa                                                                     >  2:7127247/1‑101 (MQ=255)
                               gCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAaaggaag                                                                    >  2:7127251/1‑101 (MQ=255)
                                 gTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAaggaagga                                                                  >  1:7127254/1‑101 (MQ=255)
                                  tCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAg                                                                 >  1:7127259/1‑101 (MQ=255)
                                   cTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGc                                                                >  2:7127261/1‑101 (MQ=255)
                                    tGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGt                                                                                       <  2:7126342/77‑1 (MQ=255)
                                    tGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCt                                                               <  2:7126411/101‑1 (MQ=255)
                                    tGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCt                                                               <  2:7126526/101‑1 (MQ=255)
                                     ggCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATa                                                                                          <  1:7126336/73‑1 (MQ=255)
                                       cTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATc                                                            >  1:7127268/1‑101 (MQ=255)
                                       cTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCaaagaaag                                                                        <  2:7126266/89‑1 (MQ=255)
                                        tGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcg                                                           >  2:7127270/1‑101 (MQ=255)
                                          ccACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcgcg                                                         <  2:7126954/101‑1 (MQ=255)
                                          ccACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcgcg                                                         <  2:7127052/101‑1 (MQ=255)
                                          ccACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcgcg                                                         <  2:7126710/101‑1 (MQ=255)
                                           cACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCgcgc                                                        <  2:7126906/101‑1 (MQ=255)
                                           cACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCgcgc                                                        >  1:7127274/1‑101 (MQ=255)
                                            aCGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCa                                                       <  2:7126571/101‑1 (MQ=255)
                                             cGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcgc                                                          <  1:7126716/97‑1 (MQ=255)
                                             cGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCaa                                                      >  2:7127277/1‑101 (MQ=255)
                                              gAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAAt                                                     <  1:7126515/101‑1 (MQ=255)
                                               aaTCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATaa                                                                                         <  1:7126619/64‑1 (MQ=255)
                                                 tCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAg                                                  <  1:7126491/101‑1 (MQ=255)
                                                 tCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAg                                                  >  1:7127282/1‑101 (MQ=255)
                                                  cACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGc                                                 <  1:7126738/101‑1 (MQ=255)
                                                   aCTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAgt                                                                                                   <  2:7126539/50‑1 (MQ=255)
                                                   aCTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAATGAAAGGAAGGAGCTATCGCGCAATAAGc                                                 <  1:7126771/100‑1 (MQ=255)
                                                   aCTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCg                                                >  2:7127285/1‑101 (MQ=255)
                                                     ttCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGAt                                              >  1:7127293/1‑101 (MQ=255)
                                                      tctcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATc                                             >  2:7127294/1‑101 (MQ=255)
                                                        tcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCaa                                           <  2:7126443/101‑1 (MQ=255)
                                                         cgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAAt                                          <  2:7126656/101‑1 (MQ=255)
                                                         cgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAAt                                          <  2:7126617/101‑1 (MQ=255)
                                                         cgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAAt                                          <  2:7126675/101‑1 (MQ=255)
                                                           cgcTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAg                                        >  1:7127302/1‑101 (MQ=255)
                                                            gcTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGc                                       >  2:7127303/1‑101 (MQ=255)
                                                             cTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGcc                                      <  1:7126552/101‑1 (MQ=255)
                                                             cTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGcc                                      <  1:7126838/101‑1 (MQ=255)
                                                              ttCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCg                                     >  2:7127312/1‑101 (MQ=255)
                                                              ttCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCg                                     >  1:7127310/1‑101 (MQ=255)
                                                              ttCGTCGAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCg                                     <  1:7126605/101‑1 (MQ=255)
                                                               tcgtcgAGCATCTGACGCAACAATTACACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCgg                                    <  1:7126174/101‑1 (MQ=255)
                                                                cgtcgAGCATCTGAGGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAg                                                                 >  2:7127321/1‑71 (MQ=255)
                                                                cgtcgAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGa                                   >  1:7127316/1‑101 (MQ=255)
                                                                  tcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAg                                 <  2:7126598/101‑1 (MQ=255)
                                                                   cgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGt                                >  2:7127326/1‑101 (MQ=255)
                                                                   cgAGCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGt                                <  1:7126665/101‑1 (MQ=255)
                                                                    gAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGtt                               >  2:7127330/1‑101 (MQ=255)
                                                                      gCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGa                             >  2:7127334/1‑101 (MQ=255)
                                                                      gCATCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGa                             <  1:7126966/101‑1 (MQ=255)
                                                                       cATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGAc                            >  1:7127336/1‑101 (MQ=255)
                                                                        aTCTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACt                           <  1:7126782/101‑1 (MQ=255)
                                                                          cTGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACtat                         <  2:7126739/101‑1 (MQ=255)
                                                                           tGACGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACtata                        >  1:7127341/1‑101 (MQ=255)
                                                                            gACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATag                       <  2:7126649/101‑1 (MQ=255)
                                                                              cGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATagag                     <  2:7126591/101‑1 (MQ=255)
                                                                              cGCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATagag                     <  2:7126807/101‑1 (MQ=255)
                                                                               gCAACAATTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAgg                    >  2:7127347/1‑101 (MQ=255)
                                                                                caacaaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAggg                   >  1:7127350/1‑101 (MQ=255)
                                                                                 aacaaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGa                  <  2:7127024/101‑1 (MQ=255)
                                                                                 aacaaTTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGa                  >  2:7127352/1‑101 (MQ=255)
                                                                                   caaTTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAAt                                          >  2:7127361/1‑75 (MQ=255)
                                                                                   caaTTCCACCACGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAAc                >  1:7127360/1‑101 (MQ=255)
                                                                                    aaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAAcg               <  1:7126401/101‑1 (MQ=255)
                                                                                       tCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATcgcg                                                         >  1:7127370/1‑56 (MQ=255)
                                                                                       tCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGt            >  2:7127366/1‑101 (MQ=255)
                                                                                         caccacGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGttt          <  2:7126697/101‑1 (MQ=255)
                                                                                         caccacGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGttt          <  2:7126753/101‑1 (MQ=255)
                                                                                          accacGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTc         <  1:7126879/101‑1 (MQ=255)
                                                                                          accacGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTc         >  2:7127376/1‑101 (MQ=255)
                                                                                           ccacGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCg        <  2:7126559/101‑1 (MQ=255)
                                                                                           ccacGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCg        <  1:7126451/101‑1 (MQ=255)
                                                                                           ccacGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCg        >  1:7127380/1‑101 (MQ=255)
                                                                                            cacGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCgg       >  2:7127384/1‑101 (MQ=255)
                                                                                            cacGCGAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCgg       <  1:7126810/101‑1 (MQ=255)
                                                                                               gcgAGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTg    <  2:7126874/101‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTC  >  NC_000913/3761050‑3761247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑