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,4600TG97.9% 317.4 / ‑2.3 93intergenic (+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 (1/1);  new base G (51/40);  total (52/41)
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.12e-01

AGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAG  >  NC_000913/4352364‑4352555
                                                                                                |                                                                                               
agTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttt                                                                                              <  2:889516/100‑1 (MQ=255)
 ntACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttt                                                                                             >  2:781623/2‑100 (MQ=255)
 gTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttt                                                                                             >  1:128089/1‑100 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCGGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                            <  1:2125384/100‑1 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCCCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                            >  2:1651151/1‑100 (MQ=255)
   aaagagTCACGATATTTTTCATTATCACCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                            <  2:525979/97‑1 (MQ=255)
    cAGAGTCACGTTATTTTTCATTCTCTCCGCGATGGGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                          <  2:2169643/100‑1 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                            >  2:555229/1‑98 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                          <  2:16238/100‑1 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                          >  1:90997/1‑100 (MQ=255)
    cAGAGTCACGATATTTTTCATCCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                          <  1:1015849/100‑1 (MQ=255)
     agagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                          <  2:2510380/99‑1 (MQ=255)
      gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTc                                                                                        >  1:727130/1‑100 (MQ=255)
      gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTc                                                                                        <  2:90996/100‑1 (MQ=255)
        gTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTcc                                                                                       >  1:162977/1‑99 (MQ=255)
          cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCaa                                                                                    >  2:366193/1‑100 (MQ=255)
           aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAg                                                                                   >  1:75016/1‑100 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGt                                                                                  >  1:2074595/1‑100 (MQ=255)
                atTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAATGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATc                                                                              >  2:1660091/1‑100 (MQ=255)
                  ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                            >  2:1300265/1‑100 (MQ=255)
                  ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                            >  1:999923/1‑100 (MQ=255)
                  ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATccc                                                                            <  1:1159714/100‑1 (MQ=255)
                   tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATcccc                                                                           >  1:390276/1‑100 (MQ=255)
                     tCATTCTCTCCGAGATGTGATGCAGGAGAGCATCTGAAGGGTATGGGGATGCAAAAAGAATTGGCATAGTTCTCGGTTTTTTGTCACAAGTCATcca                                                                            <  2:300010/97‑2 (MQ=255)
                        ttCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAc                                                                      <  1:366197/100‑1 (MQ=255)
                           tctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGa                                                                   >  2:924713/1‑100 (MQ=255)
                            ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                  <  2:779078/100‑1 (MQ=255)
                             tcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                 <  2:1329796/100‑1 (MQ=255)
                             tcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                 <  1:982829/100‑1 (MQ=255)
                             tcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                 <  1:2221664/100‑1 (MQ=255)
                              ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                 >  2:856382/1‑99 (MQ=255)
                              ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                >  2:2446144/1‑100 (MQ=255)
                                gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                 <  2:2078130/97‑1 (MQ=255)
                                 cgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaa                                                             >  1:1700635/1‑100 (MQ=255)
                                  gATGTGATGCAGGAGAGAATCTGAGGGGTAGGGGGATGCCCAAAGAAGGGGCAAAGAGCGCGGTTTTTTGCCCCAGTTCATCCCCTTTACGGAGAaaaa                                                             <  2:1510395/99‑1 (MQ=255)
                                    tgtgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaa                                                          >  1:921521/1‑100 (MQ=255)
                                    tgtgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaa                                                          <  1:892022/100‑1 (MQ=255)
                                      tgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAg                                                        >  1:2289120/1‑100 (MQ=255)
                                      tgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAg                                                        <  1:2155378/100‑1 (MQ=255)
                                       gATGCAGGAGAGCATCTGAAGGGTAGGGGGATGGACAAAGAAAGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaaaaaaa                                                         >  2:1533564/1‑98 (MQ=255)
                                       gATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGa                                                       >  2:728797/1‑100 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                     >  2:881656/1‑100 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                     >  2:2020479/1‑100 (MQ=255)
                                           cAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAatat                                                   >  2:43846/1‑100 (MQ=255)
                                           cAGGAGAGCATCAGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAatat                                                   >  1:1687781/1‑100 (MQ=255)
                                               agagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGAGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTcc                                               <  2:1213383/100‑1 (MQ=255)
                                                gagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTTTTTTTTTTTTTCCCCTTTACTGAGCAAAAAAAAGAATATCTCCt                                              >  2:1081438/1‑100 (MQ=255)
                                                gagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTcc                                               <  2:935898/99‑1 (MQ=255)
                                                gagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCt                                              <  1:2272239/100‑1 (MQ=255)
                                                 agCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCta                                             >  1:1404705/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTTCCAAGTCATCCCCTTTATTTAGCAAAAAAAAGAATATCTCCtat                                            >  2:2119884/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                            >  1:1257680/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                            >  1:1162276/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                            >  1:1543052/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCta                                             <  2:1544550/99‑1 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGTGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAGAATATCTCCTatat                                          >  2:655217/1‑100 (MQ=255)
                                                  gCATCTGAAGGCTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCTCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                            >  2:1341157/1‑100 (MQ=255)
                                                   cATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtata                                           <  2:1355249/100‑1 (MQ=255)
                                                     tCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATg                                         >  2:898027/1‑100 (MQ=255)
                                                      cTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATga                                        <  1:924721/100‑1 (MQ=255)
                                                       tccgatctAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                       <  1:205355/93‑1 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                       >  1:12198/1‑100 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                       <  1:881664/100‑1 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                       <  1:1941462/100‑1 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCAAATAAAAAAAAAAAAAAAAGAATATCTCCTATATgag                                       <  2:9832/100‑1 (MQ=255)
                                                        gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATGTCCTATATgag                                       >  1:1040019/1‑99 (MQ=255)
                                                        gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATCTga                                        >  1:225791/1‑98 (MQ=255)
                                                           gggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatc                                   >  2:2059605/1‑100 (MQ=255)
                                                            ggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCCTCCCCTTTACTGAGCAAAAAAATGAATATCTCCTATATGAGAatc                                   >  2:1335552/1‑99 (MQ=255)
                                                              tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                >  1:553820/1‑100 (MQ=255)
                                                              tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                <  1:532649/100‑1 (MQ=255)
                                                              tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                >  2:224104/1‑100 (MQ=255)
                                                               aGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgagagatcgga                               >  2:205353/1‑93 (MQ=255)
                                                               aGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAAtcatca                               <  2:269854/100‑1 (MQ=255)
                                                               aGGGGGATGCACAAAGAATGGGCAGAAAGCGCGGTTTTTTGTCCCAAGTCATCCCATTTAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAAtcatca                               <  2:1696704/100‑1 (MQ=255)
                                                                gggggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCaa                              >  1:2355288/1‑100 (MQ=255)
                                                                gggggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAAtcatcan                              <  2:1924059/100‑2 (MQ=255)
                                                                 ggggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAATGAGAAAAAAAAAGAATATCTCCTATATGAGAATCATCAAt                             <  2:1811955/100‑1 (MQ=255)
                                                                  gggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGCATATCTCCTATATGAGAATCATCAAt                             >  2:1645652/1‑99 (MQ=255)
                                                                   ggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTTTTTTTTTTTTTCCCCTTAACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCCTCAAt                             >  2:649090/1‑98 (MQ=255)
                                                                   ggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCg                           >  1:1881087/1‑100 (MQ=255)
                                                                      tGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgggg                        <  1:2282196/100‑1 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                      >  1:767570/1‑100 (MQ=255)
                                                                         acaAGGAATTGGCAGAGATCTCGGTTTTTTGTCCGTAGTCATCCTCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCAGCAATCggg                         <  2:1404346/96‑1 (MQ=255)
                                                                          caAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaa                    >  1:2210784/1‑100 (MQ=255)
                                                                                 tGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGttt             >  1:104124/1‑100 (MQ=255)
                                                                                  gggCAGAGAGCGCGGTTTTTTTTTTTTTTTTATCCCCTATACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGtttt            >  2:2504143/1‑100 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg           <  1:1861556/100‑1 (MQ=255)
                                                                                    gCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATTAGCATCATCAATCGGGAGTAATAAGTTTTGc          >  2:1763383/1‑100 (MQ=255)
                                                                                    gCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg           <  2:1089927/99‑1 (MQ=255)
                                                                                    gCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCACCTTGACTAAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGc          <  2:511850/100‑1 (MQ=255)
                                                                                     cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg         <  2:2095308/100‑1 (MQ=255)
                                                                                     cAGAGAGCACAGTTTTTTGTCTCAAGTCATCACCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg         <  2:1703700/100‑1 (MQ=255)
                                                                                       tagagagCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTc       >  1:2474668/2‑100 (MQ=255)
                                                                                        agagCGCGGTTTTTTTTTCTAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATAATCAATCGGGGTTAATAAGTTTTGCGt        >  2:146711/1‑98 (MQ=255)
                                                                                         gagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTccc     <  1:801144/100‑1 (MQ=255)
                                                                                           gcgcGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa   >  1:535883/1‑100 (MQ=255)
                                                                                           gcgcGGTGTTATGTCCTGAGTCATCCCCTTTACTTAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGTGTTAATAAGTTTTGCGGCCCCa   <  2:875084/100‑1 (MQ=255)
                                                                                            cgcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCag  <  1:728804/100‑1 (MQ=255)
                                                                                            cgcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCag  >  2:170022/1‑100 (MQ=255)
                                                                                                |                                                                                               
AGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAG  >  NC_000913/4352364‑4352555

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