Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,568,720 C→T 100% A33T (GCC→ACC)  eutP ← putative P‑loop NTPase ethanolamine utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,568,7200CT97.2% 855.6 / 9.7 249A33T (GCC→ACC) eutPputative P‑loop NTPase ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base T (98/143);  total (102/147)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 249 total aligned reads displayed.

GCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGCGTATAATTTAA  >  CP009273/2568573‑2568864
                                                                                                                                                   |                                                                                                                                                
gcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCt                                                                                                                                                <  2:3733451/150‑1 (MQ=255)
 caACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtg                                                                                                                                               <  1:1985343/150‑1 (MQ=255)
 caACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtg                                                                                                                                               <  2:167315/150‑1 (MQ=255)
 caACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtg                                                                                                                                               >  2:439895/1‑150 (MQ=255)
  aaCCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtgt                                                                                                                                              <  2:914835/150‑1 (MQ=255)
  aaCCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtgt                                                                                                                                              >  2:3925999/1‑150 (MQ=255)
     cGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGtt                                                                                                                                           <  1:4622618/150‑1 (MQ=255)
     cGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGtt                                                                                                                                           >  2:60096/1‑150 (MQ=255)
     cGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGtt                                                                                                                                           <  1:2632040/150‑1 (MQ=255)
       aCTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGtttt                                                                                                                                         <  2:4236097/150‑1 (MQ=255)
       aCTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTACACGGTCTGTGtttt                                                                                                                                         <  1:2983552/150‑1 (MQ=255)
         tttCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTc                                                                                                                                       >  2:4637218/1‑150 (MQ=255)
         tttCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTAAATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTc                                                                                                                                       <  1:2844389/150‑1 (MQ=255)
          ttCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCt                                                                                                                                      <  1:4638896/150‑1 (MQ=255)
           tCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTg                                                                                                                                     <  1:1478074/150‑1 (MQ=255)
            ccGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTgg                                                                                                                                    <  1:550748/150‑1 (MQ=255)
            ccGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTgg                                                                                                                                    >  1:1643999/1‑150 (MQ=255)
              gggTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCg                                                                                                                                  >  2:1029517/1‑150 (MQ=255)
                     ttGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGtat                                                                                                                           <  1:4389307/150‑1 (MQ=255)
                        cgcCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAAt                                                                                                                        <  1:1896579/150‑1 (MQ=255)
                         gcCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAAtt                                                                                                                       >  2:2344020/1‑150 (MQ=255)
                           cATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTcc                                                                                                                     <  2:2661501/150‑1 (MQ=255)
                           cATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTcc                                                                                                                     >  2:1118975/1‑150 (MQ=255)
                            aTGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTccc                                                                                                                    >  2:1289192/1‑150 (MQ=255)
                              ggACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCt                                                                                                                  >  2:4584189/1‑150 (MQ=255)
                                aCATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGc                                                                                                                <  1:16003/150‑1 (MQ=255)
                                  aTAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCaa                                                                                                              <  1:1772013/150‑1 (MQ=255)
                                        cAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCAtt                                                                                                        >  2:3236078/1‑150 (MQ=255)
                                         aGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTa                                                                                                       >  1:1283516/1‑150 (MQ=255)
                                            aTGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAAt                                                                                                    >  2:39534/1‑150 (MQ=255)
                                             tGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATa                                                                                                   <  2:1573498/150‑1 (MQ=255)
                                             tGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATa                                                                                                   <  1:2028783/150‑1 (MQ=255)
                                                cTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACg                                                                                                <  1:25340/150‑1 (MQ=255)
                                                cTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACg                                                                                                <  2:1769312/150‑1 (MQ=255)
                                                  aCATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgtt                                                                                              >  2:3598143/1‑150 (MQ=255)
                                                  aCATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgtt                                                                                              >  1:2456921/1‑150 (MQ=255)
                                                  aCATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgtt                                                                                              <  2:4389952/150‑1 (MQ=255)
                                                   cATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgttg                                                                                             <  1:793600/150‑1 (MQ=255)
                                                   cATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgttg                                                                                             <  1:4211382/150‑1 (MQ=255)
                                                      ccTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGttt                                                                                          <  2:768629/150‑1 (MQ=255)
                                                      ccTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGttt                                                                                          >  2:2191256/1‑150 (MQ=255)
                                                      ccTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGttt                                                                                          <  1:3169022/150‑1 (MQ=255)
                                                       cTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGtttt                                                                                         <  2:2904471/150‑1 (MQ=255)
                                                         gTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTc                                                                                       >  2:652420/1‑150 (MQ=255)
                                                          tAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTcc                                                                                      >  2:3099962/1‑150 (MQ=255)
                                                           aGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCt                                                                                     <  1:4366763/150‑1 (MQ=255)
                                                           aGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCt                                                                                     <  2:801025/150‑1 (MQ=255)
                                                            gCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTg                                                                                    <  1:4522951/150‑1 (MQ=255)
                                                               tGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGccc                                                                                 <  1:1839962/150‑1 (MQ=255)
                                                                ggTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGcccc                                                                                <  2:4081628/150‑1 (MQ=255)
                                                                ggTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGcccc                                                                                <  2:1579660/150‑1 (MQ=255)
                                                                    aTTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACg                                                                            <  1:1194223/150‑1 (MQ=255)
                                                                     ttAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGc                                                                           >  1:2092095/1‑150 (MQ=255)
                                                                      tAAGGCGTGATACCAGCGGGGACGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCt                                                                          <  1:578313/150‑1 (MQ=255)
                                                                         ggCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTAcc                                                                       <  1:1363489/150‑1 (MQ=255)
                                                                           cGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCca                                                                     >  1:4496527/1‑150 (MQ=255)
                                                                                 aCCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACg                                                               >  2:4190914/1‑150 (MQ=255)
                                                                                  ccAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGc                                                              <  1:1727267/150‑1 (MQ=255)
                                                                                       ggggTTGGTTAAAATATTCACCCGGCGTATCATTATCGCCCTTATCATTATATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTTAATAACGTTGTTTTTCCTGCCCCGACGCTACCAACAAACGCAATTc                                                         <  2:1960641/150‑1 (MQ=255)
                                                                                       ggggATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTc                                                         <  1:1530285/150‑1 (MQ=255)
                                                                                       ggggATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTc                                                         <  2:3903297/150‑1 (MQ=255)
                                                                                            tGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTc                                                    <  2:4144692/150‑1 (MQ=255)
                                                                                              gTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCAt                                                  <  2:2793917/150‑1 (MQ=255)
                                                                                                tAAAATATTCATCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGa                                                >  1:4563344/1‑150 (MQ=255)
                                                                                                 aaaaTATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGaa                                               <  1:2244729/150‑1 (MQ=255)
                                                                                                  aaaTATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAAt                                              >  2:4009916/1‑150 (MQ=255)
                                                                                                  aaaTATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAAt                                              <  1:745369/150‑1 (MQ=255)
                                                                                                     tatTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGcc                                           <  1:3248438/150‑1 (MQ=255)
                                                                                                      atTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCt                                          >  2:1371270/1‑150 (MQ=255)
                                                                                                       ttaaCCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTc                                         <  1:1972352/147‑1 (MQ=255)
                                                                                                       ttCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTc                                         <  1:3210550/150‑1 (MQ=255)
                                                                                                        tCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCg                                        >  1:1419244/1‑150 (MQ=255)
                                                                                                         cACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGt                                       >  2:4475895/1‑150 (MQ=255)
                                                                                                             cGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTa                                   <  2:1466324/150‑1 (MQ=255)
                                                                                                             cGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTa                                   <  1:4186542/150‑1 (MQ=255)
                                                                                                                cGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACt                                <  2:2231089/150‑1 (MQ=255)
                                                                                                                 gTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACtt                               <  1:2330875/150‑1 (MQ=255)
                                                                                                                 gTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACtt                               >  2:3049052/1‑150 (MQ=255)
                                                                                                                   aTCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACtttt                             <  2:4043265/150‑1 (MQ=255)
                                                                                                                    tCAATATCTACCTTATCATTAAATTACACGTTCTGTGTTTTTCTGGCGAGGGTATAATTACCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACTCTAACCACAAACGCAATTCGTTTCATGATTGCCTCGTATTAACttttt                            <  1:1730506/150‑2 (MQ=255)
                                                                                                                      aaTATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGt                          <  1:3675572/150‑1 (MQ=255)
                                                                                                                      aaTATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGt                          <  1:1006059/150‑1 (MQ=255)
                                                                                                                      aaTATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGt                          >  1:2663254/1‑150 (MQ=255)
                                                                                                                         atCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCAt                       >  2:3993874/1‑150 (MQ=255)
                                                                                                                          tCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCAtt                      <  1:3712415/150‑1 (MQ=255)
                                                                                                                           cGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCAttt                     <  1:643328/150‑1 (MQ=255)
                                                                                                                             cccTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTca                   >  1:2357678/1‑150 (MQ=255)
                                                                                                                              ccTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTcac                  >  2:2933260/1‑150 (MQ=255)
                                                                                                                               cTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTcaca                 >  2:2042187/1‑150 (MQ=255)
                                                                                                                                 tATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACaaa               <  1:3495719/150‑1 (MQ=255)
                                                                                                                                  aTCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAg              <  2:1872166/150‑1 (MQ=255)
                                                                                                                                  aTCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAg              <  2:1126093/150‑1 (MQ=255)
                                                                                                                                  aTCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAg              <  2:232394/150‑1 (MQ=255)
                                                                                                                                   tCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGc             <  1:3102663/150‑1 (MQ=255)
                                                                                                                                   tCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGc             <  1:945853/150‑1 (MQ=255)
                                                                                                                                   tCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGc             >  2:2048498/1‑150 (MQ=255)
                                                                                                                                   tCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGc             >  1:1083170/1‑150 (MQ=255)
                                                                                                                                    cATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGCg            >  2:3510613/1‑150 (MQ=255)
                                                                                                                                       tAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGCGtat         <  1:1905732/150‑1 (MQ=255)
                                                                                                                                          aTTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGCGTATAAt      >  1:1417714/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTCCTGCCCCGACGCTACCCACAAACGCAATTCGTTTCATGAATGCCTCGTATTAACTTTTGGTCATTTCACAAAGCGTATAATTTAA  >  CP009273/2568573‑2568864

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