Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,009 T→C 5.6% 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 5.6% 378.1 / 9.3 125H28R (CAT→CGT) ‡mazFmRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base T (59/58);  new base C (2/5);  total (62/63)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 125 total aligned reads displayed.

CTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATG  >  NC_000913/2910861‑2911157
                                                                                                                                                    |                                                                                                                                                    
cTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTATGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATg                                                                                                                                                     <  1:1560905/150‑1 (MQ=255)
  gATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTc                                                                                                                                                   >  1:401475/1‑150 (MQ=255)
   aTCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTcc                                                                                                                                                  <  1:1410782/150‑1 (MQ=255)
    tCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGCACTCAGGACAACAGCTGGACGATGTCCa                                                                                                                                                 <  1:1844454/150‑1 (MQ=255)
       gCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCt                                                                                                                                              <  1:695265/150‑1 (MQ=255)
        cTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCtt                                                                                                                                             <  2:758554/150‑1 (MQ=255)
         tAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                            <  2:864678/150‑3 (MQ=255)
         aaaCGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgt                                                                                                                                            <  2:1013475/149‑3 (MQ=255)
            cGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTc                                                                                                                                         >  2:1423421/1‑150 (MQ=255)
             gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACATCGTCCAGCtgttgtc                                                                                                                                        <  2:1632116/150‑7 (MQ=255)
               tACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCt                                                                                                                                      <  1:447293/150‑1 (MQ=255)
               aaCGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCt                                                                                                                                      <  1:476676/149‑1 (MQ=255)
                 cGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTAc                                                                                                                                    >  1:645886/1‑150 (MQ=255)
                     aTCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCttt                                                                                                                                <  2:664334/150‑1 (MQ=255)
                        aCGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGt                                                                                                                             <  1:1032537/150‑1 (MQ=255)
                        aCGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGt                                                                                                                             <  1:798805/150‑1 (MQ=255)
                         cGTTACTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTTCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTc                                                                                                                            <  2:221432/150‑1 (MQ=255)
                           ttCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCgg                                                                                                                          <  2:1414019/150‑1 (MQ=255)
                            tCCTGACCGGATAAAACTACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGTAACACACAGACACATACCTGTTTTGTGGTTGTACATGAAAGGACTCAGGACAACAGCTTGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                         <  1:288046/150‑1 (MQ=255)
                             ccTGACCGGATAAAACAACTTCGCACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGCCGATGTCCAGCTTGCTCGCTACCTTTTGACGGGt                                                                                                                        >  2:1339299/1‑150 (MQ=255)
                             ccTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGt                                                                                                                        <  2:1705751/150‑1 (MQ=255)
                             ccTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGt                                                                                                                        <  1:1298716/150‑1 (MQ=255)
                                gACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCaa                                                                                                                     <  2:1286505/150‑1 (MQ=255)
                                gACCGCATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGATTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGAAGTCCAGCTTGCTCGCTACCTTTTGTCGTGACaa                                                                                                                     >  2:1771608/1‑150 (MQ=255)
                                   cGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACAGGAAATGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAAt                                                                                                                  <  1:210081/150‑1 (MQ=255)
                                        aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                             <  1:419783/150‑1 (MQ=255)
                                        aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                             <  2:401475/150‑1 (MQ=255)
                                         aaaCAACTTCGAACGGATATACTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTTTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTACAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcac                                                                                                            >  1:664164/1‑148 (MQ=255)
                                           acaacTTCGAACGTATCTAATTTTGATTGCGTTGTACAAGGATCACACAGACACATACCTGTTTTGTTGTTGTACATGAAATGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTGTCGGGTCAAAATCAACCCaaa                                                                                                         <  1:271237/150‑1 (MQ=255)
                                           acaacTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCaa                                                                                                          >  2:1583796/1‑150 (MQ=255)
                                              acTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATc                                                                                                       <  1:427216/150‑1 (MQ=255)
                                              aattCGAACTGATATCCTTTTGATTGAGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATc                                                                                                       <  1:1359449/148‑1 (MQ=255)
                                                ttCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGGAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAg                                                                                                     <  2:781345/150‑1 (MQ=255)
                                                   taagggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATc                                                                                                  <  2:1726606/146‑1 (MQ=255)
                                                   gcaCGGATATCCTTTTGATTGCGTTGTACAAGGATCACACAGACACATACCTGTTTGGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATc                                                                                                  <  2:110055/148‑1 (MQ=255)
                                                   gAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATc                                                                                                  <  2:1550912/150‑1 (MQ=255)
                                                   gAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATc                                                                                                  >  2:94516/1‑150 (MQ=255)
                                                    aaCGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCg                                                                                                 >  2:176347/1‑150 (MQ=255)
                                                       ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                              <  2:356570/150‑1 (MQ=255)
                                                       ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                              <  2:473337/150‑1 (MQ=255)
                                                        gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCTGATCGCCCa                                                                                             <  1:1025019/150‑1 (MQ=255)
                                                                   gATTGCGTTGTACAAGGAACACACAGACACATACCTGTATTGGTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACg                                                                                  <  1:1094201/150‑1 (MQ=255)
                                                                    aTTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGt                                                                                 >  1:898306/1‑150 (MQ=255)
                                                                          ttGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGc                                                                           <  1:1423421/150‑1 (MQ=255)
                                                                          ttGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGc                                                                           <  1:1874096/150‑1 (MQ=255)
                                                                            gTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCtt                                                                         >  2:68465/1‑150 (MQ=255)
                                                                               cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGTCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                      <  2:79331/150‑1 (MQ=255)
                                                                                      cacacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacca                                                               >  1:1774948/1‑150 (MQ=255)
                                                                                       acacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAg                                                              <  1:992085/150‑1 (MQ=255)
                                                                                         acaGACACATACCTGTTTTGTTGTTGTACATGAAACGACTCAGAACAACAGCTGGAAGATGAACAGCTTGCTAGCTACCATTAGTCTATTCAAAATCAACCCAAATCACATCGCCCAAATCGAGTACTTAACGGCATACCATTACCAGAc                                                            >  1:253421/1‑150 (MQ=255)
                                                                                          caGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACTTATCGGCTTACCAGTACCAGACt                                                           >  1:1032764/1‑150 (MQ=255)
                                                                                            gACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGTACAACAGCTGGACGATGTCCAGCTTGCGATCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCAAATCGGGTACGTATCCGCTTACCATTAACAGACTTc                                                         >  2:1543115/1‑150 (MQ=255)
                                                                                            gACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                         <  1:1856997/150‑1 (MQ=255)
                                                                                                caTACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTa                                                     >  1:1412916/1‑150 (MQ=255)
                                                                                                  tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                   >  2:1128817/1‑150 (MQ=255)
                                                                                                  tACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATc                                                   >  2:1619682/1‑150 (MQ=255)
                                                                                                    ccTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACGATTACCAGACTTCCTTATCtt                                                 >  2:795621/1‑150 (MQ=255)
                                                                                                    ccTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTAGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCACATCGCCCATATCGGGAACGTATCGGCTTACCATAACCACACTTCCTTATCtt                                                 >  2:852736/1‑150 (MQ=255)
                                                                                                     cTGATTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATCTCCAGCATGCACGCTACCTTTTGTCGGGACAAAATCACCCCAAATCACATCCCCCAAAGCGGCTACGTAACGGCATCCCATTCCCAGCCGACCATATCtta                                                >  1:142233/1‑149 (MQ=255)
                                                                                                       gTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCg                                              >  2:1098794/1‑150 (MQ=255)
                                                                                                        ttttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCgg                                             <  1:1910345/150‑1 (MQ=255)
                                                                                                        ttttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACAGACTTCCTTATCTTTCGGc                                            <  2:1058570/150‑1 (MQ=255)
                                                                                                         tttGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGc                                            <  1:94516/150‑1 (MQ=255)
                                                                                                              tgttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcc                                       <  1:1013487/150‑1 (MQ=255)
                                                                                                               gttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCACCCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTccc                                      >  1:1176954/1‑150 (MQ=255)
                                                                                                                  gtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCCACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAg                                   >  1:1462679/1‑150 (MQ=255)
                                                                                                                  gtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAg                                   >  2:684831/1‑150 (MQ=255)
                                                                                                                   tACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                                 <  2:784887/150‑1 (MQ=255)
                                                                                                                    aCATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTc                                 >  2:1527561/1‑150 (MQ=255)
                                                                                                                     cATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCg                                >  1:462126/1‑150 (MQ=255)
                                                                                                                      aTGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGa                               <  2:1185287/150‑1 (MQ=255)
                                                                                                                       tGAAAGGCCTCAGGACAACAGCTGGACGATGTCCAGCGTGCTCGCTCCCTTTTGTCGGGGCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTTCCAGACTTCCTTATCTTTCGGCTCTCCCCAGGCGat                              >  2:258702/1‑150 (MQ=255)
                                                                                                                       tGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                              >  1:414918/1‑150 (MQ=255)
                                                                                                                         aaaGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGatat                            <  1:1751964/150‑1 (MQ=255)
                                                                                                                            ggACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATtct                         <  2:1867483/150‑1 (MQ=255)
                                                                                                                            ggACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATAACGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATtct                         >  2:1299473/1‑150 (MQ=255)
                                                                                                                             gACTCAGGACAACACCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACCTTACGGCTGACCAATACCAGACTTCCTATTCGTTCGGCTCTCCCCCGTCGCTATactc                        >  2:1635807/1‑150 (MQ=255)
                                                                                                                              aCTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCg                       <  1:186715/150‑1 (MQ=255)
                                                                                                                              aCTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCg                       <  2:701739/150‑1 (MQ=255)
                                                                                                                               cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGt                      >  2:330247/1‑150 (MQ=255)
                                                                                                                               cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGGACGTATCGGCTTACCATTACCAGACTTCCTGATCTTTCGGCTCTCCCCAGTCGATTTTCTCgg                      >  1:680140/1‑149 (MQ=255)
                                                                                                                               cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTAAAAATCAACCCAAATCATATCGCCAATATCGGGTACGTATCGGCTTACCATTACCAGACATCCTTATCTTTAGTCTCTCCACATTAGATATTATCgg                      >  2:1871210/1‑149 (MQ=255)
                                                                                                                                 cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                    >  1:798629/1‑150 (MQ=255)
                                                                                                                                 cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                    >  2:1622241/1‑150 (MQ=255)
                                                                                                                                   ggACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAg                  <  2:847155/150‑1 (MQ=255)
                                                                                                                                   ggACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCGTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAg                  >  1:1815942/1‑150 (MQ=255)
                                                                                                                                    gACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAgg                 <  1:1415751/150‑1 (MQ=255)
                                                                                                                                    gACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAgg                 >  2:1061500/1‑150 (MQ=255)
                                                                                                                                     acaacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCGCCCCAGTCGATATTCGCGTGGAggg                >  1:1847844/1‑149 (MQ=255)
                                                                                                                                     acaacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATGACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGt                >  1:1664619/1‑150 (MQ=255)
                                                                                                                                      caacaGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATTTTCTCGTGGAGGtt               >  2:1094666/1‑150 (MQ=255)
                                                                                                                                      agcaaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTAGTCCGGACAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGCCTTACCATTACCAGACTACCTTATCTTTCGGCTCACCCCCGTCGAGAATCTCGTGGAGGtt               >  1:667599/5‑150 (MQ=255)
                                                                                                                                        caaGCTGGACATCGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt             >  2:1829737/3‑150 (MQ=255)
                                                                                                                                         caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGttttg            >  2:1846493/1‑149 (MQ=255)
                                                                                                                                            cTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCgg         <  1:1700645/150‑1 (MQ=255)
                                                                                                                                              ggACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCATTATCTTTCGGCTCTCCCCAGTAGATATTCTCGTGGAGGTTTTCCGGCg       >  2:140455/1‑150 (MQ=255)
                                                                                                                                              ggACGATGTCCAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg       <  2:129889/150‑1 (MQ=255)
                                                                                                                                                aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGCCTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg     >  1:1147272/1‑150 (MQ=255)
                                                                                                                                                 cGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGGTTTCCGGCGTGa    >  1:1358067/1‑150 (MQ=255)
                                                                                                                                                   aTGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTTCCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTGTTGGAGGTTTTCCGGCGTGATg  >  2:41757/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATG  >  NC_000913/2910861‑2911157

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