Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,007 G→A 8.4% R29C (CGT→TGT)  mazF ← endoribonuclease toxin MazF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,911,0070GA 8.4% 284.5 / 10.6 95R29C (CGT→TGT) mazFendoribonuclease toxin MazF
Reads supporting (aligned to +/- strand):  ref base G (40/47);  new base A (3/5);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGA  >  NC_000913/2910858‑2911155
                                                                                                                                                     |                                                                                                                                                    
tACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACg                                                                                                                                                      >  1:261635/1‑150 (MQ=255)
 aCCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGa                                                                                                                                                     >  1:1368607/1‑150 (MQ=255)
  ccTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGAt                                                                                                                                                    >  2:313790/1‑150 (MQ=255)
    tGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGt                                                                                                                                                  <  1:714328/150‑1 (MQ=255)
    tGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGt                                                                                                                                                  >  1:238868/1‑150 (MQ=255)
      aTCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTcc                                                                                                                                                <  2:1263859/150‑1 (MQ=255)
          gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                            <  1:806734/150‑1 (MQ=255)
          gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                            <  2:328756/150‑1 (MQ=255)
           cTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTAGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtg                                                                                                                                           <  2:142817/150‑2 (MQ=255)
            tAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                          <  1:270722/150‑3 (MQ=255)
            tAACGCGACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                          <  2:853474/150‑3 (MQ=255)
             aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGc                                                                                                                                         >  1:964803/1‑150 (MQ=255)
             aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgtt                                                                                                                                         <  2:1466834/150‑4 (MQ=255)
               cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgttgt                                                                                                                                       <  1:614422/150‑6 (MQ=255)
                         tCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCtttt                                                                                                                             <  1:228224/150‑1 (MQ=255)
                              ttCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCgg                                                                                                                        <  1:573/150‑1 (MQ=255)
                                ccTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGt                                                                                                                      <  2:673696/150‑1 (MQ=255)
                                ccTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGt                                                                                                                      <  1:327238/150‑1 (MQ=255)
                                 cTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTc                                                                                                                     >  2:1392564/1‑150 (MQ=255)
                                        gATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCa                                                                                                              <  2:1536767/150‑1 (MQ=255)
                                        gATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCa                                                                                                              >  1:1457201/1‑150 (MQ=255)
                                        gATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCa                                                                                                              <  2:1324180/150‑1 (MQ=255)
                                         aTAAAACAACTTCGAACGGATATGCGTTAGATTTCCTTGTACAAGGAACACACAGACACAAACCAGTAATGATGTAGTAAATGAAAGGACTCAGGAAAACAGCTGGACGATGTCCAGAAACCTCGCTACCTTTTCTCGGGTCAAAAACaa                                                                                                             >  2:1347130/1‑150 (MQ=255)
                                          tAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAc                                                                                                            <  2:619196/150‑1 (MQ=255)
                                          tAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAc                                                                                                            <  2:1320968/150‑1 (MQ=255)
                                           aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTCTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                           <  1:769111/150‑1 (MQ=255)
                                            aaaCAACTTCGAACCGATATCCTTTTGATTGCGATGAATTAGGAACACACAGACACATACATGTATAGTTGATCTACATCAAAGGACTCACGACAACAGCAGCACGAAGTCCAGAATGCTCGCTACCTTTTGTCGGGTCAAAATCAAccc                                                                                                          >  1:1390070/1‑150 (MQ=255)
                                              acaacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaa                                                                                                        >  1:228835/1‑150 (MQ=255)
                                               caacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaaa                                                                                                       >  2:411538/1‑150 (MQ=255)
                                               caacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaaa                                                                                                       >  2:595772/1‑150 (MQ=255)
                                                 acTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                      >  2:621328/1‑117 (MQ=255)
                                                 acTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                      <  1:621328/117‑1 (MQ=255)
                                                  cTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCa                                                                                                    <  1:1187274/150‑1 (MQ=255)
                                                          ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                            <  2:837585/150‑1 (MQ=255)
                                                           tatatCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATGAACCCAAATCAGATCGCCCa                                                                                           <  2:286961/149‑1 (MQ=255)
                                                           gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                           <  1:1136926/150‑1 (MQ=255)
                                                               tCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATc                                                                                       >  2:1221932/1‑150 (MQ=255)
                                                               tCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATc                                                                                       <  2:1021766/150‑1 (MQ=255)
                                                                     tGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTATGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTAc                                                                                 >  1:1487252/1‑150 (MQ=255)
                                                                      gATTGCGTTGTACAAGGAACACACAGACACATACCTTTTTTGTTGTTGTACATGAAAGGACTGAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTTTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                                <  2:1381434/150‑1 (MQ=255)
                                                                      gATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                                <  2:238868/150‑1 (MQ=255)
                                                                      gATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                                >  1:1412229/1‑150 (MQ=255)
                                                                           cGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCg                                                                           <  2:263107/150‑1 (MQ=255)
                                                                           cGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCg                                                                           >  1:254303/1‑150 (MQ=255)
                                                                               gTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCtt                                                                       >  2:988161/1‑150 (MQ=255)
                                                                               gTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCtt                                                                       <  1:647636/150‑1 (MQ=255)
                                                                                 aCAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttac                                                                     <  1:929392/150‑1 (MQ=255)
                                                                                 aCAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttac                                                                     >  2:10467/1‑150 (MQ=255)
                                                                                    aGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccat                                                                  >  1:130020/1‑150 (MQ=255)
                                                                                     ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccatt                                                                 >  1:14686/1‑150 (MQ=255)
                                                                                     ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACATATCGGCttaccatt                                                                 <  2:964803/150‑1 (MQ=255)
                                                                                       aaCACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattac                                                               >  2:352817/1‑150 (MQ=255)
                                                                                        acacacAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacc                                                              <  1:313790/150‑1 (MQ=255)
                                                                                        acacacAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacc                                                              <  1:145747/150‑1 (MQ=255)
                                                                                            acaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGAc                                                          <  1:411538/150‑1 (MQ=255)
                                                                                             caGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACt                                                         >  1:662567/1‑150 (MQ=255)
                                                                                               gACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                       <  2:261635/150‑1 (MQ=255)
                                                                                                acacATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTcc                                                      <  2:1368607/150‑1 (MQ=255)
                                                                                                   caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTa                                                   <  1:1221932/150‑1 (MQ=255)
                                                                                                   caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTa                                                   >  2:491051/1‑150 (MQ=255)
                                                                                                      aCCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCt                                                <  2:1457201/150‑1 (MQ=255)
                                                                                                       ccTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCtt                                               <  2:435513/150‑1 (MQ=255)
                                                                                                         tGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTc                                             <  2:112808/150‑1 (MQ=255)
                                                                                                          gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                  <  1:865114/112‑1 (MQ=255)
                                                                                                          gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                  >  2:865114/1‑112 (MQ=255)
                                                                                                            tttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGc                                          <  1:501979/150‑1 (MQ=255)
                                                                                                              tgttgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGctc                                        >  1:248673/1‑150 (MQ=255)
                                                                                                                 tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                     <  1:595772/150‑1 (MQ=255)
                                                                                                                   ttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcccc                                   <  1:491051/150‑1 (MQ=255)
                                                                                                                       aCATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                               >  2:452676/1‑150 (MQ=255)
                                                                                                                         aTGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGa                             >  2:310415/1‑150 (MQ=255)
                                                                                                                          tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATTGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                            <  1:413354/150‑1 (MQ=255)
                                                                                                                          tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                            >  1:1415377/1‑150 (MQ=255)
                                                                                                                          tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                            >  1:1114817/1‑150 (MQ=255)
                                                                                                                                  cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGt                    >  1:931210/1‑150 (MQ=255)
                                                                                                                                   tCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTg                   <  2:647760/150‑1 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                  >  1:1203162/1‑150 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                  >  2:1333894/1‑150 (MQ=255)
                                                                                                                                     aGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGa                 <  2:931210/150‑1 (MQ=255)
                                                                                                                                      ggACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAg                <  2:890818/150‑1 (MQ=255)
                                                                                                                                       gACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAgg               >  1:291802/1‑150 (MQ=255)
                                                                                                                                        gagcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGt              >  1:615934/6‑150 (MQ=255)
                                                                                                                                        acaacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGt              >  2:109289/1‑150 (MQ=255)
                                                                                                                                         caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATtctc                      <  1:67627/141‑1 (MQ=255)
                                                                                                                                         caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtt             >  2:866043/1‑150 (MQ=255)
                                                                                                                                         caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTAGGCTCTCCCCAGTCGATATtctc                      >  2:67627/1‑141 (MQ=255)
                                                                                                                                           caaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt           >  2:983425/3‑150 (MQ=255)
                                                                                                                                            caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTc          <  2:130020/150‑1 (MQ=255)
                                                                                                                                             aGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTcc         >  2:1388208/1‑150 (MQ=255)
                                                                                                                                              gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg        >  2:664936/1‑150 (MQ=255)
                                                                                                                                              gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg        >  1:169237/1‑150 (MQ=255)
                                                                                                                                               cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg       <  1:988161/150‑1 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg     >  2:618805/1‑150 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg     <  1:1228802/150‑1 (MQ=255)
                                                                                                                                                  gACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt    <  1:925781/150‑1 (MQ=255)
                                                                                                                                                    cGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGa  >  2:511455/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGA  >  NC_000913/2910858‑2911155

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