Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,352,460 T→G 100% 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,4600TG95.1% 347.8 / 3.7 106intergenic (+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 (3/3);  new base G (58/42);  total (61/45)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

AGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGA  >  NC_000913/4352364‑4352556
                                                                                                |                                                                                                
agTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttt                                                                                               <  2:2292779/100‑1 (MQ=255)
 gTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttt                                                                                              <  2:2164452/100‑1 (MQ=255)
 gTACAGAGTCACGATATTTTTCATTATCTCAGCGATGTGATGAAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGGAGAGAGCGCGgtt                                                                                                >  2:1454337/1‑98 (MQ=255)
  tACAGTGTCATGATATTTTTCATTCTCTCCGCGATGTTATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             <  1:2331788/100‑1 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCGGAGAGCGCGgttttt                                                                                             >  1:2471323/1‑100 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             <  1:1722239/100‑1 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             <  2:1170765/100‑1 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             <  1:1608387/100‑1 (MQ=255)
   aCAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttttt                                                                                            <  1:2213758/100‑1 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATTTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           <  2:590960/100‑1 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           <  1:1032201/100‑1 (MQ=255)
      gagTTACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           <  1:2399721/98‑1 (MQ=255)
       agTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTcc                                                                                        <  2:705444/100‑1 (MQ=255)
       agTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTcc                                                                                        >  2:699991/1‑100 (MQ=255)
          cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCaa                                                                                     <  1:2315400/100‑1 (MQ=255)
          cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCaa                                                                                     >  1:1250534/1‑100 (MQ=255)
           aCGATGTTTTTCATTCTTTCCGCGATGTGATCCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAg                                                                                    <  2:810398/100‑1 (MQ=255)
           aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAg                                                                                    <  1:1526383/100‑1 (MQ=255)
           aCGATATTTTTAATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGTGGTTTTTTGTCCCAAg                                                                                    <  2:382791/100‑1 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCTCa                                                                                      >  2:2524864/1‑97 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGt                                                                                   <  2:2267077/100‑1 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGt                                                                                   >  1:1760640/1‑100 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGATTTTTGTCCCAAGt                                                                                   >  1:2896266/1‑100 (MQ=255)
             gATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGAGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTGCCAAGt                                                                                   >  1:2604698/1‑99 (MQ=255)
               tatTTTTCATTCTCTCCGCGATGTGAGTCAGGAGAGCATCTGAAGGGTAGGGGGATGTACAAAGAATGGTCATAGAGCGTTTTTTTTTGTCCGAAGTc                                                                                  <  2:2878609/98‑1 (MQ=255)
                atTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTGTGTCCCAAGTCATc                                                                               >  2:2066481/1‑100 (MQ=255)
                  ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                             >  1:27195/1‑100 (MQ=255)
                    ttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCt                                                                           >  1:813721/1‑100 (MQ=255)
                     tCATTCTCTCCGAGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCATAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCtt                                                                          <  2:2276892/100‑1 (MQ=255)
                        ttCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAc                                                                       >  1:392752/1‑100 (MQ=255)
                        ttCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAc                                                                       <  2:861568/100‑1 (MQ=255)
                         tctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACt                                                                      >  2:851246/1‑100 (MQ=255)
                          ctctcCGCGATGTGATGCAGGAGAGTATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                     <  1:1635070/100‑1 (MQ=255)
                          ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                     <  1:2485774/100‑1 (MQ=255)
                           tttgagcgATGTGATTCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                     <  2:1772639/94‑1 (MQ=255)
                            ctcCGCTATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                   <  1:2771522/100‑1 (MQ=255)
                            ctcCGCGATGTGATGCAGGAGAGCATCTGAAGTGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                   <  2:2507119/100‑1 (MQ=255)
                            ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGa                                                                    >  2:853791/1‑99 (MQ=255)
                            ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                   >  1:2172089/1‑100 (MQ=255)
                              ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                 >  2:2587733/1‑100 (MQ=255)
                               cgcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaa                                                                >  2:2259940/1‑100 (MQ=255)
                                 cgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaa                                                              <  1:1615820/100‑1 (MQ=255)
                                    tgtgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTTTCCCAAGTCATCCCCTTTACTGAGCaaaaaaa                                                           >  2:1109476/1‑100 (MQ=255)
                                     gtgATGCAAGAGAGCATATGAAGTGTTGGTTTATGTATAAAGAATGGGCAGAGAGCGCTGTTTTTTGTCCCAAGTTATCCCTTATAAAGAGaaaaaa                                                             <  2:2468955/97‑1 (MQ=255)
                                       gATGCAGGTGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGa                                                        >  1:2199473/1‑100 (MQ=255)
                                       gATGCAGGTGAGCATCTGAAGGGTACGGGGATGCACAAATAATGTGCAGAGAGCGCGTTTTTTTGTTCCAAGTTTTTTTCTTTACTGAGCAAAAAAAAGa                                                        <  2:570028/100‑1 (MQ=255)
                                       gATGCAGCAGAGCATCTGAAAGGAAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTACCAAGTCATCCGCTTTACTGAGCAAAAAAAAGa                                                        >  2:730828/1‑100 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGaa                                                       >  1:2230688/1‑99 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                      >  2:914227/1‑100 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                      <  1:914233/100‑1 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGAACAAAGACTGGGTCGAGAGAGCGGTTTTTTGTCCCAAGTCCTCCCCTTTCCTGAGCAAAAAAAAGAat                                                      <  1:1356576/100‑1 (MQ=255)
                                          gCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAata                                                     >  1:1546577/1‑100 (MQ=255)
                                           cAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAatat                                                    >  1:2237431/1‑100 (MQ=255)
                                                gagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCt                                               >  1:1801355/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAACAAAGAATATCTCCtat                                             >  2:1268429/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                             >  2:1782050/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                             >  1:1080639/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                             >  1:2996096/1‑100 (MQ=255)
                                                   cATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtata                                            >  2:1539145/1‑100 (MQ=255)
                                                    aTCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGATCGCGTTTTTTTGTCCCAAGTCATCCCCGTTACTGAGCAAAAAAAAGAATATCTCCTatat                                           <  1:40600/100‑1 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGTCAGAGAGCGCGGTTTTTTTTTTTTTTTTTTCCCCTTTAAGGAGCAAAAAAAATAATATCGGTTATATgtg                                        >  2:3035618/1‑98 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                        <  1:2939288/100‑1 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCACCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                        >  1:2811578/1‑100 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATAACCTTTACTGAGAAAAAAAAAGAATATATACTATATga                                         >  1:1115516/1‑99 (MQ=255)
                                                         aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGaa                                      <  1:2439094/100‑1 (MQ=255)
                                                          aGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAat                                     <  2:2101584/100‑1 (MQ=255)
                                                           gggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatc                                    >  2:1256644/1‑100 (MQ=255)
                                                           gggTAGGGGGATGCACAAAGAATGGGCAGAAAGCGCGGTTTTTTATCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATAAGAatc                                    >  2:941337/1‑100 (MQ=255)
                                                            ggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatca                                   >  1:1190391/1‑100 (MQ=255)
                                                             gTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcat                                  >  2:1897237/1‑100 (MQ=255)
                                                                    gATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgg                           >  1:286365/1‑100 (MQ=255)
                                                                      tGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgggg                         <  1:2526436/100‑1 (MQ=255)
                                                                       gCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATAggg                          >  1:2009048/1‑98 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTTTTTTAATTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  2:1452094/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  1:393291/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  2:730943/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  1:506410/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  2:1904487/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAACGCGGTTTTTTGTCAAAAGTCATCCCCTTAACTGAGCAAAAAAAAGAATATATCCTATATGAGAATCATCAAt                              >  2:1040769/1‑93 (MQ=255)
                                                                           aaaGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCCTCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaat                    <  2:2238403/100‑1 (MQ=255)
                                                                           aaaGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaa                     <  1:1876635/99‑1 (MQ=255)
                                                                                   ggCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGtttt             >  1:1303157/1‑99 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATTAATCGGGGTTAATAAGTTTTg            <  2:475741/100‑1 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            <  2:203418/100‑1 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            <  1:2085915/100‑1 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            >  2:2272331/1‑100 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTAACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGttt              <  2:2967764/98‑1 (MQ=255)
                                                                                     cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  1:85743/100‑1 (MQ=255)
                                                                                     cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          >  2:2882033/1‑100 (MQ=255)
                                                                                     cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  2:392749/100‑1 (MQ=255)
                                                                                      atagagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  2:2834898/97‑1 (MQ=255)
                                                                                      agagagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAATAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  1:2629639/99‑1 (MQ=255)
                                                                                      agagaaCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  2:2565652/99‑1 (MQ=255)
                                                                                        agagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcc       <  2:1563396/100‑1 (MQ=255)
                                                                                         gagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTGCGTaccc     >  2:2494022/1‑100 (MQ=255)
                                                                                          agCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTcccc     <  1:3024307/100‑1 (MQ=255)
                                                                                           gcgcGGTTTTTTGTTCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa    >  2:1021138/1‑100 (MQ=255)
                                                                                           gcgcGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa    >  1:2972683/1‑100 (MQ=255)
                                                                                            cgcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCag   >  1:2069211/1‑100 (MQ=255)
                                                                                             gcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCaga  >  1:2312165/1‑100 (MQ=255)
                                                                                                |                                                                                                
AGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGA  >  NC_000913/4352364‑4352556

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