Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,007 G→A 7.3% 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 7.3% 210.1 / 4.0 68R29C (CGT→TGT) mazFendoribonuclease toxin MazF
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (3/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTG  >  NC_000913/2910858‑2911154
                                                                                                                                                     |                                                                                                                                                   
tACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACg                                                                                                                                                     <  2:949522/150‑1 (MQ=255)
    tGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTATATGAAAGGACTCAGGACAACAGCTGGACGATGt                                                                                                                                                 >  2:241949/1‑150 (MQ=255)
    tGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGt                                                                                                                                                 >  2:1296442/1‑150 (MQ=255)
          gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                           >  2:949326/1‑150 (MQ=255)
            tAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTg                                                                                                                                         >  2:947803/1‑150 (MQ=255)
            tAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCGGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                         <  2:461763/150‑3 (MQ=255)
             aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgtt                                                                                                                                        <  2:380623/150‑4 (MQ=255)
               cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                      >  2:542645/1‑150 (MQ=255)
               cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                      >  2:697769/1‑150 (MQ=255)
                     gCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTAcc                                                                                                                                <  2:270105/150‑1 (MQ=255)
                               tCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                      <  1:1012330/150‑1 (MQ=255)
                               tCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                      <  2:1387968/150‑1 (MQ=255)
                               tCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                      <  1:952607/150‑1 (MQ=255)
                                      cGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAAt                                                                                                               <  1:949326/150‑1 (MQ=255)
                                       ggATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATc                                                                                                              <  2:366065/150‑1 (MQ=255)
                                       ggATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATc                                                                                                              <  2:961873/150‑1 (MQ=255)
                                             aacaacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCa                                                                                                        >  1:683835/1‑150 (MQ=255)
                                              acaacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaa                                                                                                       >  2:655517/1‑150 (MQ=255)
                                                  cTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACTTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCa                                                                                                   <  1:379401/150‑1 (MQ=255)
                                                     cGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGaa                                                                                                >  2:150790/1‑149 (MQ=255)
                                                      gAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATc                                                                                               >  1:1348290/1‑150 (MQ=255)
                                                          ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                           <  1:408125/150‑1 (MQ=255)
                                                          ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                           >  1:215740/1‑150 (MQ=255)
                                                                 cTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCgg                                                                                    >  2:391430/1‑150 (MQ=255)
                                                                      gATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                               <  1:542645/150‑1 (MQ=255)
                                                                          gCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATc                                                                           <  1:697769/150‑1 (MQ=255)
                                                                                  cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                   <  1:655517/150‑1 (MQ=255)
                                                                                     ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccatt                                                                >  1:580565/1‑150 (MQ=255)
                                                                                          aaacaGACACATAACTGTTTTGTTTTTGTACAAGAAAAGAAAAAGGACAACAACTTAACAATATCCAGCTTGAAAGCTACCTTTTGTCGGGTAAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattaccat                                                           >  2:1306026/3‑149 (MQ=255)
                                                                                           cacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGa                                                          >  1:601511/1‑150 (MQ=255)
                                                                                               gACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                      >  2:1298852/1‑150 (MQ=255)
                                                                                                acacATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTcc                                                     <  1:539375/150‑1 (MQ=255)
                                                                                                 cacaTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                    >  1:368769/1‑150 (MQ=255)
                                                                                                     tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                <  1:404892/150‑1 (MQ=255)
                                                                                                     tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                >  2:571380/1‑150 (MQ=255)
                                                                                                     tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCCACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                <  1:769367/150‑1 (MQ=255)
                                                                                                      aCCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCTTATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCt                                               <  1:241123/150‑1 (MQ=255)
                                                                                                      aCCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCt                                               <  2:215740/150‑1 (MQ=255)
                                                                                                       ccTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTCGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCtt                                              >  1:1254514/1‑150 (MQ=255)
                                                                                                         tGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTc                                            <  1:385469/150‑1 (MQ=255)
                                                                                                         tGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTCGACGATGTCCAGCTTTCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTc                                            <  1:8514/150‑1 (MQ=255)
                                                                                                           ttttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCgg                                          <  2:683835/150‑1 (MQ=255)
                                                                                                            tttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGc                                         <  1:571380/150‑1 (MQ=255)
                                                                                                             ttgttgttgTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGct                                        >  1:644341/1‑150 (MQ=255)
                                                                                                               gttgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCACCCCAAATCAGATCGCCAATATCGGGTACCTATCGGCTTACCATTACCAGACTTCCTTAACTTTCGGctct                                      >  2:1308449/1‑150 (MQ=255)
                                                                                                                 tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                    >  2:567710/1‑150 (MQ=255)
                                                                                                                 tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                    <  1:947803/150‑1 (MQ=255)
                                                                                                                  gttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                   >  2:942827/1‑150 (MQ=255)
                                                                                                                    tgtACATGAAAGGACTCAGGACAACAGCTCGACGATCTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATAAGATCGCCCATATCGGGTACGAATCGGCTTACCATTACCAGACTTCCTTATCTTACCGCTCTCCCCa                                 >  1:998311/1‑150 (MQ=255)
                                                                                                                     gtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAg                                >  2:840938/1‑150 (MQ=255)
                                                                                                                     gtACATGAAAGGACTCAGGACAACAGCTGGAAGATGTCCAGCTTGCTCGCTACCTTTTCTCGCGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATATTTCGGCTCTCCCCAg                                >  1:497823/1‑150 (MQ=255)
                                                                                                                        cATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCg                             >  1:1410148/1‑150 (MQ=255)
                                                                                                                           gAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGata                          >  1:669341/1‑150 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                 <  2:1348290/150‑1 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                 >  1:5248/1‑150 (MQ=255)
                                                                                                                                        acaacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGt             <  2:1355896/150‑1 (MQ=255)
                                                                                                                                         caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtt            >  1:1049623/1‑150 (MQ=255)
                                                                                                                                         caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtt            >  2:1295968/1‑150 (MQ=255)
                                                                                                                                           caaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt          >  1:699439/3‑150 (MQ=255)
                                                                                                                                           acaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt          >  2:97294/1‑150 (MQ=255)
                                                                                                                                            cacctcgACGATGTCAAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACACAAATCAAAAAACCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTc         >  2:322153/7‑150 (MQ=255)
                                                                                                                                            caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTc         >  2:185775/1‑150 (MQ=255)
                                                                                                                                             aGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTcc        >  1:1288337/1‑150 (MQ=255)
                                                                                                                                               cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg      <  1:454562/150‑1 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg    >  1:53827/1‑150 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg    >  2:630440/1‑150 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg    <  2:5248/150‑1 (MQ=255)
                                                                                                                                                  gACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt   >  1:1009867/1‑150 (MQ=255)
                                                                                                                                                   aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg  <  2:580565/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTG  >  NC_000913/2910858‑2911154

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: ‑