Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,009 T→C 7.6% H28R (CAT→CGT) ‡ mazF ← endoribonuclease toxin MazF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,911,0090TC 7.6% 500.8 / 21.4 170H28R (CAT→CGT) ‡mazFendoribonuclease toxin MazF
Reads supporting (aligned to +/- strand):  ref base T (75/82);  new base C (8/5);  total (83/87)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

CTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGT  >  NC_000913/2910861‑2911158
                                                                                                                                                    |                                                                                                                                                     
 tGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGATGATGt                                                                                                                                                     >  2:327824/1‑150 (MQ=255)
       gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                               >  1:481939/1‑150 (MQ=255)
          aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGc                                                                                                                                            <  1:2468560/150‑1 (MQ=255)
          aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGc                                                                                                                                            <  1:2104470/150‑1 (MQ=255)
            cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                          <  1:2362548/150‑1 (MQ=255)
            cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                          >  1:800109/1‑150 (MQ=255)
             gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                         <  2:1226828/150‑1 (MQ=255)
             gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGAATGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                         <  2:283108/150‑1 (MQ=255)
              cTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGc                                                                                                                                        >  2:2013865/1‑150 (MQ=255)
              cTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGc                                                                                                                                        <  2:1615175/150‑1 (MQ=255)
              cTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGc                                                                                                                                        <  1:1859076/150‑1 (MQ=255)
                 cGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTAc                                                                                                                                     <  1:1572483/150‑1 (MQ=255)
                       cACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTg                                                                                                                               <  1:987552/150‑1 (MQ=255)
                          gTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCg                                                                                                                            <  2:1159786/150‑1 (MQ=255)
                            tCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                          <  2:1715909/150‑1 (MQ=255)
                               tGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCa                                                                                                                       >  2:1680657/1‑150 (MQ=255)
                                gACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGTAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCaa                                                                                                                      >  2:199659/1‑150 (MQ=255)
                                gACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCaa                                                                                                                      >  1:945133/1‑150 (MQ=255)
                                  catgATAAAACTACTTAGAACGGATATCCTTTTGATTGCGTTGTACAAGTAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCaaaa                                                                                                                    <  2:2448663/147‑1 (MQ=255)
                                    ggATAAAACAACTTCGAACGGATATCCTTTTGCTTGCGTTGTACAAGTAACACACAGACACTTACCTGTTTTGTTTTTGTACATGAAAGGACTCATGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTATGTCGGGTCAAAATc                                                                                                                  <  2:1297176/150‑1 (MQ=255)
                                    ggATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATc                                                                                                                  <  2:1390804/150‑1 (MQ=255)
                                    ggATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATc                                                                                                                  >  2:1479252/1‑150 (MQ=255)
                                       tAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAc                                                                                                               >  1:636690/1‑150 (MQ=255)
                                       tAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgtt                                                                                                                                            <  1:2204350/121‑4 (MQ=255)
                                        aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                              <  2:2244822/150‑1 (MQ=255)
                                        aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                              <  1:2351436/150‑1 (MQ=255)
                                         aaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAccc                                                                                                             >  2:747023/1‑150 (MQ=255)
                                            caacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaaa                                                                                                          >  2:720375/1‑150 (MQ=255)
                                                ttCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAg                                                                                                      <  1:545019/150‑1 (MQ=255)
                                                    aaCGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCg                                                                                                  <  2:309976/150‑1 (MQ=255)
                                                     aCGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGc                                                                                                 >  2:227026/1‑150 (MQ=255)
                                                      cGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCGAAATCAGATCGcc                                                                                                <  2:326146/150‑1 (MQ=255)
                                                      cGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGcc                                                                                                <  1:43078/150‑1 (MQ=255)
                                                        gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                              >  2:1195875/1‑150 (MQ=255)
                                                        gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                              <  1:2331528/150‑1 (MQ=255)
                                                        gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                              <  1:395730/150‑1 (MQ=255)
                                                          tatCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTAATGTTGTTGTACATCAAACGAATCAGGACAACACCTGGACGATGTCCACCTTGCTCGCTACCTTTTCTCCGGTCAAAATCAACCCAAATCAGATCGCCCata                                                                                            >  2:2167289/1‑150 (MQ=255)
                                                            tCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATc                                                                                          >  1:2262974/1‑150 (MQ=255)
                                                             ccTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCg                                                                                         >  1:197306/1‑150 (MQ=255)
                                                              cTTTTGATTGCTTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGATTGCTCGCTACCCTTTGTCGGGTCAAAATCAACACAAATCAGATCGCCCATATCgg                                                                                        >  2:1129736/1‑150 (MQ=255)
                                                                 ttGATTGCGTTGTACAAGGCACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGCCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCACCCCAAATCAGATCGCCCATATCGGGTa                                                                                     >  2:257705/1‑150 (MQ=255)
                                                                 ttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                     <  2:2412019/150‑1 (MQ=255)
                                                                 ttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                     <  2:1829866/150‑1 (MQ=255)
                                                                 ttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                     <  1:1636853/150‑1 (MQ=255)
                                                                   gATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                                   >  2:1065615/1‑150 (MQ=255)
                                                                    aTTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGt                                                                                  <  1:654186/150‑1 (MQ=255)
                                                                     ttGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTa                                                                                 >  2:81163/1‑150 (MQ=255)
                                                                         gTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCgg                                                                             <  2:1284113/150‑1 (MQ=255)
                                                                         gTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCgg                                                                             <  2:451693/150‑1 (MQ=255)
                                                                             tACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTTGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCtta                                                                         >  1:2030812/1‑150 (MQ=255)
                                                                              aCAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttac                                                                        >  1:428897/1‑150 (MQ=255)
                                                                              aCAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttac                                                                        <  2:291000/150‑1 (MQ=255)
                                                                               cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                       <  1:720375/150‑1 (MQ=255)
                                                                               cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                       >  1:2443870/1‑150 (MQ=255)
                                                                               cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                       <  1:1050136/150‑1 (MQ=255)
                                                                                aaGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacca                                                                      >  1:1350130/1‑150 (MQ=255)
                                                                                 aGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccat                                                                     <  1:1059296/150‑1 (MQ=255)
                                                                                  ggAACACACAGACACATAGGTGTTTTGTTGTTGTAGATGAAAGGACTCAGGACAAGAGCGGGACGATGTCGAGCTTGCTCGCTAGCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccatt                                                                    <  2:1005868/150‑1 (MQ=255)
                                                                                  ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGt                                                                                      <  2:921080/132‑1 (MQ=255)
                                                                                  ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGt                                                                                      >  1:921080/1‑132 (MQ=255)
                                                                                     acacacAGACACATACCTGTTTAGTTGATGTACATGAAACTACTCAGGACAACAGATGAACGATGTACAGCTTGCTCGATACCTTTTGTCGGGTAAAAATCAAACAAAATAAGATCGCCCATATCCGGTACATATCGGCttaccattaac                                                                 >  2:240126/1‑148 (MQ=255)
                                                                                      cacacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCATGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacca                                                                >  2:1487315/1‑150 (MQ=255)
                                                                                      cacacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacca                                                                >  2:1578679/1‑150 (MQ=255)
                                                                                       acacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAg                                                               <  1:741157/150‑1 (MQ=255)
                                                                                        cacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTAGCTCGCTACCTTTTGTCGGGTAAAAATCAACCCAAATCAGATAGCCCATATCGGGTACGTATAGGCTTACCATTACCAGa                                                              >  2:1417447/1‑150 (MQ=255)
                                                                                        cacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCACCTTGCTCGCTACCTTTTGTCGCGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGa                                                              >  2:826613/1‑150 (MQ=255)
                                                                                              cacaTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                        <  2:117782/150‑1 (MQ=255)
                                                                                                caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCt                                                                           <  1:1142435/129‑1 (MQ=255)
                                                                                                caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTcc                                                         <  2:2284938/147‑1 (MQ=255)
                                                                                                  tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                    <  2:437606/150‑1 (MQ=255)
                                                                                                  tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                    <  2:2453141/150‑1 (MQ=255)
                                                                                                  tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                    >  1:2343279/1‑150 (MQ=255)
                                                                                                        ttttGTTGTTGTACATTAAAGGACTGAGGACAACAGCTGGACGATGTGCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTTTCGGCTTACCATTACCAGACTTCCTTATCTTTCgg                                              <  1:160112/150‑1 (MQ=255)
                                                                                                        ttttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCgg                                              <  1:826613/150‑1 (MQ=255)
                                                                                                           tgttgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGctc                                           <  2:945133/150‑1 (MQ=255)
                                                                                                            gttgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGctct                                          >  1:122537/1‑150 (MQ=255)
                                                                                                              tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                        <  2:353145/150‑1 (MQ=255)
                                                                                                              tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                        >  2:183682/1‑150 (MQ=255)
                                                                                                              tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                        >  1:2412599/1‑150 (MQ=255)
                                                                                                               gttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                       >  1:593040/1‑150 (MQ=255)
                                                                                                                      aTGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGa                                >  1:894580/1‑150 (MQ=255)
                                                                                                                       tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                               >  1:318733/1‑150 (MQ=255)
                                                                                                                          aaGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATAtt                            >  2:481564/1‑150 (MQ=255)
                                                                                                                            ggACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATtct                          <  1:1479252/150‑1 (MQ=255)
                                                                                                                                 cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                     <  2:602316/150‑1 (MQ=255)
                                                                                                                                    gACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAgg                  <  1:199659/150‑1 (MQ=255)
                                                                                                                                       gcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttt               >  1:768626/4‑150 (MQ=255)
                                                                                                                                       gcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttt               >  1:2257883/4‑150 (MQ=255)
                                                                                                                                       gcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttt               >  2:2152752/4‑150 (MQ=255)
                                                                                                                                        caaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt              >  2:1082846/3‑150 (MQ=255)
                                                                                                                                           gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg           <  1:1597398/150‑1 (MQ=255)
                                                                                                                                           gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg           >  1:1224402/1‑150 (MQ=255)
                                                                                                                                            cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg          <  2:428897/150‑1 (MQ=255)
                                                                                                                                            cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg          >  1:1286826/1‑150 (MQ=255)
                                                                                                                                              ggACGATGTCCAGCTTGCTCGCTACCTTTTGTAGGGTCAAAATAAACCCAAATCAGCTCGACCATATCGGGTACGTATCGGATTAACATTACGAGACTTCCTTATCTTTCGGCTGTCCCCAGACGATATACTCGGGGAGGTTTTccgccg        >  1:1746430/1‑150 (MQ=255)
                                                                                                                                               gACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt       <  1:1458391/150‑1 (MQ=255)
                                                                                                                                               gACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATAAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt       >  2:2339450/1‑150 (MQ=255)
                                                                                                                                                aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg      <  2:593040/150‑1 (MQ=255)
                                                                                                                                                 cGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGa     >  2:1723080/1‑150 (MQ=255)
                                                                                                                                                  gATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGAt    <  1:1866518/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGT  >  NC_000913/2910861‑2911158

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: