Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,091 T→A 5.8% F179I (TTT→ATT)  pagP → phospholipid:lipid A palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,0910TA 5.8% 463.3 / 12.0 148F179I (TTT→ATT) pagPphospholipid:lipid A palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base T (77/62);  new base A (4/5);  total (81/67)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

Only 100 of 152 total aligned reads displayed.

GATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGT  >  NC_000913/656948‑657239
                                                                                                                                               |                                                                                                                                                    
 aTTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACCTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAAGTACACATTGcca                                                                                                                                                <  1:1530732/149‑2 (MQ=255)
    cACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAAGTACACATTGccattg                                                                                                                                             <  1:1220563/149‑5 (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)
                    gCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTg                                                                                                                             >  1:105889/1‑149 (MQ=255)
                      aCGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGACTGGATGCGCTTTCAGTTTTgag                                                                                                                           >  1:900959/1‑149 (MQ=255)
                       cgcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTgaga                                                                                                                          <  2:1199017/149‑1 (MQ=255)
                        gcgATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGAc                                                                                                                         <  2:499883/149‑1 (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:1596379/149‑1 (MQ=255)
                                   aaTTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGttt                                                                                                              >  1:29767/1‑149 (MQ=255)
                                     ttACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTa                                                                                                            <  1:491090/149‑1 (MQ=255)
                                      tACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAg                                                                                                           >  1:797757/1‑149 (MQ=255)
                                         aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATTGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                        <  1:1377530/149‑1 (MQ=255)
                                         aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                        <  1:335830/149‑1 (MQ=255)
                                         aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                        <  2:1448310/149‑1 (MQ=255)
                                         aTCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTaa                                                                                                        >  1:1094152/1‑149 (MQ=255)
                                          tCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAAc                                                                                                       >  1:220437/1‑149 (MQ=255)
                                           cccTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAAct                                                                                                      >  2:1095382/1‑149 (MQ=255)
                                            ccTCTCCCGGGTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActt                                                                                                     <  1:1390717/149‑1 (MQ=255)
                                             ctctCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtc                                                                                                    >  1:995915/1‑149 (MQ=255)
                                              tctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtct                                                                                                   >  1:773148/1‑149 (MQ=255)
                                               ctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtctt                                                                                                  <  2:1277462/149‑1 (MQ=255)
                                               ctcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAActtctt                                                                                                  <  2:523320/149‑1 (MQ=255)
                                                tcCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGCCCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCttt                                                                                                 >  1:783969/1‑149 (MQ=255)
                                                 cccGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTa                                                                                                <  1:1004112/149‑1 (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)
                                                                  tGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAAt                                                                               <  1:643109/149‑1 (MQ=255)
                                                                    gCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATaa                                                                             <  2:1195755/149‑1 (MQ=255)
                                                                     ccTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAg                                                                            <  2:1381756/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                                                                       >  2:1270500/1‑149 (MQ=255)
                                                                          gTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACa                                                                       <  1:1451825/149‑1 (MQ=255)
                                                                              gTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaa                                                                   >  1:785452/1‑149 (MQ=255)
                                                                              gTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCaa                                                                   >  2:61998/1‑149 (MQ=255)
                                                                                    gCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATaa                                                             <  2:1550826/149‑1 (MQ=255)
                                                                                     cccAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAAc                                                            >  2:1136909/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:251034/1‑149 (MQ=255)
                                                                                              ttttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACtttt                                                   <  2:428666/149‑1 (MQ=255)
                                                                                              ttttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGACTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACtttt                                                   >  1:36027/1‑149 (MQ=255)
                                                                                                ttCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATAAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTAt                                                 <  2:1357588/149‑1 (MQ=255)
                                                                                                  cAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCa                                               <  1:1169607/149‑1 (MQ=255)
                                                                                                    gATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACt                                             >  2:1343420/1‑149 (MQ=255)
                                                                                                       gACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACtttt                                          <  2:269133/149‑1 (MQ=255)
                                                                                                         ccTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCATTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTa                                        <  2:1063272/149‑1 (MQ=255)
                                                                                                          cTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAg                                       >  1:684161/1‑149 (MQ=255)
                                                                                                             cATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTCTTAGTaa                                    >  1:440461/1‑149 (MQ=255)
                                                                                                              aTTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAACTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTaaa                                   <  1:1499831/149‑1 (MQ=255)
                                                                                                               ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCTATATCTAAAATAAGTAACATCAAAAATAACGCGTCTTTTATCACTTTTTAGTAAAg                                  <  2:105889/149‑1 (MQ=255)
                                                                                                               ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAg                                  >  1:867222/1‑149 (MQ=255)
                                                                                                               ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAg                                  >  1:1480281/1‑149 (MQ=255)
                                                                                                               ttCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTAAGTAAAg                                  >  2:751737/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)
                                                                                                                       aCCTACAACAATGGCAATGTGTACTTTGCCTTGATGAGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTg                          <  1:1121033/149‑1 (MQ=255)
                                                                                                                       aCCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTg                          <  2:773148/149‑1 (MQ=255)
                                                                                                                        ccTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTgg                         >  2:809521/1‑149 (MQ=255)
                                                                                                                         cTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                        <  2:628376/149‑1 (MQ=255)
                                                                                                                         cTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGa                        >  1:1099670/1‑149 (MQ=255)
                                                                                                                          tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGaa                       >  2:1366890/1‑148 (MQ=255)
                                                                                                                          tACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGAc                       >  2:1507336/1‑149 (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: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)
                                                                                                                                aaTGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCg                 >  1:1042577/1‑149 (MQ=255)
                                                                                                                                  aggCAAAGTACACATTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTa               >  2:408410/2‑149 (MQ=255)
                                                                                                                                   ggCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTAc              <  1:1194954/149‑1 (MQ=255)
                                                                                                                                     cAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACca            <  1:1472720/149‑1 (MQ=255)
                                                                                                                                           gtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTg      <  1:1095382/149‑1 (MQ=255)
                                                                                                                                           gtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTg      <  2:867222/149‑1 (MQ=255)
                                                                                                                                           gtACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTNAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTg      >  1:1595412/1‑149 (MQ=255)
                                                                                                                                             aCTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgt    >  1:931065/1‑149 (MQ=255)
                                                                                                                                               tttGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgtgt  >  2:1559316/1‑149 (MQ=255)
                                                                                                                                               tttGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgtgt  >  1:1434716/1‑149 (MQ=255)
                                                                                                                                               tttGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGgtgt  >  2:849951/1‑149 (MQ=255)
                                                                                                                                               tttGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTATTATCACTTTTTAGTAAAGTTACACTGGAAAAAGCGTACCACAATTTgtgt  >  1:770419/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTTTGCCTGGATGCGCTTTCAGTTTTGAGACAAATGAAGTTTTAGTAACTTCTTTAAAATCAATAGCTAAAATAAGTAACATCAAAAATAACGCGACTTTTATCACTTTTTAGTAAAGTTACACTGGACAAAGCGTACCACAATTGGTGT  >  NC_000913/656948‑657239

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: