Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,088 T→C 5.6% Y178H (TAC→CAC)  pagP → phospholipid:lipid A palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,0880TC 5.6% 444.8 / 11.1 142Y178H (TAC→CAC) pagPphospholipid:lipid A palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base T (74/60);  new base C (3/5);  total (77/65)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 146 total aligned reads displayed.

GTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTG  >  NC_000913/656942‑657235
                                                                                                                                                  |                                                                                                                                                    
gTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTAc                                                                                                                                                    <  1:448503/149‑1 (MQ=255)
      gATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGcc                                                                                                                                              >  2:884655/1‑149 (MQ=255)
       aTTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAAGTACACATTGcca                                                                                                                                             <  2:1554232/149‑2 (MQ=255)
        ttCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAAGTACACATTGccat                                                                                                                                            <  1:1325787/149‑3 (MQ=255)
        ttCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCAGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAAGTACACATTGccat                                                                                                                                            <  2:1070837/149‑3 (MQ=255)
           aCCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGACGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGAt                                                                                                                                         >  2:605434/1‑149 (MQ=255)
            ccGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATg                                                                                                                                        <  1:587799/149‑1 (MQ=255)
             cGCTGGCGTAACGGCACGCGATTACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgc                                                                                                                                       <  2:185650/149‑1 (MQ=255)
               cTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATgcgc                                                                                                                                     <  1:1279854/149‑1 (MQ=255)
                  gCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCttt                                                                                                                                  >  1:1515988/1‑149 (MQ=255)
                     tAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGCCCTACATTCCGGGTACCTACAACAATGGCAATGTGTACATTGCCTGGATGCGCTTTCAg                                                                                                                               >  2:71785/1‑149 (MQ=255)
                          gCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTg                                                                                                                          >  1:105889/1‑149 (MQ=255)
                             cgcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTgaga                                                                                                                       <  2:1199017/149‑1 (MQ=255)
                              gcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGAc                                                                                                                      >  1:1392382/1‑149 (MQ=255)
                                gATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACaa                                                                                                                    >  1:454091/1‑149 (MQ=255)
                                 aTAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCTGTGACTTTTCAGATGACCTTCATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGGCTGGATGGGCTTTCAGTTTTGAGTCaaa                                                                                                                   >  2:114865/1‑149 (MQ=255)
                                 aTAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACaaa                                                                                                                   <  1:700008/149‑1 (MQ=255)
                                      tGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAg                                                                                                              >  1:1566717/1‑149 (MQ=255)
                                         aaTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGttt                                                                                                           >  1:29767/1‑149 (MQ=255)
                                         aaTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGttt                                                                                                           <  1:1596379/149‑1 (MQ=255)
                                           ttACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTa                                                                                                         <  1:491090/149‑1 (MQ=255)
                                               aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATTGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                     <  1:1377530/149‑1 (MQ=255)
                                               aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                     <  1:335830/149‑1 (MQ=255)
                                               aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                     <  2:1448310/149‑1 (MQ=255)
                                                tCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGCCTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGTCTGCGCTATCAGTTTTGAGACAAATGAAGTTTTAGTAAc                                                                                                    >  1:78102/1‑149 (MQ=255)
                                                tCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAAc                                                                                                    >  2:1097958/1‑149 (MQ=255)
                                                tCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAAc                                                                                                    >  1:220437/1‑149 (MQ=255)
                                                 cccTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAAct                                                                                                   >  2:458190/1‑149 (MQ=255)
                                                    tctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtct                                                                                                >  1:773148/1‑149 (MQ=255)
                                                     ctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtctt                                                                                               <  2:523320/149‑1 (MQ=255)
                                                     ctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtctt                                                                                               <  2:1277462/149‑1 (MQ=255)
                                                       cccGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTa                                                                                             <  1:1004112/149‑1 (MQ=255)
                                                       cccGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTa                                                                                             >  1:1080105/1‑149 (MQ=255)
                                                        ccGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTaa                                                                                            >  2:160466/1‑149 (MQ=255)
                                                           gTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAAt                                                                                         >  2:380579/1‑149 (MQ=255)
                                                            ttCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGTTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATc                                                                                        <  1:986947/149‑1 (MQ=255)
                                                              cttcttccATTTGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCaa                                                                                      <  2:1526929/143‑1 (MQ=255)
                                                               tACTGCCATTGGCCTCCGTGGGTGATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAAt                                                                                     >  1:148127/1‑149 (MQ=255)
                                                                   gCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCt                                                                                 <  1:349118/149‑1 (MQ=255)
                                                                    ccATTGGCCTCCGTGGGTTATGGCCAAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTa                                                                               >  1:1007444/1‑149 (MQ=255)
                                                                      aTTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATTACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaa                                                                             >  2:1121033/1‑149 (MQ=255)
                                                                      aTTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaa                                                                             <  1:1445109/149‑1 (MQ=255)
                                                                       ttGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTaaaa                                                                            >  2:788184/1‑149 (MQ=255)
                                                                        tGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAAt                                                                           <  1:643109/149‑1 (MQ=255)
                                                                          gCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATaa                                                                         <  2:1195755/149‑1 (MQ=255)
                                                                           ccTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        <  2:1381756/149‑1 (MQ=255)
                                                                           ccTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        <  2:1515988/149‑1 (MQ=255)
                                                                           ccTCCGTGGGTTATGGCCCAGTGACTTTTCAGACGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        <  1:605434/149‑1 (MQ=255)
                                                                            ttCCGTTGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCTTGTACCTACAACAATGGCAATGTTTACTTTGCCTGGATTCGCTTTAAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGt                                                                       <  1:207955/148‑1 (MQ=255)
                                                                               cGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAAc                                                                    <  2:356803/149‑1 (MQ=255)
                                                                                gTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACa                                                                   <  1:1451825/149‑1 (MQ=255)
                                                                                gTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACa                                                                   >  2:1270500/1‑149 (MQ=255)
                                                                                  gggTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATc                                                                 >  1:1555653/1‑149 (MQ=255)
                                                                                    gTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaa                                                               >  1:785452/1‑149 (MQ=255)
                                                                                      tATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTAGTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaaaa                                                             >  1:788421/1‑149 (MQ=255)
                                                                                            ccAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAAcg                                                       <  1:404411/149‑1 (MQ=255)
                                                                                             cAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGAAAAATGAAGTTTTAGTAACTTCTTTAAAATAAATAGCTAAAATAAGTAACATCAAAAATAAcgc                                                      >  2:1079328/1‑149 (MQ=255)
                                                                                                   cTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACttt                                                <  1:831930/149‑1 (MQ=255)
                                                                                                    ttttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACtttt                                               <  2:428666/149‑1 (MQ=255)
                                                                                                    ttttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACtttt                                               >  2:251034/1‑149 (MQ=255)
                                                                                                    ttttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGACTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACtttt                                               >  1:36027/1‑149 (MQ=255)
                                                                                                      ttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATAAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTAt                                             <  2:1357588/149‑1 (MQ=255)
                                                                                                          gATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACt                                         >  2:1343420/1‑149 (MQ=255)
                                                                                                             gACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACtttt                                      <  2:269133/149‑1 (MQ=255)
                                                                                                                cTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAg                                   >  1:684161/1‑149 (MQ=255)
                                                                                                                   cATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTCTTAGTaa                                >  1:440461/1‑149 (MQ=255)
                                                                                                                    aTTCCTGGTACCTACAACAATGGCAATGTGTACTTTGCATGTATTCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTaaa                               <  1:1097958/149‑1 (MQ=255)
                                                                                                                    aTTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAACTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTaaa                               <  1:1499831/149‑1 (MQ=255)
                                                                                                                     ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCTATATCTAAAATAAGTAACATCAAAAATAACGCGTCTTTTATCACTTTTTAGTAAAg                              <  2:105889/149‑1 (MQ=255)
                                                                                                                     ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAg                              >  1:1480281/1‑149 (MQ=255)
                                                                                                                     ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAg                              >  1:867222/1‑149 (MQ=255)
                                                                                                                        cGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGAGACTTTTATCACTTTTTAGTAAAGTTa                           >  2:262767/1‑149 (MQ=255)
                                                                                                                         gggTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTac                          >  2:1564958/1‑149 (MQ=255)
                                                                                                                           gTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        >  2:818281/1‑101 (MQ=255)
                                                                                                                           gTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        >  2:957138/1‑101 (MQ=255)
                                                                                                                           gTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        <  1:957138/101‑1 (MQ=255)
                                                                                                                           gTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                        <  1:818281/101‑1 (MQ=255)
                                                                                                                             aCCTACAACAATGGCAATGTGTACTTTGCCTTGATGAGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTg                      <  1:1121033/149‑1 (MQ=255)
                                                                                                                             aCCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTg                      <  2:773148/149‑1 (MQ=255)
                                                                                                                               cTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                    <  2:628376/149‑1 (MQ=255)
                                                                                                                                tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                    >  1:1168878/1‑149 (MQ=255)
                                                                                                                                tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGaa                   >  2:1366890/1‑148 (MQ=255)
                                                                                                                                tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGAc                   >  2:1507336/1‑149 (MQ=255)
                                                                                                                                tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAAAATAAAAAATAACGAGACTTTTATAACTTTTTAGTAAATTTACACTTGaa                   >  1:103842/1‑148 (MQ=255)
                                                                                                                                 acaacaATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACa                  >  2:117626/1‑149 (MQ=255)
                                                                                                                                  caacaaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaa                 <  2:220437/149‑1 (MQ=255)
                                                                                                                                   aacaaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaaa                <  2:1371986/149‑1 (MQ=255)
                                                                                                                                   aacaaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaaa                >  1:1212411/1‑149 (MQ=255)
                                                                                                                                   aacaaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaaa                <  1:1564958/149‑1 (MQ=255)
                                                                                                                                   aacaaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACaaa                >  1:301982/1‑149 (MQ=255)
                                                                                                                                    atccaggCAAAGTACACATTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAg               >  1:1447773/6‑149 (MQ=255)
                                                                                                                                     tccaggCAAAGTACACATTGCCTGGAAGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGc              >  1:1395497/5‑149 (MQ=255)
                                                                                                                                       aTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGt            <  2:1007444/149‑1 (MQ=255)
                                                                                                                                        aggCAAAGTACACATTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTa           >  2:408410/2‑149 (MQ=255)
                                                                                                                                         ggCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTAc          <  1:1194954/149‑1 (MQ=255)
                                                                                                                                          gCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTAcc         >  2:827863/1‑149 (MQ=255)
                                                                                                                                           cAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACca        <  1:1472720/149‑1 (MQ=255)
                                                                                                                                               gtgtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAAt    >  2:1136924/1‑149 (MQ=255)
                                                                                                                                                 gtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTg  <  1:1095382/149‑1 (MQ=255)
                                                                                                                                                 gtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTNAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTg  >  1:1595412/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTG  >  NC_000913/656942‑657235

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