Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 351,447 T→A 19.5% L272Q (CTG→CAG)  codB → cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273351,4470TA19.5% 774.7 / 189.1 318L272Q (CTG→CAG) codBcytosine transporter
Reads supporting (aligned to +/- strand):  ref base T (113/138);  new base A (26/36);  total (144/174)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

ACTGGCGGTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTG  >  CP009273/351298‑351596
                                                                                                                                                     |                                                                                                                                                     
 cTGGCGGTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTg                                                                                                                                                      <  2:3520601/150‑1 (MQ=255)
  tGGCGGTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTgg                                                                                                                                                     >  2:3046264/1‑150 (MQ=255)
  tGGCGGTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGAGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTgg                                                                                                                                                     >  2:1025552/1‑150 (MQ=255)
      ggTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCt                                                                                                                                                 >  1:3262460/1‑150 (MQ=255)
      ggTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCt                                                                                                                                                 >  2:253570/1‑150 (MQ=255)
        tGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGa                                                                                                                                               <  2:2346707/150‑1 (MQ=255)
         gCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGaa                                                                                                                                              <  1:5049094/150‑1 (MQ=255)
          cTGGTGGCGATGGTGGCCTTTGTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAat                                                                                                                                             <  2:2998252/150‑1 (MQ=255)
           tggtggCTATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAata                                                                                                                                            <  1:2192222/150‑1 (MQ=255)
           tggtggCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAata                                                                                                                                            <  1:4064637/150‑1 (MQ=255)
           tggtggCGATGGTGGCCTATTTCCTCGGCAACTCGTTGATGTTTATTTTCGGGGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAata                                                                                                                                            <  2:4668634/150‑1 (MQ=255)
             gtggCGATGGGGGCCTTTTTCCTCGGCAACTCGTGGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATc                                                                                                                                          <  2:2671053/150‑1 (MQ=255)
                 cGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGa                                                                                                                                      >  1:3606034/1‑150 (MQ=255)
                   aTGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGacc                                                                                                                                    >  1:4118045/1‑150 (MQ=255)
                       tggCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGTATATCTGGAccacca                                                                                                                                <  2:2117548/150‑1 (MQ=255)
                         ggcTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAAc                                                                                                                              <  2:2993976/148‑1 (MQ=255)
                          ccTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACg                                                                                                                             <  1:628101/150‑1 (MQ=255)
                           cTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGa                                                                                                                            <  2:1700286/150‑1 (MQ=255)
                           cTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGa                                                                                                                            >  1:2173087/1‑150 (MQ=255)
                             ttttCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATa                                                                                                                          >  2:1799312/1‑150 (MQ=255)
                             ttttCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATa                                                                                                                          >  1:1066901/1‑150 (MQ=255)
                               ttCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAAc                                                                                                                        <  1:2302188/150‑1 (MQ=255)
                                tCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACg                                                                                                                       >  1:299073/1‑150 (MQ=255)
                                   tCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCAc                                                                                                                    >  1:267479/1‑150 (MQ=255)
                                     ggCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCActc                                                                                                                  >  1:3809821/1‑150 (MQ=255)
                                       cAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTa                                                                                                                >  1:3324518/1‑150 (MQ=255)
                                       cAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTa                                                                                                                >  2:734133/1‑150 (MQ=255)
                                         aCTCGTTGATGTTTATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATg                                                                                                              <  2:4663350/150‑1 (MQ=255)
                                          cTCGTTGATGTTTATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGc                                                                                                             >  2:1468612/1‑150 (MQ=255)
                                          cTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGc                                                                                                             >  2:4839308/1‑150 (MQ=255)
                                            cGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGt                                                                                                           >  2:4036506/1‑150 (MQ=255)
                                             gTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTc                                                                                                          <  1:5155422/150‑1 (MQ=255)
                                             gTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTc                                                                                                          >  2:4464186/1‑150 (MQ=255)
                                               tGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCgg                                                                                                        >  1:1944559/1‑150 (MQ=255)
                                                 aTGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGt                                                                                                      <  1:1893788/150‑1 (MQ=255)
                                                    tttATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTa                                                                                                   <  1:2567090/150‑1 (MQ=255)
                                                     ttATTTTCTGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAg                                                                                                  <  1:4485769/150‑1 (MQ=255)
                                                          ttCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTc                                                                                             <  1:2461477/150‑1 (MQ=255)
                                                             ggTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGcc                                                                                          >  1:3162443/1‑150 (MQ=255)
                                                               tGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCaa                                                                                        >  1:4247394/1‑150 (MQ=255)
                                                                gCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCGGATGTGTTGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAAc                                                                                       <  1:4676412/150‑1 (MQ=255)
                                                                     gggCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTAc                                                                                  >  1:810542/1‑150 (MQ=255)
                                                                      ggCGCTGCGGCACTGGGCCTGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTAcc                                                                                 <  1:3566487/150‑1 (MQ=255)
                                                                          cTGCGGCACTGGGCATGACGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGa                                                                             <  2:1100856/150‑1 (MQ=255)
                                                                          cTGCGGCACTGGGCATGACGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGa                                                                             <  2:1449898/150‑1 (MQ=255)
                                                                           tGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGAt                                                                            <  1:452333/150‑1 (MQ=255)
                                                                              ggCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTc                                                                         >  1:3750351/1‑150 (MQ=255)
                                                                              ggCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTc                                                                         <  2:4099986/150‑1 (MQ=255)
                                                                                 aCTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAg                                                                      <  1:2714926/150‑1 (MQ=255)
                                                                                 aCTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAg                                                                      <  1:1304252/150‑1 (MQ=255)
                                                                                  cTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGc                                                                     >  2:3671720/1‑150 (MQ=255)
                                                                                  cTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGc                                                                     >  2:2440336/1‑150 (MQ=255)
                                                                                   tGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCa                                                                    >  1:2033361/1‑150 (MQ=255)
                                                                                     ggCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCaaa                                                                  >  1:5119563/1‑150 (MQ=255)
                                                                                     ggCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCCGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCaaa                                                                  >  1:4332671/1‑150 (MQ=255)
                                                                                           gCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCtt                                                            <  2:948625/150‑1 (MQ=255)
                                                                                             ggATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTc                                                          <  2:1780769/150‑1 (MQ=255)
                                                                                               atatCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCgg                                                        >  1:3718783/1‑150 (MQ=255)
                                                                                                 atCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGtc                                                      >  2:5116196/1‑149 (MQ=255)
                                                                                                 atCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTa                                                      <  2:4238227/150‑1 (MQ=255)
                                                                                                    tctGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATc                                                   >  2:4392129/1‑150 (MQ=255)
                                                                                                     ctGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCa                                                  >  2:5154827/1‑150 (MQ=255)
                                                                                                     ctGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCa                                                  >  1:86993/1‑150 (MQ=255)
                                                                                                       gatgtgatgATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAAc                                                <  1:1605518/150‑1 (MQ=255)
                                                                                                          gtgatgATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCCGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGt                                             >  1:4208374/1‑150 (MQ=255)
                                                                                                          gtgatgATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCCGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGt                                             >  1:5098865/1‑150 (MQ=255)
                                                                                                            gttgtttGCTCAGTGCCTGTTGATGCCTGCGATTGTGGTGCTGGGGCTGATTATCTGGACCACCAACGATAACGCACTCTATGCGTCGTTTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtat                                           <  2:2771948/145‑1 (MQ=255)
                                                                                                            gatgatTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtat                                           >  2:4937098/1‑150 (MQ=255)
                                                                                                            gatgatTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtat                                           <  2:4481875/150‑1 (MQ=255)
                                                                                                             atgatTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtatt                                          >  1:1323874/1‑150 (MQ=255)
                                                                                                               gatTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtattat                                        >  2:1245369/1‑150 (MQ=255)
                                                                                                               gatTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGtattat                                        >  2:4016996/1‑150 (MQ=255)
                                                                                                                atTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATc                                       >  2:3389351/1‑150 (MQ=255)
                                                                                                                 ttGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCg                                      <  2:2203436/150‑1 (MQ=255)
                                                                                                                 ttGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCg                                      <  1:723433/150‑1 (MQ=255)
                                                                                                                 ttGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCg                                      <  1:5080115/150‑1 (MQ=255)
                                                                                                                  tGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCgg                                     <  2:296098/150‑1 (MQ=255)
                                                                                                                     tCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTAc                                  <  1:2385334/150‑1 (MQ=255)
                                                                                                                      ctgggCCTGCTGCTGCCTGCGTTTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACg                                 <  2:4219108/148‑1 (MQ=255)
                                                                                                                      cAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGTTACg                                 >  2:3116247/1‑150 (MQ=255)
                                                                                                                      cAGGGCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACg                                 <  2:684042/150‑1 (MQ=255)
                                                                                                                        gggCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGt                               <  1:3996897/150‑1 (MQ=255)
                                                                                                                         ggCCTGCTGCTGCCTGCGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTc                              <  1:4434369/150‑1 (MQ=255)
                                                                                                                          gCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCt                             >  1:3418533/1‑150 (MQ=255)
                                                                                                                          gCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCt                             <  1:950606/150‑1 (MQ=255)
                                                                                                                           ccTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTg                            <  1:5205952/150‑1 (MQ=255)
                                                                                                                             tgctgctgcCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTgcg                          >  2:713086/1‑150 (MQ=255)
                                                                                                                               ctgctgcCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCa                        <  1:2906542/150‑1 (MQ=255)
                                                                                                                               ctgctgcCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCa                        >  1:1240883/1‑150 (MQ=255)
                                                                                                                                 gctgcCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCAtt                      <  1:4470102/150‑1 (MQ=255)
                                                                                                                                 gctgcCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCAtt                      <  1:3707103/150‑1 (MQ=255)
                                                                                                                                   tgcctgcGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTAt                    >  2:3837784/1‑150 (MQ=255)
                                                                                                                                   tgcctgcGATTGTGGTGCTGGGGATGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTAt                    >  2:1970533/1‑150 (MQ=255)
                                                                                                                                      ctgcGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGc                 >  1:489382/1‑150 (MQ=255)
                                                                                                                                      ctgcGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGc                 <  2:2971021/150‑1 (MQ=255)
                                                                                                                                       tgcGATTGTGGTGCAGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACAGTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCt                <  2:3606034/150‑1 (MQ=255)
                                                                                                                                          caTTGTGTTGATGGGGCTGAAAATCTGGACCACCAAAGATAACTCACTCTATGCCTCGGGTTTAGGATTCGCCTACATTACCGGGATGTCGAGCAAACCCCTTTCGGTAATCACCGGTATTATCGTTACGGTCTGCGCATTATGGCTTta             <  2:3235785/149‑1 (MQ=255)
                                                                                                                                           aTTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGtat            >  1:4234082/1‑150 (MQ=255)
                                                                                                                                                ggtgCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCAGTATGGCTGTATAACaa       >  2:1277374/1‑150 (MQ=255)
                                                                                                                                                     tGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTg  <  2:2606388/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGGCGGTGCTGGTGGCGATGGTGGCCTTTTTCCTCGGCAACTCGTTGATGTTTATTTTCGGTGCAGCGGGCGCTGCGGCACTGGGCATGGCGGATATCTCTGATGTGATGATTGCTCAGGGCCTGCTGCTGCCTGCGATTGTGGTGCTGGGGCTGAATATCTGGACCACCAACGATAACGCACTCTATGCGTCGGGTTTAGGTTTCGCCAACATTACCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTG  >  CP009273/351298‑351596

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: