Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,352,460 T→G 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,4600TG90.0% 297.4 / 15.5 91intergenic (+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 (2/7);  new base G (39/42);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCAGA  >  NC_000913/4352364‑4352556
                                                                                                |                                                                                                
agTACAGAGTCACGATATTTTTCTTTTTTTTCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttt                                                                                               >  1:1440056/1‑100 (MQ=255)
agTACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGtttt                                                                                               >  2:1285205/1‑100 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             <  1:193265/100‑1 (MQ=255)
  tACAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgttttt                                                                                             >  1:456317/1‑100 (MQ=255)
   aCAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttttt                                                                                            >  1:951801/1‑100 (MQ=255)
   aCAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGgtttttt                                                                                            >  2:1426278/1‑100 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           <  1:166462/100‑1 (MQ=255)
    cAGAGTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           >  2:1117614/1‑100 (MQ=255)
      gagTCGCGATATTTTGCATGGTCTCCGGGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTg                                                                                           <  1:1570725/98‑1 (MQ=255)
      gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTc                                                                                         <  2:862422/100‑1 (MQ=255)
      gagTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTc                                                                                         <  2:795770/100‑1 (MQ=255)
      gagTCACGATATTTTTCAATCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGTATGCACAAAGAATGGGCTGAGAGCGCGTTTTTTTGTc                                                                                         <  1:1827660/100‑1 (MQ=255)
        gTCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTccc                                                                                       <  2:100060/100‑1 (MQ=255)
         tCACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCa                                                                                      <  2:442659/100‑1 (MQ=255)
          cACGATATTTTTTATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGAGCGGTTTTTTGTcc                                                                                        >  1:267417/1‑97 (MQ=255)
          cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCaa                                                                                     <  2:411832/100‑1 (MQ=255)
          cACGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCaa                                                                                     >  1:200508/1‑100 (MQ=255)
           aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAg                                                                                    <  1:371715/100‑1 (MQ=255)
           aCGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAg                                                                                    <  1:930341/100‑1 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGt                                                                                   >  2:318228/1‑100 (MQ=255)
            cGATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGt                                                                                   >  2:914770/1‑100 (MQ=255)
             gATATTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTc                                                                                  <  2:684982/100‑1 (MQ=255)
               tatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCAt                                                                                >  2:1014571/1‑100 (MQ=255)
               tatTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCAt                                                                                >  2:446406/1‑100 (MQ=255)
                atTTTTCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATc                                                                               >  2:746677/1‑100 (MQ=255)
                  ttttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCAt                                                                                >  2:967621/1‑97 (MQ=255)
                   tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATcccc                                                                            <  1:238535/100‑1 (MQ=255)
                   tttCATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATcccc                                                                            <  2:617940/100‑1 (MQ=255)
                      cATTCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCttt                                                                         <  2:1504681/100‑1 (MQ=255)
                        ttCTCTCCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAc                                                                       >  2:447549/1‑100 (MQ=255)
                         tctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACt                                                                      <  2:470582/100‑1 (MQ=255)
                          ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCTCATGTCATCCCCTTTACTg                                                                     >  1:1119051/1‑100 (MQ=255)
                          ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                     >  2:1583979/1‑100 (MQ=255)
                          ctctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTg                                                                     <  1:1165536/100‑1 (MQ=255)
                            ctcCGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAg                                                                   <  2:1033624/100‑1 (MQ=255)
                             ttcgcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGc                                                                  <  1:1184398/98‑1 (MQ=255)
                              ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCCCAAAGACTGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                 <  1:160630/100‑1 (MQ=255)
                              ccGCGATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCa                                                                 >  1:23219/1‑100 (MQ=255)
                                gcgATGTGATTCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTCTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaa                                                               <  1:1280793/100‑1 (MQ=255)
                                gcgATGTGATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCaaa                                                               <  2:920894/100‑1 (MQ=255)
                                      tgATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAg                                                         >  1:1226139/1‑100 (MQ=255)
                                      acatcgcggAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAg                                                         <  1:1671890/93‑1 (MQ=255)
                                       gATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGa                                                        <  1:333832/100‑1 (MQ=255)
                                       gATGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGa                                                        >  2:1184705/1‑100 (MQ=255)
                                        aTGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGaa                                                       >  1:480560/1‑100 (MQ=255)
                                        aTGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGaa                                                       >  2:755976/1‑100 (MQ=255)
                                         tGCAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAat                                                      <  1:901008/100‑1 (MQ=255)
                                           cAGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAatat                                                    >  2:271324/1‑100 (MQ=255)
                                            aGGAGAGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATAtc                                                   >  2:1097444/1‑100 (MQ=255)
                                              gagaGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAATAAAAGAATAtctc                                                 <  2:1306736/100‑1 (MQ=255)
                                              gagaGCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATAtctc                                                 <  2:1871473/100‑1 (MQ=255)
                                               agagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTcc                                                <  1:332615/100‑1 (MQ=255)
                                               agagCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTcc                                                <  2:1092172/100‑1 (MQ=255)
                                                 agCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCta                                              <  1:271322/100‑1 (MQ=255)
                                                 agCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCta                                              <  2:402955/100‑1 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                             >  1:1478743/1‑100 (MQ=255)
                                                  gCATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtat                                             >  2:550151/1‑100 (MQ=255)
                                                   cATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtata                                            <  1:318226/100‑1 (MQ=255)
                                                   cATCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCtata                                            <  2:697243/100‑1 (MQ=255)
                                                    aTCTGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTatat                                           >  1:1653732/1‑100 (MQ=255)
                                                       tGAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgag                                        <  2:1093698/100‑1 (MQ=255)
                                                        gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgaga                                       >  1:1430578/1‑100 (MQ=255)
                                                        gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgaga                                       >  2:1510072/1‑100 (MQ=255)
                                                        gAAGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATgaga                                       >  2:1641020/1‑100 (MQ=255)
                                                         aaGGGTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGaa                                      >  1:1102935/1‑100 (MQ=255)
                                                           gggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAATTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatc                                    <  1:1708284/100‑1 (MQ=255)
                                                            ggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCGCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatca                                   >  2:1878133/1‑100 (MQ=255)
                                                            ggTAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatca                                   >  2:600607/1‑100 (MQ=255)
                                                              tAGGGGGATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAatcatc                                 <  1:250292/100‑1 (MQ=255)
                                                               aaggggATGCACAAAGAATGGACAGAAAAAAAGATTTTTTGTCCCAAGTCATCCCTTAAAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAatcatc                                 <  1:472414/97‑1 (MQ=255)
                                                                   ggATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCg                            <  2:772446/100‑1 (MQ=255)
                                                                    gATGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTCACTGAGCAAAAAAAAGAATATATCCTATATGAGAATAATAAATCgg                           >  2:962526/1‑100 (MQ=255)
                                                                      tGCACAAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgggg                         <  2:926126/100‑1 (MQ=255)
                                                                      tGCACAAAGAATGGGAAGAGAGCGCGTTTTTTTTTACCAAGTCATCCACTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCgg                           >  1:1159734/1‑98 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGTTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGTATATCTCCTATATGAGAATCATCAATCGGGGtt                       <  1:1285198/100‑1 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGACTCATCAATCGGGGtt                       >  1:314305/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  2:907171/1‑100 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       <  2:1457568/100‑1 (MQ=255)
                                                                        cacaAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGtt                       >  1:1426847/1‑100 (MQ=255)
                                                                          caAAGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAATGAATAGCTCCTATACGAGAATCATCAATCGGGGTtaa                     >  1:1314741/1‑100 (MQ=255)
                                                                            aaGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaata                   >  1:404741/1‑100 (MQ=255)
                                                                            aaGAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAATAAGAAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTtaata                   <  1:446404/100‑1 (MQ=255)
                                                                              gAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAg                 <  1:1432449/100‑1 (MQ=255)
                                                                              gAATGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTAAAAAAAAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAg                 <  1:1097443/100‑1 (MQ=255)
                                                                               aaTGGGCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGt                <  2:491960/100‑1 (MQ=255)
                                                                                  gggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGtttt             >  2:1224217/1‑100 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            >  1:803348/1‑100 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            <  2:1902455/100‑1 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            >  2:198980/1‑100 (MQ=255)
                                                                                   ggCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            <  1:1867116/100‑1 (MQ=255)
                                                                                    gCAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTg            >  1:654730/1‑99 (MQ=255)
                                                                                     cAGAGAGCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCg          <  1:550149/100‑1 (MQ=255)
                                                                                      agagagCGCGGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGt         >  2:716309/1‑100 (MQ=255)
                                                                                            cgcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCag   >  2:1360720/1‑100 (MQ=255)
                                                                                            cgcgGTTTTTTGTCCCAAGTCATCCCCTTTACTGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCag   <  2:1533655/100‑1 (MQ=255)
                                                                                             gcgTTTTTTTGTCCCAAGTCATCCCATTTAATGAGCAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCa  <  1:451606/99‑2 (MQ=255)
                                                                                             gcgGTTTTTTGTCACAAGTCCTCCCCTTTACAGAAAAAAAAAAAGAATATCTCCTATATGAGAATCATCAATCGGGGTTAATAAGTTTTGCGTCCCCaga  <  1:363772/100‑1 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: