Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,632 T→G 9.6% intergenic (+372/+48) 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,6320TG 9.6% 200.6 / 11.9 71intergenic (+372/+48)yjjV/yjjWputative DNase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base G (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

CGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTG  >  NC_000913/4614483‑4614781
                                                                                                                                                     |                                                                                                                                                     
cgcATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCt                                                                                                                                                       <  2:269646/150‑1 (MQ=255)
    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:111419/1‑150 (MQ=255)
                                                                                 gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      <  2:15391/150‑1 (MQ=255)
                                                                                 gCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTc                                                                      >  2:205722/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:714682/150‑1 (MQ=255)
                                                                                                           cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGc                                            >  1:714682/1‑150 (MQ=255)
                                                                                                           cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGc                                            >  2:194902/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                              >  2:1408106/1‑150 (MQ=255)
                                                                                                                         caATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGt                              >  1:1195936/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                          >  2:948107/1‑150 (MQ=255)
                                                                                                                             gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTc                          >  1:1346430/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)
                                                                                                                                                     |                                                                                                                                                     
CGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTG  >  NC_000913/4614483‑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: