Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,252 A→G 8.6% intergenic (+27/+34) 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,2520AG 8.6% 582.8 / 29.5 190intergenic (+27/+34)rpmE/yiiX50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base A (85/88);  new base G (9/8);  total (94/96)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

ACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCT  >  NC_000913/4127104‑4127400
                                                                                                                                                    |                                                                                                                                                    
         tCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttt                                                                                                                                             >  1:991507/1‑149 (MQ=255)
            aCGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttttt                                                                                                                                          <  2:742490/149‑1 (MQ=255)
               tgtgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                       <  1:1262536/149‑1 (MQ=255)
               tgtgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                       <  1:1521087/149‑1 (MQ=255)
                 tgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGc                                                                                                                                     <  2:991507/149‑1 (MQ=255)
                  gcagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGcc                                                                                                                                    <  1:801799/149‑1 (MQ=255)
                   cagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCt                                                                                                                                   >  3:1715334/1‑149 (MQ=255)
                     gcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGa                                                                                                                                 >  4:1262647/1‑149 (MQ=255)
                      caAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGaa                                                                                                                                >  1:118903/1‑149 (MQ=255)
                           gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                           <  4:97126/149‑1 (MQ=255)
                           gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                           <  3:1314704/149‑1 (MQ=255)
                           gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                           <  3:994449/149‑1 (MQ=255)
                           gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                           <  4:1581046/149‑1 (MQ=255)
                            ccACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTgg                                                                                                                          >  3:1419939/1‑149 (MQ=255)
                             cACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                         >  2:1649400/1‑149 (MQ=255)
                             cACCCGCTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCTCCGCACGGCGCTTTTTTTGTGCCTGAATTTTTGt                                                                                                                         <  3:1149332/149‑1 (MQ=255)
                              aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttt                                                                                                                                            >  1:1703338/1‑129 (MQ=255)
                              aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttt                                                                                                                                            <  2:1703338/129‑1 (MQ=255)
                              aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                        >  3:1492292/1‑149 (MQ=255)
                               cccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCg                                                                                                                       >  1:43792/1‑149 (MQ=255)
                                  gTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGTCCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAGTTTGGTCGTAt                                                                                                                    <  1:1431534/149‑1 (MQ=255)
                                     cttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAc                                                                                                                 <  2:1474359/149‑1 (MQ=255)
                                      ttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAcc                                                                                                                <  1:1649400/149‑1 (MQ=255)
                                              cAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttttttt                                                                                                                                         >  2:879982/1‑116 (MQ=255)
                                              cAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttttttt                                                                                                                                         <  1:879982/116‑1 (MQ=255)
                                                aaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAg                                                                                                      >  4:1416068/1‑149 (MQ=255)
                                                aaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAg                                                                                                      >  2:1208649/1‑149 (MQ=255)
                                                aaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAg                                                                                                      <  3:1580456/149‑1 (MQ=255)
                                                    gCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGc                                                                                                  <  4:69157/149‑1 (MQ=255)
                                                    gCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGc                                                                                                  >  1:252058/1‑149 (MQ=255)
                                                    gCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGc                                                                                                  >  1:222109/1‑149 (MQ=255)
                                                           gTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCtt                                                                                           >  2:1602759/1‑149 (MQ=255)
                                                            ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCa                                                                                                >  2:638127/1‑143 (MQ=255)
                                                            ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCa                                                                                                <  1:638127/143‑1 (MQ=255)
                                                             tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTg                                                                                         <  3:1444179/149‑1 (MQ=255)
                                                              gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                               <  4:1509998/126‑1 (MQ=255)
                                                               cTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGc                                                                                       >  2:1045117/1‑149 (MQ=255)
                                                                tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAt                                                                                               >  2:17985/1‑140 (MQ=255)
                                                                tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAt                                                                                               <  1:17985/140‑1 (MQ=255)
                                                                  ccGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAc                                                                                    >  3:442219/1‑149 (MQ=255)
                                                                      tGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtg                                                                                <  1:1602759/149‑1 (MQ=255)
                                                                      tGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtg                                                                                >  4:602138/1‑149 (MQ=255)
                                                                        gCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                              <  2:43792/149‑1 (MQ=255)
                                                                        gCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                              <  3:853233/149‑1 (MQ=255)
                                                                         ccGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggtg                                                                             >  1:283220/1‑149 (MQ=255)
                                                                           gtgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAg                                                                           <  4:1165291/149‑1 (MQ=255)
                                                                             gtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGtt                                                                         >  4:816629/1‑149 (MQ=255)
                                                                                gACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTtgtg                                                                      >  2:25433/1‑149 (MQ=255)
                                                                                gACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTtgtg                                                                      <  2:1392498/149‑1 (MQ=255)
                                                                                 aCCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                     >  4:421661/1‑149 (MQ=255)
                                                                                 aCCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                     >  3:660251/1‑149 (MQ=255)
                                                                                  ccGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgc                                                                    >  4:1592046/1‑149 (MQ=255)
                                                                                   cGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgcg                                                                   >  1:974510/1‑149 (MQ=255)
                                                                                    gCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGa                                                                                     >  3:1514177/1‑130 (MQ=255)
                                                                                     cTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCa                                                                 >  2:5061/1‑149 (MQ=255)
                                                                                      ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGGGTTGTGGCGCAt                                                                >  4:1562092/1‑149 (MQ=255)
                                                                                         aacaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGa                                                             >  4:1394359/1‑149 (MQ=255)
                                                                                           caaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaa                                                           >  2:1299046/1‑149 (MQ=255)
                                                                                            aaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                          >  3:1088817/1‑149 (MQ=255)
                                                                                            aaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                          >  1:528863/1‑149 (MQ=255)
                                                                                                 tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                     >  3:1297997/1‑149 (MQ=255)
                                                                                                 tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                     >  4:1627837/1‑149 (MQ=255)
                                                                                                  ttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCtt                                                    <  2:860685/149‑1 (MQ=255)
                                                                                                      aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                               >  3:503359/1‑86 (MQ=255)
                                                                                                      aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                               <  4:503359/86‑1 (MQ=255)
                                                                                                        aTCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCg                                              >  2:840230/1‑149 (MQ=255)
                                                                                                          cccGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                            >  1:410710/1‑117 (MQ=255)
                                                                                                          cccGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                            <  2:410710/117‑1 (MQ=255)
                                                                                                                cagcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGGTGTGGCGCATCGAAAACGGCTTGTGGCGAGaccacc                                      >  4:1466495/1‑149 (MQ=255)
                                                                                                                     aaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTc                                 <  3:480292/149‑1 (MQ=255)
                                                                                                                      aTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                                <  2:151319/149‑1 (MQ=255)
                                                                                                                       tAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCg                                                                                                    >  4:1304410/1‑80 (MQ=255)
                                                                                                                       tAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCt                               <  3:157411/149‑1 (MQ=255)
                                                                                                                        aaGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTc                              <  4:1089049/149‑1 (MQ=255)
                                                                                                                         aGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCt                               >  1:844069/1‑147 (MQ=255)
                                                                                                                          gACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCaa                            <  2:234932/149‑1 (MQ=255)
                                                                                                                           aCCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAg                           >  4:533894/1‑149 (MQ=255)
                                                                                                                               gATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGa                       >  1:1554511/1‑149 (MQ=255)
                                                                                                                                 tttCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAat                     >  3:241211/1‑149 (MQ=255)
                                                                                                                                     cGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTc                 <  3:623501/149‑1 (MQ=255)
                                                                                                                                     cGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTc                 <  4:241211/149‑1 (MQ=255)
                                                                                                                                      gAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAtttt                                                                                                                            <  2:660897/41‑1 (MQ=255)
                                                                                                                                      gAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCAt                                                                >  3:1400072/1‑101 (MQ=255)
                                                                                                                                      gAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt                >  1:1007834/1‑149 (MQ=255)
                                                                                                                                       aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt               <  3:1627837/149‑1 (MQ=255)
                                                                                                                                        aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGCAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTg              >  1:698597/1‑149 (MQ=255)
                                                                                                                                        aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGc                          >  1:1236379/1‑137 (MQ=255)
                                                                                                                                        aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGc                          <  2:1236379/137‑1 (MQ=255)
                                                                                                                                        aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTg              >  4:559106/1‑149 (MQ=255)
                                                                                                                                         aaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc             >  3:766050/1‑149 (MQ=255)
                                                                                                                                         aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc             <  1:661983/149‑1 (MQ=255)
                                                                                                                                         aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc             >  2:904253/1‑149 (MQ=255)
                                                                                                                                          aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCaa                            <  4:196878/133‑1 (MQ=255)
                                                                                                                                          aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCaa                            >  3:196878/1‑133 (MQ=255)
                                                                                                                                          aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAAAGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc            >  4:648133/1‑149 (MQ=255)
                                                                                                                                          aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc            <  1:295631/149‑1 (MQ=255)
                                                                                                                                            aGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                              >  2:291602/1‑81 (MQ=255)
                                                                                                                                               gccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc       >  2:1190157/1‑149 (MQ=255)
                                                                                                                                                  gcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                            <  3:1634790/77‑1 (MQ=255)
                                                                                                                                                    aCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCt  >  1:1049375/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCT  >  NC_000913/4127104‑4127400

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