Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,071,745 A→G P7P (CCT→CCC rutD ← putative reactive intermediate detoxifying aminoacrylate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,071,7450AG98.1% 388.1 / ‑0.5 103coding (21/801 nt)rutDputative reactive intermediate detoxifying aminoacrylate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (46/55);  minor base . (1/1);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGT  >  NC_000913/1071600‑1071891
                                                                                                                                                 |                                                                                                                                                  
agtCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTg                                                                                                                                                >  4:1357512/2‑150 (MQ=255)
 gtCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTgag                                                                                                                                              <  4:291707/151‑1 (MQ=255)
    gggATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTCCCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTgaga                                                                                                                                             <  2:543705/149‑1 (MQ=255)
    gggATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGagag                                                                                                                                            <  4:1443837/150‑1 (MQ=255)
    gggATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGagag                                                                                                                                            >  4:222215/1‑150 (MQ=255)
    gggATTATTGCCGGTGCCGCGCTGGTCGAAACAGACTCCCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGagag                                                                                                                                            <  1:319906/150‑1 (MQ=255)
       attattGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGa                                                                                                                                         >  2:784583/1‑150 (MQ=255)
        ttattGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCCCCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGa                                                                                                                                         <  4:1300661/149‑1 (MQ=255)
          attGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGt                                                                                                                                     >  3:164476/1‑151 (MQ=255)
          attGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGGGAGAGTGaaa                                                                                                                                       >  2:234793/1‑149 (MQ=255)
           ttGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGtt                                                                                                                                    >  2:982920/1‑151 (MQ=255)
           ttGCCGGTGCCGCGCGGGTCGTAACAGACTCCCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGt                                                                                                                                     <  4:707058/150‑1 (MQ=255)
            tGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGttt                                                                                                                                   >  4:320300/1‑151 (MQ=255)
            tGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGtt                                                                                                                                    <  4:301654/150‑1 (MQ=255)
            tGCCGGTGCCACGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGttt                                                                                                                                   >  4:76774/1‑151 (MQ=255)
             gccgttccggcgcTGGTGTTAACAAACCCCCAGATACCCCTGCTCCAGCTCCGCCAGTGGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGt                                                                                                                                     <  2:1037135/139‑1 (MQ=255)
                ggTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATc                                                                                                                               <  4:1080009/151‑1 (MQ=255)
                    ccGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATcgcg                                                                                                                            <  4:416484/150‑1 (MQ=255)
                      gcgcTGGTCGTAACAGACTCCCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGTGAGAGTGAAAGTTTCATCGCGGcc                                                                                                                         <  4:1058873/150‑1 (MQ=255)
                      gcgcTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGcc                                                                                                                         >  2:1577885/1‑151 (MQ=255)
                      gcgcTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGcc                                                                                                                         <  1:98005/151‑1 (MQ=255)
                       cgcTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCt                                                                                                                        >  4:545094/1‑151 (MQ=255)
                         cTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCt                                                                                                                        >  4:37317/1‑149 (MQ=255)
                           ggTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCAc                                                                                                                     <  4:1644894/150‑1 (MQ=255)
                           ggTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACt                                                                                                                    >  2:811029/1‑151 (MQ=255)
                            gTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACtt                                                                                                                   <  4:785617/151‑1 (MQ=255)
                              cGTAACAGCCTGCCTGATCCTCCTCCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTg                                                                                                                  <  2:713387/150‑1 (MQ=255)
                              cGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTgg                                                                                                                 >  2:1553446/1‑151 (MQ=255)
                               gTAACAGTCTCCCTGTTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTACCTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTgg                                                                                                                 <  4:368575/150‑1 (MQ=255)
                                tAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCg                                                                                                               >  2:215427/1‑151 (MQ=255)
                                 aaCAGACTCCCTGATACCCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCg                                                                                                               <  2:86913/150‑1 (MQ=255)
                                 aaCAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCg                                                                                                               <  2:769941/150‑1 (MQ=255)
                                     gACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATg                                                                                                          <  4:669463/151‑1 (MQ=255)
                                      aCTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATgc                                                                                                         >  4:1250349/1‑151 (MQ=255)
                                          ccTGATCCTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCa                                                                                                      <  2:128189/150‑1 (MQ=255)
                                          ccTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATgcg                                                                                                        <  4:767242/148‑1 (MQ=255)
                                          ccTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCaa                                                                                                     <  2:414070/151‑1 (MQ=255)
                                             gATCCTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAAt                                                                                                    <  1:121572/149‑1 (MQ=255)
                                              aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATtgt                                                                                                 >  2:38953/1‑151 (MQ=255)
                                              aTACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATtgt                                                                                                 <  4:395154/151‑1 (MQ=255)
                                               tACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATtgt                                                                                                 <  2:997705/150‑1 (MQ=255)
                                                   ccTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGc                                                                                              <  1:399372/149‑1 (MQ=255)
                                                   ccTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCg                                                                                             <  2:249627/150‑1 (MQ=255)
                                                     tGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGc                                                                                              <  4:1036245/147‑1 (MQ=255)
                                                         ccAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGGCCCGAAATCAACCCCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTcc                                                                                       <  2:484321/150‑1 (MQ=255)
                                                         ccAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTcc                                                                                       >  1:44750/1‑150 (MQ=255)
                                                         ccAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTc                                                                                        <  4:553929/149‑1 (MQ=255)
                                                          cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTccac                                                                                     <  3:473793/151‑1 (MQ=255)
                                                          cAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTccac                                                                                     >  2:1204773/1‑151 (MQ=255)
                                                           aGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTccacc                                                                                    >  2:1578113/1‑151 (MQ=255)
                                                             cACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAg                                                                                  <  2:1231851/151‑1 (MQ=255)
                                                              aCCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAgc                                                                                 >  4:1460165/1‑151 (MQ=255)
                                                                     gTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCa                                                                           <  4:1132847/150‑1 (MQ=255)
                                                                           gTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTa                                                                    <  4:1349981/151‑1 (MQ=255)
                                                                             aaCCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTAcc                                                                  >  2:134510/1‑151 (MQ=255)
                                                                             aaCCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTAc                                                                   <  2:67258/150‑1 (MQ=255)
                                                                               ccAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCa                                                                 <  2:1451467/150‑1 (MQ=255)
                                                                                cAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGt                                                               <  4:367214/151‑1 (MQ=255)
                                                                                   tAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCg                                                            >  2:610208/1‑151 (MQ=255)
                                                                                    aaCTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCg                                                            <  4:1481630/150‑1 (MQ=255)
                                                                                      cTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCa                                                          <  2:698535/150‑1 (MQ=255)
                                                                                      cTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCAc                                                         >  3:20112/1‑151 (MQ=255)
                                                                                           gCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAAt                                                    >  4:88484/1‑151 (MQ=255)
                                                                                           gCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAAt                                                    >  4:494253/1‑151 (MQ=255)
                                                                                            cTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATg                                                   >  2:669688/1‑151 (MQ=255)
                                                                                            cTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATg                                                   >  3:125727/1‑151 (MQ=255)
                                                                                              aCCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGc                                                  >  4:420221/1‑150 (MQ=255)
                                                                                               cccccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAg                                                <  4:1183753/151‑1 (MQ=255)
                                                                                                ccccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGa                                               <  4:701929/151‑1 (MQ=255)
                                                                                                ccccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGa                                               <  3:37704/151‑1 (MQ=255)
                                                                                                  ccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAAt                                             <  4:780334/151‑1 (MQ=255)
                                                                                                  ccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAAt                                             <  4:737644/151‑1 (MQ=255)
                                                                                                              aTCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTa                                 <  3:159516/151‑1 (MQ=255)
                                                                                                                 aaCACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATc                               >  3:408399/1‑150 (MQ=255)
                                                                                                                  acacCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCAc                             >  4:33720/1‑151 (MQ=255)
                                                                                                                   caccacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCAc                             <  2:535116/150‑1 (MQ=255)
                                                                                                                    accacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCAccc                           <  4:1060538/151‑1 (MQ=255)
                                                                                                                    accacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCAccc                           >  4:874091/1‑151 (MQ=255)
                                                                                                                     ccacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCg                          >  4:7649/1‑151 (MQ=255)
                                                                                                                     ccacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCg                          >  4:53688/1‑151 (MQ=255)
                                                                                                                      cacTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCgg                         <  3:131276/151‑1 (MQ=255)
                                                                                                                        cTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGa                        >  4:107091/1‑150 (MQ=255)
                                                                                                                          aCGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGaaaa                     <  3:444621/151‑1 (MQ=255)
                                                                                                                          aCGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGaa                       <  4:216386/149‑1 (MQ=255)
                                                                                                                          aCGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGaa                       <  4:816193/149‑1 (MQ=255)
                                                                                                                            gggCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGaaaaaa                   <  1:218811/151‑1 (MQ=255)
                                                                                                                            gggCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGaaaaa                    >  4:1381220/1‑150 (MQ=255)
                                                                                                                               cgcATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTc                >  3:218946/1‑151 (MQ=255)
                                                                                                                               cgcATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTc                >  2:203632/1‑151 (MQ=255)
                                                                                                                               cgcATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTc                <  4:812104/151‑1 (MQ=255)
                                                                                                                               cgcATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTc                <  4:1495428/151‑1 (MQ=255)
                                                                                                                                gcATCAGCATAAGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCgg              >  4:839854/1‑151 (MQ=255)
                                                                                                                                 cATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCgg              <  4:1481218/151‑1 (MQ=255)
                                                                                                                                     aGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGt           >  2:495280/1‑150 (MQ=255)
                                                                                                                                     aGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGt           >  4:831528/1‑150 (MQ=255)
                                                                                                                                       cATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGt           >  3:471878/1‑148 (MQ=255)
                                                                                                                                         tAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGAtt      >  2:559226/1‑151 (MQ=255)
                                                                                                                                           aGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTc    >  4:285535/1‑151 (MQ=255)
                                                                                                                                             ggggggGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGt  <  4:495974/151‑1 (MQ=255)
                                                                                                                                             ggggggGCGAGAGCGAAAGTTTTATCGCGCACTCACTGCCCGGTATGCAAACTTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCGGCACTCAATGCAGAAGCGCGACGGTTCACCAACAGGACAAACACCGTAGTAGAtttag   >  1:412902/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGTCGGGATTATTGCCGGTGCCGCGCTGGTCGTAACAGACTACCTGATACTCCTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGT  >  NC_000913/1071600‑1071891

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