Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 652,442 T→A 15.0% Y144N (TAC→AAC)  pagP → phospholipid:lipid A palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273652,4420TA15.0% 809.1 / 121.0 286Y144N (TAC→AAC) pagPphospholipid:lipid A palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base T (95/148);  new base A (27/16);  total (122/164)
Fisher's exact test for biased strand distribution p-value = 4.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 289 total aligned reads displayed.

TGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGT  >  CP009273/652293‑652591
                                                                                                                                                     |                                                                                                                                                     
  tATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTAc                                                                                                                                                     >  1:1935039/1‑150 (MQ=255)
   aTGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACa                                                                                                                                                    >  1:1202189/1‑150 (MQ=255)
    tGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACAt                                                                                                                                                   >  1:1005947/1‑150 (MQ=255)
       cATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATccc                                                                                                                                                >  1:522074/1‑150 (MQ=255)
        aTGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCct                                                                                                                                               <  1:1279662/150‑1 (MQ=255)
           gCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCtctc                                                                                                                                            <  1:2762141/150‑1 (MQ=255)
           gCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCtctc                                                                                                                                            <  1:3398495/150‑1 (MQ=255)
            cATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTcc                                                                                                                                           >  1:993709/1‑150 (MQ=255)
            cATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTcc                                                                                                                                           <  2:2608596/150‑1 (MQ=255)
             aTTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTccc                                                                                                                                          >  2:626337/1‑150 (MQ=255)
             aTTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTccc                                                                                                                                          <  1:1115141/150‑1 (MQ=255)
              tttAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCg                                                                                                                                         >  1:4511963/1‑150 (MQ=255)
                   ggACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCt                                                                                                                                    <  1:2788030/150‑1 (MQ=255)
                      cTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACt                                                                                                                                 <  1:3079595/150‑1 (MQ=255)
                          tGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCa                                                                                                                             >  1:885798/1‑150 (MQ=255)
                            gAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCAtt                                                                                                                           <  2:3923701/150‑1 (MQ=255)
                                aaaTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGcc                                                                                                                       >  2:2850786/1‑150 (MQ=255)
                                 aaTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCt                                                                                                                      >  1:3218156/1‑150 (MQ=255)
                                  aTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTc                                                                                                                     <  2:3643385/150‑1 (MQ=255)
                                  aTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTc                                                                                                                     >  2:408967/1‑150 (MQ=255)
                                  aTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTc                                                                                                                     >  2:4554518/1‑150 (MQ=255)
                                   tGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTcc                                                                                                                    <  1:659432/150‑1 (MQ=255)
                                    gggAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCg                                                                                                                   >  2:1508413/1‑150 (MQ=255)
                                       aaCCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTgg                                                                                                                >  1:4575594/1‑150 (MQ=255)
                                          cGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGtt                                                                                                             >  1:2228352/1‑150 (MQ=255)
                                           gATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTa                                                                                                            >  1:933442/1‑150 (MQ=255)
                                           gATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTa                                                                                                            >  1:3599311/1‑150 (MQ=255)
                                           gATTGCCGGATAAGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGAACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTa                                                                                                            >  2:769292/1‑150 (MQ=255)
                                            aTTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTAt                                                                                                           <  1:1582185/150‑1 (MQ=255)
                                             ttGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATg                                                                                                          >  2:4247742/1‑150 (MQ=255)
                                                 cGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGccc                                                                                                      >  1:1513348/1‑150 (MQ=255)
                                                  ggATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCa                                                                                                     <  2:2666485/150‑1 (MQ=255)
                                                       cGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGAc                                                                                                <  1:3590207/150‑1 (MQ=255)
                                                       cGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGAc                                                                                                <  1:1996830/150‑1 (MQ=255)
                                                        ggATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACt                                                                                               <  1:1156702/150‑1 (MQ=255)
                                                         gATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACtt                                                                                              >  2:1926061/1‑150 (MQ=255)
                                                             ggAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCa                                                                                          <  2:1708550/150‑1 (MQ=255)
                                                                 aGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATg                                                                                      >  1:3953016/1‑150 (MQ=255)
                                                                 aGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATg                                                                                      <  1:3988024/150‑1 (MQ=255)
                                                                   tACCTGGCGACCGCTGGCTGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTTGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGAc                                                                                    <  2:1314809/150‑1 (MQ=255)
                                                                   tACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGAc                                                                                    <  2:1314494/150‑1 (MQ=255)
                                                                   tACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGAc                                                                                    >  2:840864/1‑150 (MQ=255)
                                                                     ccTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCt                                                                                  <  2:3236955/150‑1 (MQ=255)
                                                                      cTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTa                                                                                 >  1:2626612/1‑150 (MQ=255)
                                                                      cTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTa                                                                                 >  1:250636/1‑150 (MQ=255)
                                                                       tGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTAc                                                                                >  2:2279047/1‑150 (MQ=255)
                                                                       tGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTAc                                                                                >  2:2804985/1‑150 (MQ=255)
                                                                        ggCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACa                                                                               >  1:2841745/1‑150 (MQ=255)
                                                                         gCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACAt                                                                              <  2:1156139/150‑1 (MQ=255)
                                                                            aCCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTcc                                                                           <  2:3869393/150‑1 (MQ=255)
                                                                                 tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTa                                                                      >  2:160771/1‑150 (MQ=255)
                                                                                 tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTa                                                                      >  2:2209700/1‑150 (MQ=255)
                                                                                    cGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCt                                                                   >  2:2382594/1‑150 (MQ=255)
                                                                                     ggATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTa                                                                  <  2:3034765/150‑1 (MQ=255)
                                                                                        tGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaa                                                               <  1:160533/150‑1 (MQ=255)
                                                                                        tGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaa                                                               <  1:205002/150‑1 (MQ=255)
                                                                                             aTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATg                                                          >  2:3655145/1‑150 (MQ=255)
                                                                                             aTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATg                                                          <  2:1634451/150‑1 (MQ=255)
                                                                                              ttttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATgg                                                         >  2:3094327/1‑150 (MQ=255)
                                                                                                ttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCa                                                       >  1:1012435/1‑150 (MQ=255)
                                                                                                 tCATTTAGGTCCGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCaa                                                      >  2:4608191/1‑150 (MQ=255)
                                                                                                  cATTAAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTTGCCTCCGTGGGTTATGGCCCAGTGACTTTTCATATGACATACATTCCGGGTACCTACAACAATGGCAAt                                                     >  1:2685741/1‑150 (MQ=255)
                                                                                                    tttAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAtgt                                                   <  2:4219190/150‑1 (MQ=255)
                                                                                                       aGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTa                                                >  1:1453462/1‑150 (MQ=255)
                                                                                                        ggTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTAc                                               <  2:3710056/150‑1 (MQ=255)
                                                                                                          tCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACtt                                             <  1:3720894/150‑1 (MQ=255)
                                                                                                            tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTg                                           <  2:561990/150‑1 (MQ=255)
                                                                                                            tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTg                                           <  1:744149/150‑1 (MQ=255)
                                                                                                            tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTg                                           <  2:2570136/150‑1 (MQ=255)
                                                                                                              ggATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGcc                                         <  2:1206652/150‑1 (MQ=255)
                                                                                                                 ttCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTgg                                      <  1:3747196/150‑1 (MQ=255)
                                                                                                                    aCCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATg                                   <  2:1876936/150‑1 (MQ=255)
                                                                                                                     ccGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgc                                  <  2:3127379/150‑1 (MQ=255)
                                                                                                                     ccGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgc                                  <  2:4458180/150‑1 (MQ=255)
                                                                                                                          ggCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCttt                             <  1:891945/150‑1 (MQ=255)
                                                                                                                          ggCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCttt                             <  1:4337588/150‑1 (MQ=255)
                                                                                                                            cGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCa                           <  2:361324/150‑1 (MQ=255)
                                                                                                                              tAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCGTTCAGt                         <  2:776419/150‑1 (MQ=255)
                                                                                                                               aaCGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGtt                        <  1:1432191/150‑1 (MQ=255)
                                                                                                                                aCGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACACTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGttt                       <  1:305306/150‑1 (MQ=255)
                                                                                                                                aCGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGttt                       >  1:4142841/1‑150 (MQ=255)
                                                                                                                                aCGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGttt                       <  2:1998880/150‑1 (MQ=255)
                                                                                                                                 cGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGtttt                      <  2:1764190/150‑1 (MQ=255)
                                                                                                                                  ggCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTg                     <  2:1473606/150‑1 (MQ=255)
                                                                                                                                  ggCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTg                     <  2:2522965/150‑1 (MQ=255)
                                                                                                                                    cACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTgag                   <  2:2001675/150‑1 (MQ=255)
                                                                                                                                     aCGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTgaga                  <  2:4325364/150‑1 (MQ=255)
                                                                                                                                       gcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACa                <  1:164504/150‑1 (MQ=255)
                                                                                                                                        cgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACaa               <  1:678307/150‑1 (MQ=255)
                                                                                                                                        cgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACaa               >  2:1247819/1‑150 (MQ=255)
                                                                                                                                         gATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACaaa              >  2:336657/1‑150 (MQ=255)
                                                                                                                                           tAACTGGAAGTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATg            <  2:3888067/150‑1 (MQ=255)
                                                                                                                                              cTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg         >  1:2922164/1‑150 (MQ=255)
                                                                                                                                              cTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg         >  1:2449863/1‑150 (MQ=255)
                                                                                                                                              cTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg         <  1:3541495/150‑1 (MQ=255)
                                                                                                                                               tGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGt        <  2:1027412/150‑1 (MQ=255)
                                                                                                                                               tGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGt        <  1:4418031/150‑1 (MQ=255)
                                                                                                                                                  aaTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGtttt     >  2:3908750/1‑150 (MQ=255)
                                                                                                                                                   aTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTa    <  1:342320/150‑1 (MQ=255)
                                                                                                                                                     tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  >  2:3252672/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGT  >  CP009273/652293‑652591

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: