Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,007 G→A 5.7% R29C (CGT→TGT)  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,0070GA 5.7% 385.4 / 7.5 122R29C (CGT→TGT) mazFmRNA interferase toxin, antitoxin is MazE
Reads supporting (aligned to +/- strand):  ref base G (57/58);  new base A (2/5);  total (59/63)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 123 total aligned reads displayed.

TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGA  >  NC_000913/2910858‑2911155
                                                                                                                                                     |                                                                                                                                                    
   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)
                gCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACACGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCg                                                                                                                                      >  2:1086339/1‑150 (MQ=255)
                  tACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCt                                                                                                                                    <  1:447293/150‑1 (MQ=255)
                  aaCGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCt                                                                                                                                    <  1:476676/149‑1 (MQ=255)
                    cGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTAc                                                                                                                                  >  1:645886/1‑150 (MQ=255)
                     gCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTAcc                                                                                                                                 <  2:681175/150‑1 (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)
                               tCCTGACCGGATAAAACAACTTCGAACGGAAATCCTTTTTATTGCTTTGTACAAGTACCACACAGACACATACCTGTTTTGTTGTTGTACATCAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCggg                                                                                                                       <  2:1609544/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)
                                   gACCGCATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGATTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGAAGTCCAGCTTGCTCGCTACCTTTTGTCGTGACaa                                                                                                                   >  2:1771608/1‑150 (MQ=255)
                                      cGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACAGGAAATGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAAt                                                                                                                <  1:210081/150‑1 (MQ=255)
                                           aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                           <  2:401475/150‑1 (MQ=255)
                                           aaaaCAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAAcc                                                                                                           <  1:419783/150‑1 (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:94516/1‑150 (MQ=255)
                                                          ggttatCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCCTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                            <  2:360411/147‑1 (MQ=255)
                                                          ggATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGccc                                                                                            <  2:473337/150‑1 (MQ=255)
                                                           gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCTGATCGCCCa                                                                                           <  1:1025019/150‑1 (MQ=255)
                                                           gATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCa                                                                                           <  1:663986/150‑1 (MQ=255)
                                                              atCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCatat                                                                                        >  1:962469/1‑150 (MQ=255)
                                                                ccTTTTGATTGCGTTATACAAGGAACACACACACACATACCTGTATTGTTGTTGTACATGAAAGGACTCAGGACAACAGCAGGACGATGTCCAGCTTGCTCCCTACCTATTGTCGGCTCAAAATCAACCCAAATCAGATCGCCCAAATCg                                                                                      >  1:1040936/1‑150 (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)
                                                                                  cAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGGCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttacc                                                                    <  2:1162134/150‑1 (MQ=255)
                                                                                         cacacaGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacca                                                             >  1:1774948/1‑150 (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)
                                                                                                        cTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCttt                                              >  1:701739/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)
                                                                                                                   ttgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAAACAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTcccc                                   <  2:1324494/150‑1 (MQ=255)
                                                                                                                    tgtACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTTCCAGACTTCCTTATCTTTCGGCTCTCCCCa                                  >  2:1774700/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:198222/150‑1 (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                     <  2:701739/150‑1 (MQ=255)
                                                                                                                                 aCTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCg                     <  1:186715/150‑1 (MQ=255)
                                                                                                                                  cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGt                    >  2:330247/1‑150 (MQ=255)
                                                                                                                                  cTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGGACGTATCGGCTTACCATTACCAGACTTCCTGATCTTTCGGCTCTCCCCAGTCGATTTTCTCgg                    >  1:680140/1‑149 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                  >  2:1622241/1‑150 (MQ=255)
                                                                                                                                    cAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTgg                  >  2:1452236/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)
                                                                                                                                        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)
                                                                                                                                            caGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTc          >  1:1308068/1‑150 (MQ=255)
                                                                                                                                             aGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTcc         <  1:1849010/150‑1 (MQ=255)
                                                                                                                                              gCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCg        >  2:1676098/1‑150 (MQ=255)
                                                                                                                                                 ggACGATGTCCAGCTTGCTAGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCg     <  2:129889/150‑1 (MQ=255)
                                                                                                                                                   aCGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGCCTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTg   >  1:1147272/1‑150 (MQ=255)
                                                                                                                                                    cGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTTTCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGGTTTCCGGCGTGa  >  1:1414776/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGA  >  NC_000913/2910858‑2911155

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