Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,009 T→C 7.2% H28R (CAT→CGT) ‡ mazF ← mRNA interferase toxin, antitoxin is MazE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,911,0090TC 7.2% 371.4 / 14.2 127H28R (CAT→CGT) ‡mazFmRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base T (65/52);  new base C (3/6);  total (69/58)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 127 total aligned reads displayed.

TGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGAT  >  NC_000913/2910862‑2911156
                                                                                                                                                   |                                                                                                                                                   
tGATCAGCTAACGCTACGACATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCGTTTGATTGCGTTGTCCAAGGAACACACCGACGCATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACCACAGCTGCACGCTGt                                                                                                                                                   >  2:70240/1‑150 (MQ=255)
 gATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTc                                                                                                                                                  >  1:1324777/1‑150 (MQ=255)
      gCTAACGCTACGCCATCACGTTCCTGACCTGATAAAACAACTTCGAACGGATATCCTTTTGATTGAGTTGTAGAATGAACACACAGACTAAGACCAGTTTTGTTGTTGTACATGTATGTACTCATGTCCACAGCTGGACGATGTcctggt                                                                                                                                             >  2:890341/1‑146 (MQ=255)
      gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACATCTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                             >  1:268702/1‑150 (MQ=255)
      gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCt                                                                                                                                             >  1:1480217/1‑150 (MQ=255)
      gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTGTGTTGTTGTACATTAAAGGCCTCAGGACAAACGCTTGAAGATGTCCAGCt                                                                                                                                             >  1:912769/1‑150 (MQ=255)
        tAACGCTGCGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                           <  1:1749242/150‑3 (MQ=255)
        tAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                           <  1:1734694/150‑3 (MQ=255)
         aaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGc                                                                                                                                          >  1:1622168/1‑150 (MQ=255)
         aaCGCTACGCCATCAAGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgtt                                                                                                                                          <  2:1415563/150‑4 (MQ=255)
          aCGCTACGCCATCACTTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCt                                                                                                                                         <  1:307389/150‑1 (MQ=255)
           cGCTACGCCATCACGTTCCTGCCCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                        <  2:758544/150‑1 (MQ=255)
            gcaaCGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                       <  1:816254/147‑1 (MQ=255)
            gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                       <  2:1392307/150‑1 (MQ=255)
            gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                       <  2:1673057/150‑1 (MQ=255)
            gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                       <  1:858915/150‑1 (MQ=255)
             cTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGc                                                                                                                                      >  1:1544881/1‑150 (MQ=255)
             cTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGATTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGc                                                                                                                                      >  2:1106097/1‑150 (MQ=255)
                  ccATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACCCACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCt                                                                                                                                 >  1:1086383/1‑150 (MQ=255)
                  ccATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTCTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCt                                                                                                                                 >  2:1079828/1‑150 (MQ=255)
                    aTCTCGTTCCTGACCGGATAAAACAACTTCGAACGTATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCttt                                                                                                                               <  1:388952/150‑1 (MQ=255)
                      cACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTg                                                                                                                             <  1:929505/150‑1 (MQ=255)
                           tCCTGCCCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                        <  1:1778566/150‑1 (MQ=255)
                           tCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                        <  2:154842/150‑1 (MQ=255)
                             cTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTc                                                                                                                      <  2:1618681/150‑1 (MQ=255)
                              tGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCa                                                                                                                     >  1:436917/1‑150 (MQ=255)
                              tGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCa                                                                                                                     <  1:791812/150‑1 (MQ=255)
                               gACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTAGGGTCaa                                                                                                                    >  2:942709/1‑150 (MQ=255)
                                 ccGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCaaaa                                                                                                                  >  2:353410/1‑150 (MQ=255)
                                 ccGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGGCaaaa                                                                                                                  >  2:104315/1‑150 (MQ=255)
                                 ccGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCGTTTGTCGGGTCaaaa                                                                                                                  >  2:211353/1‑150 (MQ=255)
                                      tAAAACAACTTCGAACGGATATCCTTTTGAATGCGGTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAc                                                                                                             <  1:986841/150‑1 (MQ=255)
                                             acTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATc                                                                                                      <  2:1612703/150‑1 (MQ=255)
                                             acTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATc                                                                                                      >  2:756084/1‑150 (MQ=255)
                                              cTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGATTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCa                                                                                                     <  1:1502773/150‑1 (MQ=255)
                                              cTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACAAACAGACACATACCTGTTTTGTTGTTGTACATGATAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCa                                                                                                     >  2:809968/1‑150 (MQ=255)
                                               ttCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAg                                                                                                    <  2:1614304/150‑1 (MQ=255)
                                                    aCGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGc                                                                                               <  1:431166/150‑1 (MQ=255)
                                                     cGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGcc                                                                                              <  1:1083667/150‑1 (MQ=255)
                                                     cGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGcc                                                                                              <  1:1724008/150‑1 (MQ=255)
                                                     cGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGcc                                                                                              <  1:429037/150‑1 (MQ=255)
                                                       gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                            >  1:425997/1‑150 (MQ=255)
                                                        atatCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCat                                                                                           <  1:884369/150‑1 (MQ=255)
                                                             cTTTTGATTGCGTTGTACAAGGAACACACAGAAACATACCTGTTTTTTTGTTGTACATGAAAGGACTCAGGACAACAGCTGTACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCgg                                                                                      <  1:302759/150‑1 (MQ=255)
                                                              ttttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCggg                                                                                     >  1:1043102/1‑150 (MQ=255)
                                                               tttGATTGCGTTGTACAACGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGt                                                                                    <  1:700272/150‑1 (MQ=255)
                                                                ttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                   <  2:1805647/150‑1 (MQ=255)
                                                                ttGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTa                                                                                   <  1:1775319/150‑1 (MQ=255)
                                                                      gCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATc                                                                             >  1:1791457/1‑150 (MQ=255)
                                                                      gCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATc                                                                             >  1:232038/1‑150 (MQ=255)
                                                                      gCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATc                                                                             <  1:147653/150‑1 (MQ=255)
                                                                        gTTGTACAAGGAACACACAGACACATATCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCgg                                                                           >  1:453731/1‑150 (MQ=255)
                                                                            tACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCtta                                                                       <  2:1204829/150‑1 (MQ=255)
                                                                              cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCCAGCCAAATCAGATCGCCCAGATCGGGTACGTATCGGcatccc                                                                     >  2:1319467/1‑145 (MQ=255)
                                                                              cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                     >  2:1309047/1‑150 (MQ=255)
                                                                                aGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccat                                                                   >  2:211035/1‑150 (MQ=255)
                                                                                aGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccat                                                                   >  2:546573/1‑150 (MQ=255)
                                                                                 ggAACACACAGACACATACCTGTTTTGTTGTTGTACATGATAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCTGGTCAAAATCAACCCAAATCAGATCGCCCAAATCGGGTACGTATCGGCttaccatt                                                                  >  2:768415/1‑150 (MQ=255)
                                                                                       cacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGa                                                            >  2:68500/1‑150 (MQ=255)
                                                                                           gACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                        >  2:252430/1‑150 (MQ=255)
                                                                                            acacATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGGCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTcc                                                       <  2:1690554/150‑1 (MQ=255)
                                                                                             cacaTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCt                                                      <  2:1665440/150‑1 (MQ=255)
                                                                                              acaTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCtt                                                     <  1:723459/150‑1 (MQ=255)
                                                                                               caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTa                                                    >  2:770314/1‑150 (MQ=255)
                                                                                                 tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                  >  2:1131158/1‑150 (MQ=255)
                                                                                                   ccTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCtt                                                >  1:975238/1‑150 (MQ=255)
                                                                                                      gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGCCTTCCTTATCTTTCg                                             >  1:596288/1‑150 (MQ=255)
                                                                                                      gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCg                                             <  2:1108435/150‑1 (MQ=255)
                                                                                                        tttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGc                                           <  1:1536265/150‑1 (MQ=255)
                                                                                                            ttgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCtctc                                       >  1:679644/1‑150 (MQ=255)
                                                                                                            ttgttgtACATGAAAGGACTCAGGACAACAGATGGACCAAGTCAAGCTTGCTCGCTACCTTTTGTCGGGTCAAAAACAACCCAAATCAGATCGCCCATAACGGGTACGAATCGGCTTACCAATACCAGAATTACTTATCTTTAGGCtctc                                       >  1:781658/1‑150 (MQ=255)
                                                                                                                 gtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAg                                  <  1:279101/150‑1 (MQ=255)
                                                                                                                   aCATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                                <  2:1130680/150‑1 (MQ=255)
                                                                                                                    cATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCg                               <  1:599533/150‑1 (MQ=255)
                                                                                                                     aTGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGa                              <  1:1106097/150‑1 (MQ=255)
                                                                                                                       gAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTGACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGatt                            >  1:1382543/1‑149 (MQ=255)
                                                                                                                        aaaGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGatat                           >  2:617894/1‑150 (MQ=255)
                                                                                                                         aaGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATAtt                          >  2:428725/1‑150 (MQ=255)
                                                                                                                              cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCCACCCAAATCAGCACGCCCAAATCGGGTACGTATCGGCTTACCATTACCAGACTTCCATATCTTTCGGCTCACCCCAGTCGATATACTCGt                     >  1:372511/1‑150 (MQ=255)
                                                                                                                               tCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTg                    <  2:1536226/150‑1 (MQ=255)
                                                                                                                                cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                   >  1:1676471/1‑150 (MQ=255)
                                                                                                                                cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                   >  1:242777/1‑150 (MQ=255)
                                                                                                                                cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCAATATCGGGTTCGTATAGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCAAGTCGATATTCTCGTgg                   >  2:100843/1‑150 (MQ=255)
                                                                                                                                      gcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCAAATCGGGTACGTATCGGCTTACGATTCCCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGAGATTCTCGTGGAGGttt             >  1:777830/4‑150 (MQ=255)
                                                                                                                                      aacaGCTGGCCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCCAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttt             >  1:1738767/1‑150 (MQ=255)
                                                                                                                                      aacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttt             >  2:1501014/1‑150 (MQ=255)
                                                                                                                                       acaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCAATATTCTCGTGGAGGtttt            >  1:146293/1‑150 (MQ=255)
                                                                                                                                        caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGTTACGTATCGGCTTACCATTACCAGACTTCCTTATCTGTCGGCTCTCCCCAGTCGATATTCACGAGGAGGtttta           >  2:303741/1‑149 (MQ=255)
                                                                                                                                         aGCTGGACGATGTCCAGCTTGCGCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTcc          >  2:71425/1‑150 (MQ=255)
                                                                                                                                          gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg         >  2:1385149/1‑150 (MQ=255)
                                                                                                                                          gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg         <  2:715813/150‑1 (MQ=255)
                                                                                                                                           cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg        <  1:1616375/150‑1 (MQ=255)
                                                                                                                                             ggACGATGTCCATCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGGCGATATTCTCGTGGAGGTTTTCCGGCg      <  2:507971/150‑1 (MQ=255)
                                                                                                                                             ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg      >  1:1417420/1‑150 (MQ=255)
                                                                                                                                             ggAAGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATATTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg      >  1:1502279/1‑150 (MQ=255)
                                                                                                                                              gACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTCCCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt     >  1:708285/1‑150 (MQ=255)
                                                                                                                                               aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg    <  2:1423177/150‑1 (MQ=255)
                                                                                                                                               aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg    >  1:1641012/1‑150 (MQ=255)
                                                                                                                                               aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGGCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCGCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg    <  2:370106/150‑1 (MQ=255)
                                                                                                                                                 gATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGAt  >  1:1574423/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGAT  >  NC_000913/2910862‑2911156

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