Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,144,697 C→A 10.2% A27S (GCA→TCA)  argE ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,144,6970CA10.2% 1162.6 / 71.9 328A27S (GCA→TCA) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (138/156);  new base A (20/14);  total (158/170)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

ACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTT  >  CP009273/4144550‑4144846
                                                                                                                                                   |                                                                                                                                                     
acaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCct                                                                                                                                                     <  1:3960846/150‑1 (MQ=255)
acaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCct                                                                                                                                                     >  1:4157805/1‑150 (MQ=255)
acaaGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCct                                                                                                                                                     >  1:344382/1‑150 (MQ=255)
   aGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCtt                                                                                                                                                  >  2:3867138/1‑150 (MQ=255)
     ccgccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTcc                                                                                                                                                <  2:284005/150‑1 (MQ=255)
        ccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTg                                                                                                                                             >  2:3756945/1‑150 (MQ=255)
          aGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGgc                                                                                                                                           >  1:5109333/1‑150 (MQ=255)
           gCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGgcg                                                                                                                                          >  2:4333328/1‑150 (MQ=255)
           gCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGgcg                                                                                                                                          <  1:3116723/150‑1 (MQ=255)
             ccccTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCGGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCt                                                                                                                                        >  2:2327740/1‑150 (MQ=255)
                   tCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGa                                                                                                                                  <  2:1381883/150‑1 (MQ=255)
                   tCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGa                                                                                                                                  <  2:2046380/150‑1 (MQ=255)
                    ccGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGaa                                                                                                                                 <  2:3153097/150‑1 (MQ=255)
                      gATACGTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAgg                                                                                                                               <  1:2576706/150‑1 (MQ=255)
                        tACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGgtg                                                                                                                             >  1:3180311/1‑150 (MQ=255)
                        tACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGgtg                                                                                                                             >  1:3015397/1‑150 (MQ=255)
                          cTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACACTTCCGTGGCGCTTATTGAAGGtgtg                                                                                                                           <  2:3030232/150‑1 (MQ=255)
                          cTTGCCAGCAGATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCGTCCGTGGCGCTTATTGAAGGtgtg                                                                                                                           <  2:4226693/150‑1 (MQ=255)
                              ccAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCaa                                                                                                                       >  2:1877899/1‑150 (MQ=255)
                                    tatTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCagag                                                                                                                 >  2:4327397/1‑150 (MQ=255)
                                       tGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAgcgc                                                                                                              >  2:3958862/1‑150 (MQ=255)
                                          aTTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGt                                                                                                           <  1:1952522/150‑1 (MQ=255)
                                          aTTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGt                                                                                                           <  2:5108833/150‑1 (MQ=255)
                                          aTTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCACAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGt                                                                                                           >  2:1219509/1‑150 (MQ=255)
                                             tgttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTaaa                                                                                                        >  1:3375866/1‑150 (MQ=255)
                                              gttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAt                                                                                                       <  2:687115/150‑1 (MQ=255)
                                              gttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAt                                                                                                       >  1:496543/1‑150 (MQ=255)
                                              gttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAt                                                                                                       >  1:3287962/1‑150 (MQ=255)
                                              gttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAt                                                                                                       <  1:1568149/150‑1 (MQ=255)
                                               ttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAtc                                                                                                      >  2:5010191/1‑150 (MQ=255)
                                               ttgCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAtc                                                                                                      >  1:4915657/1‑150 (MQ=255)
                                                 gCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAAtctc                                                                                                    >  1:2839624/1‑150 (MQ=255)
                                                         cTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATg                                                                                            >  1:2652869/1‑150 (MQ=255)
                                                          tGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATgg                                                                                           <  2:2111602/150‑1 (MQ=255)
                                                           ggAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGc                                                                                          <  2:4157672/150‑1 (MQ=255)
                                                            gAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCg                                                                                         <  2:4838015/150‑1 (MQ=255)
                                                              aCAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGt                                                                                       <  1:3827068/150‑1 (MQ=255)
                                                               cAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTa                                                                                      <  2:1560936/150‑1 (MQ=255)
                                                                   cTGCACTTCCACATTGAAGCCCAAATCTTTAAACCCGTCCGCCAGCAGAGTGATTAACTCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAAttt                                                                                  <  2:4295335/150‑1 (MQ=255)
                                                                   cTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAAttt                                                                                  <  1:1610814/150‑1 (MQ=255)
                                                                     gCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGt                                                                                >  2:4387436/1‑150 (MQ=255)
                                                                      cACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGtt                                                                               >  2:4525728/1‑150 (MQ=255)
                                                                      cACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGtt                                                                               >  2:1973245/1‑150 (MQ=255)
                                                                       aCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGttt                                                                              >  1:3303981/1‑150 (MQ=255)
                                                                       aCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGttt                                                                              >  1:99469/1‑150 (MQ=255)
                                                                       aCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGGAATTTGttt                                                                              >  1:645050/1‑150 (MQ=255)
                                                                           ccACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCTGAGCGCGGTAACTCTCGATAAATGGCGGTAATTTGTTTTTCa                                                                          <  2:395004/150‑1 (MQ=255)
                                                                               aTTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCAttgt                                                                      <  1:4386507/150‑1 (MQ=255)
                                                                                 tGAAGCCCAACTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCAttgttg                                                                    <  2:3540030/150‑1 (MQ=255)
                                                                                   aaGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGac                                                                  >  1:913059/1‑150 (MQ=255)
                                                                                        cAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACAcc                                                             <  1:1135623/150‑1 (MQ=255)
                                                                                        cAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACAcc                                                             <  1:1678373/150‑1 (MQ=255)
                                                                                         aaaTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACct                                                            <  2:1737659/150‑1 (MQ=255)
                                                                                          aaTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACctc                                                           <  1:1815454/150‑1 (MQ=255)
                                                                                           aTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACctct                                                          <  2:2637398/150‑1 (MQ=255)
                                                                                           aTCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACctct                                                          <  1:1722392/150‑1 (MQ=255)
                                                                                               ttAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTc                                                      <  1:4432930/150‑1 (MQ=255)
                                                                                                tAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGACCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCa                                                     >  1:4035461/1‑150 (MQ=255)
                                                                                                 aaaCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCAt                                                    <  2:2674237/150‑1 (MQ=255)
                                                                                                 aaaCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCAt                                                    <  1:3886255/150‑1 (MQ=255)
                                                                                                   aCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                                                  <  2:3078253/150‑1 (MQ=255)
                                                                                                   aCCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                                                  >  1:3693516/1‑150 (MQ=255)
                                                                                                     cAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATa                                                >  2:816544/1‑150 (MQ=255)
                                                                                                      aGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAg                                               >  1:5062910/1‑150 (MQ=255)
                                                                                                      aGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAg                                               >  2:1068178/1‑150 (MQ=255)
                                                                                                       gTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGt                                              >  1:2062248/1‑150 (MQ=255)
                                                                                                        tCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTa                                             >  2:4911480/1‑150 (MQ=255)
                                                                                                           gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCa                                          >  2:1843688/1‑150 (MQ=255)
                                                                                                           gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCa                                          <  2:3326626/150‑1 (MQ=255)
                                                                                                           gCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCa                                          >  2:2137283/1‑150 (MQ=255)
                                                                                                            ccAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCaa                                         >  1:2157329/1‑150 (MQ=255)
                                                                                                            ccAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCaa                                         >  1:474740/1‑150 (MQ=255)
                                                                                                             cagcagAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAat                                        <  1:2765227/150‑1 (MQ=255)
                                                                                                                cagAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATAtt                                     <  1:2327740/150‑1 (MQ=255)
                                                                                                                     tGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGc                                <  2:3303961/150‑1 (MQ=255)
                                                                                                                       aTTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAg                              >  1:1456023/1‑150 (MQ=255)
                                                                                                                         tAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTa                            <  1:78549/150‑1 (MQ=255)
                                                                                                                         tAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTa                            <  1:3084184/150‑1 (MQ=255)
                                                                                                                          aaaTCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAt                           >  1:1191887/1‑150 (MQ=255)
                                                                                                                           aaTCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAtt                          >  1:108583/1‑150 (MQ=255)
                                                                                                                           aaTCTGCATTGCTTTGATCGAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTAtt                          <  2:277728/150‑1 (MQ=255)
                                                                                                                             tCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTa                        <  2:2339607/150‑1 (MQ=255)
                                                                                                                             tCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTa                        <  1:4628630/150‑1 (MQ=255)
                                                                                                                              cTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTAt                       <  2:3393698/150‑1 (MQ=255)
                                                                                                                                gCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGa                     <  1:1231864/150‑1 (MQ=255)
                                                                                                                                  aTTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACATCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAAt                   <  1:3243606/150‑1 (MQ=255)
                                                                                                                                   ttGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATa                  <  2:1493723/150‑1 (MQ=255)
                                                                                                                                   ttGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACATCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATa                  <  2:495358/150‑1 (MQ=255)
                                                                                                                                    tGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaa                 <  2:2541387/150‑1 (MQ=255)
                                                                                                                                     gCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaaa                >  2:3224921/1‑150 (MQ=255)
                                                                                                                                      cTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATaaaa               >  1:249356/1‑150 (MQ=255)
                                                                                                                                        ttGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAATATTTATGAATAAAAAt             <  2:1753429/150‑1 (MQ=255)
                                                                                                                                         tGATCGAGTGCCTCTTCCGTGGGGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTATTTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATAAATATTCATGAAGTATTTATGAATAAAAATa            <  1:4777772/150‑1 (MQ=255)
                                                                                                                                           aTCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATaca          <  1:705077/150‑1 (MQ=255)
                                                                                                                                             cGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACt        <  2:57557/150‑1 (MQ=255)
                                                                                                                                              gAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATacacaa       >  2:789292/1‑148 (MQ=255)
                                                                                                                                              gAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTa       >  2:292736/1‑150 (MQ=255)
                                                                                                                                              gAGTGACTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTa       >  1:3334524/1‑150 (MQ=255)
                                                                                                                                               aGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTaa      <  1:844508/150‑1 (MQ=255)
                                                                                                                                                 tGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACATCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACg    <  2:1330617/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTT  >  CP009273/4144550‑4144846

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: