Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,633 T→C 9.8% intergenic (+373/+47) yjjV → / ← yjjW putative DNase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,6330TC 9.8% 190.7 / 9.6 71intergenic (+373/+47)yjjV/yjjWputative DNase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base C (5/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTG  >  NC_000913/4614487‑4614781
                                                                                                                                                  |                                                                                                                                                    
tCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTGTCGTCTTATCATGCCTACAGGTTTGTTACGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACTCTTGACGCGTCTTATc                                                                                                                                                   <  2:732439/150‑1 (MQ=14)
tCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                                                                                   <  2:384682/150‑1 (MQ=255)
  cGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                                                                                 <  1:1179880/150‑1 (MQ=255)
  cGGCAGTCGGCGCATAATGCCTGATGAGACGCTTGTCTCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                                                                                 >  2:941140/1‑150 (MQ=255)
     cAGTGGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTGTGTGACGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  2:994727/150‑1 (MQ=255)
     cAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  1:1363935/150‑1 (MQ=255)
     cAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  1:1396433/150‑1 (MQ=255)
     cAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  1:294703/150‑1 (MQ=255)
        tCGGCGCATAATGGATGATGCGACGCTTGTCGAGTCTTATCATGCGTAGAGGTTTGTGCCGATGCGTAGGCCGGATAAGGCGTTGGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACTCTTGACGCGTCTTATGTGGCCTAc                                                                                                                                           <  2:483150/150‑1 (MQ=14)
          ggCGCACAATGCCTGATGCGACGCTTTTCGCGTCTTATCTTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcattgtg                                                                                                                                         <  2:410538/150‑7 (MQ=11)
           gcgcATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                                        <  1:1312803/150‑1 (MQ=255)
              cATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCt                                                                                                                                     <  1:613609/150‑1 (MQ=255)
                       tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAAc                                                                                                                            <  1:929585/150‑1 (MQ=255)
                        gATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAGGGGTCTTATAAGGCCTACAGGTGTGTGCCGAAcc                                                                                                                           >  2:272057/1‑150 (MQ=14)
                           gCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTa                                                                                                                        <  1:272057/150‑1 (MQ=255)
                                gCTTGTCGCGGCTTATCCTGCCTACAGGTGTGGGCCGACCCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                                                                                   <  1:622880/150‑1 (MQ=14)
                                        cgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                                                                           <  1:228904/150‑1 (MQ=255)
                                          tCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCAt                                                                                                         <  2:1461175/150‑1 (MQ=255)
                                             tATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCAGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACTCTGTGCCGAACAGTAGGCCGTATCCGGCATGTc                                                                                                      >  2:472775/1‑150 (MQ=255)
                                                cATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTcaca                                                                                                   >  2:1359228/1‑150 (MQ=255)
                                                    ccTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATa                                                                                               >  2:282574/1‑150 (MQ=255)
                                                      tACAGGTTTGTGCCGATCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCCGCAGTTGGCGCAAAATGCCTGATGCGACGCTAGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATagc                                                                                             >  1:986184/1‑149 (MQ=255)
                                                                tgCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGaaa                                                                                   >  2:1268629/1‑150 (MQ=255)
                                                                 gCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAat                                                                                  >  2:1099615/1‑150 (MQ=255)
                                                                   cGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAatat                                                                                <  1:976793/150‑1 (MQ=255)
                                                                    gAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAatata                                                                               >  1:806326/1‑149 (MQ=255)
                                                                      aCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCaa                                                                             <  1:67725/150‑1 (MQ=255)
                                                                      aCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCaa                                                                             >  2:41400/1‑150 (MQ=255)
                                                                         gTAGGCCGTATTAGGCGTTCGCGCCGCATCGGTCAGTTGGCGCACAATGCGTGATTCGACGATTGACGAGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCg                                                                          <  2:42887/150‑1 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTAGCAGGCCTACAAGTCTGTGCCGAGCCGTAGGCCGTATCCGGCATGTCCCACATAGAGCGCCGGCAATATCAACCGGCt                                                                       >  2:508900/1‑150 (MQ=18)
                                                                             gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      <  2:15391/150‑1 (MQ=255)
                                                                             gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      >  2:205722/1‑150 (MQ=255)
                                                                             gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      >  2:111419/1‑150 (MQ=255)
                                                                             gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      >  1:524216/1‑150 (MQ=255)
                                                                                ggATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCAcc                                                                   <  2:334962/150‑1 (MQ=255)
                                                                                   tAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                                                                                                            <  1:875744/74‑1 (MQ=255)
                                                                                      ggCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCgcgc                                                             <  1:1400105/150‑1 (MQ=255)
                                                                                       gCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCa                                                            <  2:1372075/150‑1 (MQ=255)
                                                                                        cGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                              <  1:723716/67‑1 (MQ=25)
                                                                                         gTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGcc                                                                                                                                <  1:354774/80‑1 (MQ=255)
                                                                                            cgcgCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCttt                                                       >  1:805296/1‑150 (MQ=255)
                                                                                             gcgcCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTa                                                      >  2:1532184/1‑150 (MQ=255)
                                                                                             gcgcCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTa                                                      <  2:155911/150‑1 (MQ=255)
                                                                                               gccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAAc                                                    <  1:1278400/150‑1 (MQ=255)
                                                                                                ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACg                                                   <  1:282264/150‑1 (MQ=255)
                                                                                                    aTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCAta                                               >  2:145670/1‑149 (MQ=255)
                                                                                                    aTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATc                                               <  1:946866/150‑1 (MQ=255)
                                                                                                       cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGc                                            >  2:194902/1‑150 (MQ=255)
                                                                                                       cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGc                                            <  2:714682/150‑1 (MQ=255)
                                                                                                       cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGc                                            >  1:714682/1‑150 (MQ=255)
                                                                                                              tGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCATGACTACAAGTCTGTGCAGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGc                                     <  1:1482722/150‑1 (MQ=255)
                                                                                                              tGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGc                                     <  1:662606/150‑1 (MQ=255)
                                                                                                                   cacaATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAAc                                >  1:479918/1‑150 (MQ=255)
                                                                                                                    acaATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACg                               >  2:822125/1‑150 (MQ=255)
                                                                                                                    acaATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACg                               <  1:1058919/150‑1 (MQ=255)
                                                                                                                     caATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGt                              >  1:1195936/1‑150 (MQ=255)
                                                                                                                     caATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGt                              >  2:1408106/1‑150 (MQ=255)
                                                                                                                      aaTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTc                             <  2:805296/150‑1 (MQ=255)
                                                                                                                       aggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCt                            >  2:644234/3‑150 (MQ=255)
                                                                                                                        ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCtt                           >  1:1100567/2‑150 (MQ=255)
                                                                                                                        ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCtt                           >  1:1439921/2‑150 (MQ=255)
                                                                                                                         gacTGATGCGACGCTTGACGCGTCTTATGAGGCCTACAAGTGTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGTGCGACGGAAATATCAACCGGCTCAGGCCTCGCACCTTTAACGCATCAGCCTACGGCTCAACGTCTTc                          <  2:986184/148‑1 (MQ=38)
                                                                                                                         gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTc                          >  1:1346430/1‑150 (MQ=255)
                                                                                                                         gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTc                          >  2:948107/1‑150 (MQ=255)
                                                                                                                          ccTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTcc                         >  1:1369387/1‑150 (MQ=255)
                                                                                                                          ccTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTcc                         >  1:550807/1‑150 (MQ=255)
                                                                                                                             gATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGc                      >  1:321759/1‑150 (MQ=255)
                                                                                                                                 cGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAAAGTCTTCCGGCGTgg                  >  2:1464988/1‑150 (MQ=255)
                                                                                                                                  gACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGgc                 <  2:1346430/150‑1 (MQ=255)
                                                                                                                                             cgTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGc      >  2:1007587/1‑150 (MQ=255)
                                                                                                                                              gTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCt     >  2:1210823/1‑150 (MQ=255)
                                                                                                                                                 ttATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTg  <  2:287407/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTG  >  NC_000913/4614487‑4614781

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