Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,251 C→T 8.1% intergenic (+26/+35) rpmE → / ← yiiX 50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,2510CT 8.1% 257.1 / 9.7 85intergenic (+26/+35)rpmE/yiiX50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base C (33/45);  new base T (4/3);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAC  >  NC_000913/4127108‑4127395
                                                                                                                                               |                                                                                                                                                
gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGgc                                                                                                                                             >  4:675720/1‑149 (MQ=255)
    cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtt                                                                                                                                         >  2:862596/1‑149 (MQ=255)
    cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtt                                                                                                                                         >  1:408572/1‑149 (MQ=255)
     tggACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttt                                                                                                                                        >  4:606028/3‑149 (MQ=255)
         cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCttttttc                                                                                                                                    <  1:865708/149‑2 (MQ=255)
         cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCttttttc                                                                                                                                    <  1:862596/149‑2 (MQ=255)
           tgtgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                  <  2:235736/149‑1 (MQ=255)
           tgcgcagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                  <  2:656981/146‑1 (MQ=255)
               cagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCt                                                                                                                              >  2:497791/1‑149 (MQ=255)
                     gTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAttt                                                                                                                        <  3:653868/149‑1 (MQ=255)
                      tGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAtttt                                                                                                                       <  1:383341/149‑1 (MQ=255)
                       gCCACCCGTTCTTCACTGGCGAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  2:682724/149‑1 (MQ=255)
                       gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  1:110899/149‑1 (MQ=255)
                       gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  2:661280/149‑1 (MQ=255)
                       gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  4:103040/149‑1 (MQ=255)
                       gCCACCCGTTCTTCACTGGCAAACAGCGTGATGATGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  4:583301/149‑1 (MQ=255)
                        ccACCCGTTCTCCAATGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGAGCTTTTTTTGTGCCTGAATTTTgg                                                                                                                     >  3:774870/1‑149 (MQ=255)
                          aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      >  1:784589/1‑146 (MQ=255)
                          aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  2:784589/146‑1 (MQ=255)
                          aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                   >  4:193990/1‑149 (MQ=255)
                          aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                   >  3:104492/1‑149 (MQ=255)
                           cccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCg                                                                                                                  >  1:14007/1‑149 (MQ=255)
                            ccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGt                                                                                                                 <  2:110292/149‑1 (MQ=255)
                            ccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGt                                                                                                                 >  3:687409/1‑149 (MQ=255)
                                 cttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAc                                                                                                            <  3:347789/149‑1 (MQ=255)
                                   tcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                          <  1:367781/149‑1 (MQ=255)
                                    cACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCAc                                                                                                         >  4:460755/1‑149 (MQ=255)
                                         gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAttt                                                                                                                        <  4:460997/129‑1 (MQ=255)
                                         gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAttt                                                                                                                        >  3:460997/1‑129 (MQ=255)
                                           aaaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCGACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGATAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACa                                                                                                  <  1:731300/149‑1 (MQ=255)
                                                      tGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCt                                                                                       >  1:210854/1‑149 (MQ=255)
                                                      tGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCt                                                                                       >  4:196253/1‑149 (MQ=255)
                                                        ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                     >  2:301827/1‑149 (MQ=255)
                                                        ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                     <  2:210854/149‑1 (MQ=255)
                                                        ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                     <  3:273111/149‑1 (MQ=255)
                                                        ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                     <  1:832525/149‑1 (MQ=255)
                                                         tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTg                                                                                    <  2:235218/149‑1 (MQ=255)
                                                          gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                   <  1:455384/149‑1 (MQ=255)
                                                            tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                 >  3:510200/1‑149 (MQ=255)
                                                            tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                 >  2:521054/1‑149 (MQ=255)
                                                                 gTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgt                                                                            >  3:755411/1‑149 (MQ=255)
                                                                   ggCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtgg                                                                          <  1:338733/149‑1 (MQ=255)
                                                                    gCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAGAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                         <  4:166012/149‑1 (MQ=255)
                                                                      cGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                       <  4:104492/149‑1 (MQ=255)
                                                                         gtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGtt                                                                    <  4:367335/149‑1 (MQ=255)
                                                                           tGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTtgt                                                                  <  4:813210/149‑1 (MQ=255)
                                                                            gtctgCTTCAACACGCGTTTCACCCTCCCGGGTAGCAAATAAGACCAGATTTCCGAAAAAAGCGCCGTTCGGCGCTTCTTTTGTGCCTGTGTTTTGGTCGTATTACCACGAAAAAAGCGGCCATTCTTTGGCGACCGTGGGGAGTtgtg                                                                 <  4:57072/145‑1 (MQ=255)
                                                                               cGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgcg                                                              >  2:16361/1‑149 (MQ=255)
                                                                                  ttCAACAAGGGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCAt                                                           <  2:365244/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGTCGCATc                                                          <  4:444372/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                          <  4:672494/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                          <  4:687409/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                          <  4:755411/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                          <  1:107675/149‑1 (MQ=255)
                                                                                   tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                          <  2:385216/149‑1 (MQ=255)
                                                                                      acaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaa                                                       >  1:1027/1‑149 (MQ=255)
                                                                                        aaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                     >  3:293864/1‑149 (MQ=255)
                                                                                        aaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                     >  2:323961/1‑149 (MQ=255)
                                                                                         aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAAc                                                    <  3:664841/149‑1 (MQ=255)
                                                                                         aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAAc                                                    >  2:308439/1‑149 (MQ=255)
                                                                                          gCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                       <  4:476953/129‑1 (MQ=255)
                                                                                          gCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                       >  3:476953/1‑129 (MQ=255)
                                                                                             tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                >  2:334264/1‑149 (MQ=255)
                                                                                                     tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                        <  3:463381/149‑1 (MQ=255)
                                                                                                             agcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCg                                >  4:406894/1‑149 (MQ=255)
                                                                                                             agcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCg                                >  3:287581/1‑149 (MQ=255)
                                                                                                               caAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGtt                              >  1:97910/1‑149 (MQ=255)
                                                                                                                 aaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTc                            >  2:60783/1‑149 (MQ=255)
                                                                                                                  aTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                           <  4:121582/149‑1 (MQ=255)
                                                                                                                   tAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                           >  2:475427/1‑148 (MQ=255)
                                                                                                                   tAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                           <  1:475427/148‑1 (MQ=255)
                                                                                                                      gACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCaa                       <  3:739194/149‑1 (MQ=255)
                                                                                                                       aCCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      <  2:412432/53‑1 (MQ=255)
                                                                                                                       aCCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                      >  1:412432/1‑53 (MQ=255)
                                                                                                                         cGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCg                    >  3:872247/1‑149 (MQ=255)
                                                                                                                           gATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGa                  <  3:675720/149‑1 (MQ=255)
                                                                                                                                acaaaaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATAtt             >  3:139889/4‑149 (MQ=255)
                                                                                                                                 cGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTc            <  4:76790/149‑1 (MQ=255)
                                                                                                                                   aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt          <  4:853953/149‑1 (MQ=255)
                                                                                                                                   aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt          >  2:209729/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt          <  3:193990/149‑1 (MQ=255)
                                                                                                                                    aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTg         >  4:645027/1‑149 (MQ=255)
                                                                                                                                      aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc       >  4:368827/1‑149 (MQ=255)
                                                                                                                                      aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc       <  1:323961/149‑1 (MQ=255)
                                                                                                                                      aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc       <  1:308439/149‑1 (MQ=255)
                                                                                                                                           gccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc  <  1:301827/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAC  >  NC_000913/4127108‑4127395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: