Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 652,442 T→A 15.6% 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.6% 777.9 / 122.4 276Y144N (TAC→AAC) pagPphospholipid:lipid A palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base T (97/136);  new base A (19/24);  total (116/160)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 284 total aligned reads displayed.

GTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGT  >  CP009273/652294‑652591
                                                                                                                                                    |                                                                                                                                                     
gTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTa                                                                                                                                                      >  2:2454860/1‑150 (MQ=255)
  aTGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACa                                                                                                                                                    >  2:1968478/1‑150 (MQ=255)
     ccATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATcc                                                                                                                                                 <  1:3668479/150‑1 (MQ=255)
      cATGGCATTTAGGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATccc                                                                                                                                                <  1:1481467/150‑1 (MQ=255)
        tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCctc                                                                                                                                              <  1:3773134/150‑1 (MQ=255)
        tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCctc                                                                                                                                              >  2:3424374/1‑150 (MQ=255)
        tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCctc                                                                                                                                              <  1:337896/150‑1 (MQ=255)
        tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCctc                                                                                                                                              >  2:3249106/1‑150 (MQ=255)
          gCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCtctc                                                                                                                                            <  2:4249905/150‑1 (MQ=255)
            aTTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTccc                                                                                                                                          <  1:968041/150‑1 (MQ=255)
               tAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGt                                                                                                                                       <  1:4228221/150‑1 (MQ=255)
                  ggACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCt                                                                                                                                    <  1:662267/150‑1 (MQ=255)
                  ggACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCt                                                                                                                                    <  1:1936100/150‑1 (MQ=255)
                    aCTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTAc                                                                                                                                  >  1:1339903/1‑150 (MQ=255)
                     cTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACt                                                                                                                                 <  1:4100325/150‑1 (MQ=255)
                      tCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTg                                                                                                                                <  2:1295177/150‑1 (MQ=255)
                         tGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGACa                                                                                                                             <  1:3193814/150‑1 (MQ=255)
                         tGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGACa                                                                                                                             <  1:3193717/150‑1 (MQ=255)
                             aCAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTgg                                                                                                                         <  1:2672154/150‑1 (MQ=255)
                              cAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGc                                                                                                                        <  1:1642184/150‑1 (MQ=255)
                               aaaTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGcc                                                                                                                       >  1:3710108/1‑150 (MQ=255)
                               aaaTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGcc                                                                                                                       <  1:1443233/150‑1 (MQ=255)
                                 aTGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTc                                                                                                                     >  2:3409694/1‑150 (MQ=255)
                                      aaCCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTgg                                                                                                                >  2:206722/1‑150 (MQ=255)
                                         cGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGtt                                                                                                             >  2:2134517/1‑150 (MQ=255)
                                          gATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTa                                                                                                            >  2:1277286/1‑150 (MQ=255)
                                           aTTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTAt                                                                                                           <  2:2240697/150‑1 (MQ=255)
                                                cGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGccc                                                                                                      <  1:3051167/150‑1 (MQ=255)
                                                 ggATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCa                                                                                                     <  2:4388303/150‑1 (MQ=255)
                                                  gATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAg                                                                                                    >  2:702613/1‑150 (MQ=255)
                                                      cGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGAc                                                                                                <  1:2423809/150‑1 (MQ=255)
                                                       ggATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACt                                                                                               >  2:36687/1‑150 (MQ=255)
                                                       ggATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACt                                                                                               >  2:2781349/1‑150 (MQ=255)
                                                       ggATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACt                                                                                               <  1:1042715/150‑1 (MQ=255)
                                                       ggATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACt                                                                                               <  1:1458632/150‑1 (MQ=255)
                                                           gggAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTc                                                                                           <  1:4512664/150‑1 (MQ=255)
                                                              aaaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGa                                                                                        >  2:3922097/1‑150 (MQ=255)
                                                               aaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGAt                                                                                       <  2:4603941/150‑1 (MQ=255)
                                                                aGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATg                                                                                      <  1:3774378/150‑1 (MQ=255)
                                                                    ccTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCt                                                                                  >  1:39103/1‑150 (MQ=255)
                                                                    ccTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCt                                                                                  <  2:771470/150‑1 (MQ=255)
                                                                     cTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTa                                                                                 <  2:1190823/150‑1 (MQ=255)
                                                                       ggCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACa                                                                               <  2:472581/150‑1 (MQ=255)
                                                                       ggCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACa                                                                               <  2:4263325/150‑1 (MQ=255)
                                                                          gACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                                                            >  2:2738393/1‑150 (MQ=255)
                                                                          gACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                                                            >  2:3093188/1‑150 (MQ=255)
                                                                          gACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                                                            >  1:132181/1‑150 (MQ=255)
                                                                          gACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                                                            >  1:1268131/1‑150 (MQ=255)
                                                                                tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTa                                                                      >  2:2664183/1‑150 (MQ=255)
                                                                                tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTa                                                                      >  1:3300234/1‑150 (MQ=255)
                                                                                 ggCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAc                                                                     <  2:3417487/150‑1 (MQ=255)
                                                                                      aTGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTaca                                                                >  1:3888567/1‑150 (MQ=255)
                                                                                        gAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaac                                                              >  1:722175/1‑150 (MQ=255)
                                                                                        gAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaac                                                              >  1:2088484/1‑150 (MQ=255)
                                                                                         aaaaTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaaca                                                             <  1:272765/150‑1 (MQ=255)
                                                                                         aaaaTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaaca                                                             <  1:4609240/150‑1 (MQ=255)
                                                                                            aTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATg                                                          <  1:3430344/150‑1 (MQ=255)
                                                                                             ttttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATgg                                                         >  1:2547070/1‑150 (MQ=255)
                                                                                               ttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCa                                                       >  2:3705719/1‑150 (MQ=255)
                                                                                               ttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCa                                                       >  1:1605413/1‑150 (MQ=255)
                                                                                                tCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCaa                                                      <  2:407372/150‑1 (MQ=255)
                                                                                                  aTTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAtg                                                    >  1:871862/1‑150 (MQ=255)
                                                                                                      aGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTa                                                >  1:3045465/1‑150 (MQ=255)
                                                                                                      aGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGCTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTa                                                <  2:3995727/150‑1 (MQ=255)
                                                                                                       ggTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTAc                                               >  2:2490219/1‑150 (MQ=255)
                                                                                                        gTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACt                                              <  2:4639301/150‑1 (MQ=255)
                                                                                                        gTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACt                                              >  2:3868475/1‑150 (MQ=255)
                                                                                                          cTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACttt                                            <  2:3622294/150‑1 (MQ=255)
                                                                                                           tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTg                                           <  2:831574/150‑1 (MQ=255)
                                                                                                           tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTg                                           <  2:1049362/150‑1 (MQ=255)
                                                                                                            gggATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGc                                          >  1:3497634/1‑150 (MQ=255)
                                                                                                                  cACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGAt                                    <  1:1467559/150‑1 (MQ=255)
                                                                                                                   aCCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATg                                   >  1:332709/1‑150 (MQ=255)
                                                                                                                   aCCGCTGGCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATg                                   >  1:4079652/1‑150 (MQ=255)
                                                                                                                    ccGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgc                                  <  1:163147/150‑1 (MQ=255)
                                                                                                                     cGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgcg                                 >  2:1935482/1‑150 (MQ=255)
                                                                                                                         ggCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCttt                             <  2:3463736/150‑1 (MQ=255)
                                                                                                                         ggCGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCttt                             <  1:206722/150‑1 (MQ=255)
                                                                                                                           cGTAACGGCACGCGATAACTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCa                           <  2:4157375/150‑1 (MQ=255)
                                                                                                                            gTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAg                          <  1:1554412/150‑1 (MQ=255)
                                                                                                                                cGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGtttt                      >  1:2622617/1‑150 (MQ=255)
                                                                                                                                 ggCACGCGATTACTGGAAGTTCATCCCTCTCCCGGTTCTACTGCCATGGGCCTCCGTGGGTTATGGCCCAGTGTCTTTTCAGAGGAGCTACATTCCGGGTACCTCCAACAATGGCAATGTGTACTTTGCATGGATGCGCTGTCAGTTTTg                     <  2:1494851/150‑1 (MQ=255)
                                                                                                                                   cACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTgag                   <  1:2781393/150‑1 (MQ=255)
                                                                                                                                     cgcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGAc                 >  1:1727638/1‑150 (MQ=255)
                                                                                                                                      gcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACa                >  2:641773/1‑150 (MQ=255)
                                                                                                                                         aTAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAAt             <  2:1015418/150‑1 (MQ=255)
                                                                                                                                          tAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATg            <  1:427198/150‑1 (MQ=255)
                                                                                                                                            aCTGGAATAACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGaa          <  2:1639634/150‑1 (MQ=255)
                                                                                                                                             cTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg         >  1:2163390/1‑150 (MQ=255)
                                                                                                                                             cTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg         >  1:4454892/1‑150 (MQ=255)
                                                                                                                                               ggAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGtt       <  1:4334262/150‑1 (MQ=255)
                                                                                                                                                gAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGttt      >  1:2488265/1‑150 (MQ=255)
                                                                                                                                                gAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGttt      >  2:1759965/1‑150 (MQ=255)
                                                                                                                                                  aTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTa    <  1:4344178/150‑1 (MQ=255)
                                                                                                                                                   ttACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAg   >  1:1108457/1‑150 (MQ=255)
                                                                                                                                                    tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  <  1:2999249/150‑1 (MQ=255)
                                                                                                                                                    tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  <  2:3032530/150‑1 (MQ=255)
                                                                                                                                                    tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  >  2:3539502/1‑150 (MQ=255)
                                                                                                                                                    tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  <  2:234481/150‑1 (MQ=255)
                                                                                                                                                    tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAAGGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGt  <  2:2820540/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGT  >  CP009273/652294‑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: