Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,352,460 T→G 25.8% intergenic (+103/‑124) yjdJ → / → ghoS GNAT family putative N‑acetyltransferase/antitoxin of GhoTS toxin‑antitoxin pair; endonuclease for ghoT mRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,352,4600TG25.8% 260.2 / 121.4 162intergenic (+103/‑124)yjdJ/ghoSGNAT family putative N‑acetyltransferase/antitoxin of GhoTS toxin‑antitoxin pair; endonuclease for ghoT mRNA
Reads supporting (aligned to +/- strand):  ref base T (65/55);  new base G (25/17);  total (90/72)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-02

Only 100 of 174 total aligned reads displayed.

GAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCG  >  NC_000913/4352361‑4352559
                                                                                                   |                                                                                                   
 agagTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGtt                                                                                                    >  1:4117179/1‑100 (MQ=255)
   agTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGtttt                                                                                                  <  1:3522599/100‑1 (MQ=255)
     tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGtttttt                                                                                                <  2:1527671/100‑1 (MQ=255)
     tACAGAGTCACGATATTTTTCATTCGCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGTGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                                <  2:4058331/100‑1 (MQ=255)
      aCAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGttttttt                                                                                               <  2:3822964/100‑1 (MQ=255)
       cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTg                                                                                              <  2:267328/100‑1 (MQ=255)
       cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                              <  1:1562618/100‑1 (MQ=255)
       cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                              >  1:2578092/1‑100 (MQ=255)
        agagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGt                                                                                             >  1:1045262/1‑100 (MQ=255)
        agagTCACGATATTTTTCATTCGATCCGCGATGTGATTCAGGAGAGCATCTGAAGGGTAGGGGTATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTg                                                                                              <  2:1133912/99‑1 (MQ=255)
         gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTc                                                                                            <  1:815217/100‑1 (MQ=255)
           gTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTccc                                                                                          <  1:2142353/100‑1 (MQ=255)
             cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCaa                                                                                        >  1:1596142/1‑100 (MQ=255)
              aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAg                                                                                       <  1:2995060/100‑1 (MQ=255)
              aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAg                                                                                       <  1:2628969/100‑1 (MQ=255)
                gATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTc                                                                                     >  1:2726007/1‑100 (MQ=255)
                gATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTc                                                                                     <  1:519396/100‑1 (MQ=255)
                 atatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCa                                                                                    <  1:3763159/100‑1 (MQ=255)
                 atatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCa                                                                                    >  1:1050260/1‑100 (MQ=255)
                  tatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCAt                                                                                   >  1:3129102/1‑100 (MQ=255)
                   atTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATc                                                                                  <  2:1184344/100‑1 (MQ=255)
                     ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATccc                                                                                <  2:625855/100‑1 (MQ=255)
                     ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATccc                                                                                <  1:1609471/100‑1 (MQ=255)
                     ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                                >  1:483297/1‑100 (MQ=255)
                     ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                                <  1:4106442/100‑1 (MQ=255)
                      tttCATTCTCTCCGTGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATcccc                                                                               <  2:2817885/100‑1 (MQ=255)
                      tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGACTGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATcccc                                                                               <  1:2086821/100‑1 (MQ=255)
                      tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATcccc                                                                               <  2:698560/100‑1 (MQ=255)
                      tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATcccc                                                                               <  2:9588/100‑1 (MQ=255)
                      tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATcccc                                                                               <  2:3587988/100‑1 (MQ=255)
                         cATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCttt                                                                            >  1:2812586/1‑100 (MQ=255)
                          aTTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTa                                                                           >  1:96916/1‑100 (MQ=255)
                           ttCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTAc                                                                          <  1:3403962/100‑1 (MQ=255)
                             ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                        <  1:2495219/100‑1 (MQ=255)
                             ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                        <  1:385797/100‑1 (MQ=255)
                             ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                        <  2:3029613/100‑1 (MQ=255)
                               ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                      <  2:3719156/100‑1 (MQ=255)
                                tcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                     >  2:3661875/1‑100 (MQ=255)
                                 ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                    >  1:2201834/1‑100 (MQ=255)
                                  cgcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTTTCCCAAGTCATCCCCTTTACTGAGCaa                                                                   >  2:1556235/1‑100 (MQ=255)
                                   gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaa                                                                  <  1:3029809/100‑1 (MQ=255)
                                         tgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAg                                                            <  1:203515/100‑1 (MQ=255)
                                         tgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATACCCTTTACTGAGCAAAAAAAAg                                                            >  2:1906097/1‑100 (MQ=255)
                                          gATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGa                                                           <  1:3272871/100‑1 (MQ=255)
                                           aTGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGaa                                                          >  2:3675076/1‑100 (MQ=255)
                                            tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                         >  2:1885760/1‑100 (MQ=255)
                                            tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                         <  1:3207822/100‑1 (MQ=255)
                                             gCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAata                                                        <  2:566201/100‑1 (MQ=255)
                                              cAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAACAAAAGAata                                                        >  2:384769/1‑99 (MQ=255)
                                               agaagagCATCTGAAGGATAGGGGGAAGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCACGTCATCACCTTTACTGAGCAAAAAAACAAatat                                                       >  2:1094739/4‑99 (MQ=255)
                                                ggAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATAtct                                                     >  1:2604387/1‑100 (MQ=255)
                                                  agagCATCTGAAGGGTAAGGGGAAGCACAAAGAATGCGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCTTAAAAAAAAAAAAAAAAAGAATATCTcc                                                   <  2:2289916/100‑1 (MQ=255)
                                                     gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGTAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                                >  1:1894778/1‑100 (MQ=255)
                                                     gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                                >  2:85011/1‑100 (MQ=255)
                                                     gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                                >  1:3619125/1‑100 (MQ=255)
                                                     gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                                >  1:380305/1‑100 (MQ=255)
                                                      cATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtata                                               >  1:348812/1‑100 (MQ=255)
                                                         cTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTTTTTTTTTTTTTCCCATTAACTGAGCAAAAAAATTAATATCTCCttta                                               >  2:738237/1‑97 (MQ=255)
                                                         cTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATga                                            >  1:3754467/1‑100 (MQ=255)
                                                          tgtcgggTGGGGGGACACACAAAGAGGGGGAAGAGAGCGAGTTTTTTTGTCCTAAGTCATCCTCTTGATTGAAAAAAAAAAAGAATATCTCCTATATgag                                           <  2:2429750/96‑1 (MQ=255)
                                                          tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTTTTTTTTTTTTTCCCCTTTAATGAGCAAAAAAAAGAATATCTCCTATATg                                             >  2:1953249/1‑98 (MQ=255)
                                                          tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                           >  1:231379/1‑100 (MQ=255)
                                                            aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGaa                                         >  1:184433/1‑100 (MQ=255)
                                                             aGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCGCAAGTCATCCCCTGTAGTGAGCAAAAAAAAGAATAACTCCTATATGAGAat                                        >  2:3825130/1‑100 (MQ=255)
                                                              gggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatc                                       >  2:2617102/1‑100 (MQ=255)
                                                              gggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTTTTTTATTTCTTCCCCTTTACTGAGCAAAAAAAACAATATCTCCTATATga                                            >  2:2017023/1‑95 (MQ=255)
                                                                gTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcat                                     >  1:1678350/1‑100 (MQ=255)
                                                                 tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                    >  1:2526980/1‑100 (MQ=255)
                                                                 tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                    <  1:3888190/100‑1 (MQ=255)
                                                                    ggggATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAAt                                 >  1:3670399/1‑100 (MQ=255)
                                                                     gggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATc                                >  1:1800860/1‑100 (MQ=255)
                                                                       gATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgg                              <  1:149485/100‑1 (MQ=255)
                                                                        aTGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCggg                             >  2:1071367/1‑100 (MQ=255)
                                                                        aTGCACAAAGAATGCGCAGAGAGCGCGGTTTTTTGTCCAAAGTCATCACCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgg                              <  1:2026734/99‑1 (MQ=255)
                                                                          gCAAAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCATTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgggg                            <  2:53019/99‑1 (MQ=255)
                                                                           cacaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                          >  2:2244026/1‑100 (MQ=255)
                                                                           cacaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                          >  2:2827985/1‑100 (MQ=255)
                                                                           cacaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                          >  2:390340/1‑100 (MQ=255)
                                                                           cacaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                          >  1:3905945/1‑100 (MQ=255)
                                                                           cacaAAGAATGGGCAGAGAGCGCGGTTTTTTTTTTTTTTTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                          >  2:645849/1‑100 (MQ=255)
                                                                            acaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTta                         >  1:2847620/1‑100 (MQ=255)
                                                                             caAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaa                        >  1:2693032/1‑100 (MQ=255)
                                                                              aaaGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaat                       >  1:1830020/1‑100 (MQ=255)
                                                                               aaGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaata                      >  1:1139187/1‑100 (MQ=255)
                                                                                aGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaataa                     >  1:3985283/1‑100 (MQ=255)
                                                                                  aaTGTGCAGAGAGCGCGTTTTTTTGTCTCAAGTCATCTCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATTGGGGTTAATAAg                    <  2:4013188/99‑1 (MQ=255)
                                                                                   aTGGGCAGAGAGAGCGTTTTTTTGTCCCAAGTCATCCCCTTTAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGtt                  <  2:3784299/100‑1 (MQ=255)
                                                                                      ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg               <  1:1189224/100‑1 (MQ=255)
                                                                                      ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg               <  2:3268714/100‑1 (MQ=255)
                                                                                      ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg               <  2:2361597/100‑1 (MQ=255)
                                                                                      ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg               <  1:807674/100‑1 (MQ=255)
                                                                                      ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg               >  2:1468370/1‑100 (MQ=255)
                                                                                        cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg             >  1:2541034/1‑100 (MQ=255)
                                                                                         agagagCACGTGTGTTTGACTCATGTCAACTCCTTAAATGTTAAAAAAAATGAATAGCTCCTCTGTGAGAATCATCCATCTGTGCTAATAAGTTTTGCg             <  2:777348/99‑1 (MQ=255)
                                                                                           agagCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcc          <  2:2997329/100‑1 (MQ=255)
                                                                                             agCGCGTTTTTTTGTCCCAAGTCATCCCCTTAAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcccc        <  2:4012773/100‑1 (MQ=255)
                                                                                                gcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCaga     >  1:2203063/1‑100 (MQ=255)
                                                                                                 cgTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCagag    >  1:1354474/1‑100 (MQ=255)
                                                                                                  gTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGc   >  1:1965521/1‑100 (MQ=255)
                                                                                                   tttttttGTCCCAATTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGATTAATAAGTTTTGCGTCCCCAGAGCg  <  2:3029783/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAGAGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGAGCG  >  NC_000913/4352361‑4352559

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