Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,057,226 G→T 14.3% E839* (GAA→TAA)  torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,057,2260GT14.3% 683.1 / 60.7 221E839* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (98/91);  new base T (15/16);  total (114/107)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 221 total aligned reads displayed.

AACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057082‑1057372
                                                                                                                                                |                                                                                                                                                  
aaCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAata                                                                                                                                                <  2:4075280/149‑1 (MQ=255)
 aCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAatat                                                                                                                                               <  2:819235/149‑1 (MQ=255)
  cGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATg                                                                                                                                              <  1:607179/149‑1 (MQ=255)
    tgtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGtt                                                                                                                                            <  1:910312/149‑1 (MQ=255)
     gtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                                                                                                                           <  2:3558598/149‑1 (MQ=255)
     gtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                                                                                                                           <  1:2709416/149‑1 (MQ=255)
       tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                                                                                                                         >  2:1441571/1‑149 (MQ=255)
       tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                                                                                                                         <  1:4093422/149‑1 (MQ=255)
       tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                                                                                                                         <  2:619088/149‑1 (MQ=255)
          ccATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCcgcg                                                                                                                                      >  2:1194158/1‑149 (MQ=255)
          ccATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCcgcg                                                                                                                                      >  1:352169/1‑149 (MQ=255)
           cATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGt                                                                                                                                     <  2:264245/149‑1 (MQ=255)
            aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTc                                                                                                                                    <  2:1989147/149‑1 (MQ=255)
            aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTc                                                                                                                                    <  2:1491899/149‑1 (MQ=255)
            aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTc                                                                                                                                    <  2:3539197/149‑1 (MQ=255)
             tCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCg                                                                                                                                   >  2:1958275/1‑149 (MQ=255)
                aCATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                                                                                                                                >  2:2215057/1‑149 (MQ=255)
                 cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCCCACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               <  1:1178068/149‑1 (MQ=255)
                 cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               >  2:2438812/1‑149 (MQ=255)
                      gTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaaa                                                                                                                          >  1:2246425/1‑149 (MQ=255)
                        aCATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATc                                                                                                                        <  2:823239/149‑1 (MQ=255)
                         cATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCa                                                                                                                       <  1:1202793/149‑1 (MQ=255)
                         cATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCa                                                                                                                       >  1:3373523/1‑149 (MQ=255)
                          aTCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCAt                                                                                                                      <  2:3851270/149‑1 (MQ=255)
                             gCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGAc                                                                                                                   >  2:2383435/1‑149 (MQ=255)
                              cAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGAcc                                                                                                                  <  2:1454705/149‑1 (MQ=255)
                               aGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCa                                                                                                                 >  1:241609/1‑149 (MQ=255)
                                  tGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGc                                                                                                              <  1:3305144/149‑1 (MQ=255)
                                  tGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGc                                                                                                              <  1:3658685/149‑1 (MQ=255)
                                  tGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGc                                                                                                              <  2:3620466/149‑1 (MQ=255)
                                         ggCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGc                                                                                                       >  1:1400408/1‑149 (MQ=255)
                                          gCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGca                                                                                                      >  1:1462937/1‑149 (MQ=255)
                                              ccAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCacaac                                                                                                  <  1:3458341/149‑1 (MQ=255)
                                                aGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAg                                                                                                >  2:1630117/1‑149 (MQ=255)
                                                aGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAg                                                                                                >  1:591077/1‑149 (MQ=255)
                                                 gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGa                                                                                               <  1:2446231/149‑1 (MQ=255)
                                                  tGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGAt                                                                                              <  2:2541398/149‑1 (MQ=255)
                                                   gcgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACTCTGACAGCACAACAGAtt                                                                                             >  1:2316648/1‑149 (MQ=255)
                                                      cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGcc                                                                                          <  1:3808586/149‑1 (MQ=255)
                                                      cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGcc                                                                                          <  1:2164632/149‑1 (MQ=255)
                                                      cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGcc                                                                                          >  2:2206962/1‑149 (MQ=255)
                                                         acTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgt                                                                                       >  1:267462/1‑149 (MQ=255)
                                                         acTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgt                                                                                       <  2:190403/149‑1 (MQ=255)
                                                           tACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTgtgt                                                                                     <  2:3347969/149‑1 (MQ=255)
                                                              gCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTa                                                                                  <  2:2784393/149‑1 (MQ=255)
                                                                    ggAAATTGAGCAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTg                                                                            <  1:713522/149‑1 (MQ=255)
                                                                     gAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTgg                                                                           <  1:820580/149‑1 (MQ=255)
                                                                     gAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTgg                                                                           >  2:542450/1‑149 (MQ=255)
                                                                      aaaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAAGGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGc                                                                          >  2:3257896/1‑149 (MQ=255)
                                                                      aaaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGc                                                                          <  2:1616807/149‑1 (MQ=255)
                                                                       aaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCt                                                                         >  2:91657/1‑149 (MQ=255)
                                                                             gaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCa                                                                   >  2:2674102/1‑149 (MQ=255)
                                                                             gaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCa                                                                   >  1:3185509/1‑149 (MQ=255)
                                                                               aGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgt                                                                 <  1:1441571/149‑1 (MQ=255)
                                                                                gTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgtt                                                                <  2:933232/149‑1 (MQ=255)
                                                                                gTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgtt                                                                >  2:2781357/1‑149 (MQ=255)
                                                                                    aaCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtc                                                            >  2:3781342/1‑149 (MQ=255)
                                                                                        gAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTccc                                                        >  1:2600256/1‑149 (MQ=255)
                                                                                        gAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTccc                                                        >  1:3771080/1‑149 (MQ=255)
                                                                                          aCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGt                                                      >  1:1589916/1‑149 (MQ=255)
                                                                                           cAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGCATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     >  1:2260043/1‑149 (MQ=255)
                                                                                           cAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     >  1:2923156/1‑149 (MQ=255)
                                                                                           cAGTGGAGCAGGTGACGGCGTTTAACGGACCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     >  2:2269671/1‑149 (MQ=255)
                                                                                            aGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGa                                                    >  1:1794970/1‑149 (MQ=255)
                                                                                                gAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTg                                                >  1:159411/1‑149 (MQ=255)
                                                                                                 aGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTgg                                               <  1:35421/149‑1 (MQ=255)
                                                                                                  gCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGa                                              <  2:3111697/149‑1 (MQ=255)
                                                                                                  gCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGa                                              <  2:1256748/149‑1 (MQ=255)
                                                                                                     ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGa                                           >  1:3704171/1‑149 (MQ=255)
                                                                                                       tGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATg                                         >  1:2623605/1‑149 (MQ=255)
                                                                                                         aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaa                                       <  2:2313344/149‑1 (MQ=255)
                                                                                                         aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaa                                       <  2:249036/149‑1 (MQ=255)
                                                                                                          cGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaac                                      <  2:600239/149‑1 (MQ=255)
                                                                                                          cGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaac                                      >  1:1632568/1‑149 (MQ=255)
                                                                                                           ggCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaaca                                     >  1:104228/1‑149 (MQ=255)
                                                                                                            gCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaacaa                                    >  2:2971537/1‑149 (MQ=255)
                                                                                                              gTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACt                                  >  1:2760323/1‑149 (MQ=255)
                                                                                                               tttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTg                                 >  2:1161387/1‑149 (MQ=255)
                                                                                                               tttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTg                                 >  2:3213104/1‑149 (MQ=255)
                                                                                                                   aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGc                             >  2:169785/1‑149 (MQ=255)
                                                                                                                   aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGc                             >  1:1897803/1‑149 (MQ=255)
                                                                                                                   aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGc                             >  1:960484/1‑149 (MQ=255)
                                                                                                                     ggCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaa                           >  1:1250196/1‑149 (MQ=255)
                                                                                                                      gCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaaa                          <  1:2840569/149‑1 (MQ=255)
                                                                                                                         ccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCg                       <  2:3711572/149‑1 (MQ=255)
                                                                                                                          cGTGGCGATGGTGGCGCAGGGCGACTCTGTGCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATGGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGc                      <  2:4017847/149‑1 (MQ=255)
                                                                                                                          cGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGc                      <  2:1980214/149‑1 (MQ=255)
                                                                                                                           gTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGcc                     <  2:1744843/149‑1 (MQ=255)
                                                                                                                             ggAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAg                   >  2:2730623/1‑149 (MQ=255)
                                                                                                                                 aTGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTgcg               >  2:3768108/1‑149 (MQ=255)
                                                                                                                                  tggtggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTgcgc              >  1:588152/1‑149 (MQ=255)
                                                                                                                                   ggtggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCa             <  1:521086/149‑1 (MQ=255)
                                                                                                                                     tggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGa           <  1:2214895/149‑1 (MQ=255)
                                                                                                                                     tggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGa           <  2:299468/149‑1 (MQ=255)
                                                                                                                                     tggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGa           <  2:299324/149‑1 (MQ=255)
                                                                                                                                      ggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGAt          <  1:2484473/149‑1 (MQ=255)
                                                                                                                                         gcAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGc       <  2:2481449/149‑1 (MQ=255)
                                                                                                                                            gTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGa    >  2:1412885/1‑149 (MQ=255)
                                                                                                                                             tGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGaa   >  2:1093461/1‑149 (MQ=255)
                                                                                                                                             tGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGaa   <  2:204047/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057082‑1057372

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: