Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,296,865 T→A 12.7% T38S (ACC→TCC)  napF ← ferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,296,8650TA12.7% 575.1 / 61.6 196T38S (ACC→TCC) napFferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA
Reads supporting (aligned to +/- strand):  ref base T (89/82);  new base A (18/7);  total (107/89)
Fisher's exact test for biased strand distribution p-value = 8.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

GCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGA  >  CP009273/2296718‑2297012
                                                                                                                                                   |                                                                                                                                                   
gcgGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGt                                                                                                                                                     <  1:490841/148‑1 (MQ=255)
 cgGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTg                                                                                                                                                    >  1:2379623/1‑148 (MQ=255)
     gaAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAAt                                                                                                                                                <  2:415393/148‑1 (MQ=255)
      aaaaTAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATg                                                                                                                                               <  1:133927/148‑1 (MQ=255)
      aaaaTAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATg                                                                                                                                               >  2:2838839/1‑148 (MQ=255)
       aaaTAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATgg                                                                                                                                              <  1:1829030/148‑1 (MQ=255)
         aTAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGt                                                                                                                                            <  1:2848136/148‑1 (MQ=255)
         aTAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGt                                                                                                                                            >  1:350561/1‑148 (MQ=255)
                ttCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGaaaa                                                                                                                                     >  2:1254621/1‑148 (MQ=255)
                ttCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGaaaa                                                                                                                                     <  1:1472515/148‑1 (MQ=255)
                   ggggCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATga                                                                                                                                  <  1:1290150/148‑1 (MQ=255)
                     ggCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATgaga                                                                                                                                <  2:1844601/148‑1 (MQ=255)
                       cAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGAtt                                                                                                                              >  1:1112366/1‑148 (MQ=255)
                       cAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGAtt                                                                                                                              >  1:2808224/1‑148 (MQ=255)
                         aGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATtca                                                                                                                            >  2:2083548/1‑148 (MQ=255)
                                 cgcAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGa                                                                                                                    <  2:811278/148‑1 (MQ=255)
                                   cAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                                                                                  <  1:385650/148‑1 (MQ=255)
                                     ggCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAg                                                                                                                <  1:1439066/148‑1 (MQ=255)
                                     ggCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGAGTCATCACCGGACCAg                                                                                                                <  1:2760695/148‑1 (MQ=255)
                                      gCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAgg                                                                                                               >  2:1521781/1‑148 (MQ=255)
                                        gTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGc                                                                                                             <  1:62512/148‑1 (MQ=255)
                                         tAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcg                                                                                                            >  2:2744140/1‑148 (MQ=255)
                                          aaCAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgg                                                                                                           <  2:2892079/148‑1 (MQ=255)
                                           aCAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgga                                                                                                          >  2:2215622/1‑148 (MQ=255)
                                            cAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcggac                                                                                                         >  2:2724696/1‑148 (MQ=255)
                                              gAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcggacgg                                                                                                       >  1:1838775/1‑148 (MQ=255)
                                               aaGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcggacgga                                                                                                      <  1:2249034/148‑1 (MQ=255)
                                                 gCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATa                                                                                                    >  2:209833/1‑148 (MQ=255)
                                                 gCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATa                                                                                                    >  2:338799/1‑148 (MQ=255)
                                                 gCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATa                                                                                                    >  1:855142/1‑148 (MQ=255)
                                                  cTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATAc                                                                                                   >  2:2304673/1‑148 (MQ=255)
                                                    gCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCg                                                                                                 >  2:138781/1‑148 (MQ=255)
                                                           ttATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTgg                                                                                          >  1:1700844/1‑148 (MQ=255)
                                                           ttATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTgg                                                                                          >  1:2986942/1‑148 (MQ=255)
                                                            tATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGc                                                                                         <  1:1962373/148‑1 (MQ=255)
                                                             aTTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCt                                                                                        <  2:2975912/148‑1 (MQ=255)
                                                                tttGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTg                                                                                     <  1:2207566/148‑1 (MQ=255)
                                                                   gAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTgcgc                                                                                  <  1:1383418/148‑1 (MQ=255)
                                                                    aaaTTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGAGGGTGCAATGGGTCAGAAAATGAGATTAATCACCGGACCAGGGCGGACGGATACCGTTAATGGCTTTGCGcc                                                                                 >  1:1176641/1‑148 (MQ=255)
                                                                      aTTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAg                                                                               >  2:108795/1‑148 (MQ=255)
                                                                      aTTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAg                                                                               >  1:1253643/1‑148 (MQ=255)
                                                                       ttAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGc                                                                              >  1:2783310/1‑148 (MQ=255)
                                                                        tAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCg                                                                             >  2:226697/1‑148 (MQ=255)
                                                                        tAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCg                                                                             <  2:1319629/148‑1 (MQ=255)
                                                                          aCGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGAc                                                                           >  1:1570555/1‑148 (MQ=255)
                                                                                  aTAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGt                                                                   <  1:2481334/148‑1 (MQ=255)
                                                                                    aGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTAt                                                                 <  2:1636509/148‑1 (MQ=255)
                                                                                      ccgcccgcGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGc                                                               >  1:2446290/1‑148 (MQ=255)
                                                                                      ccgcccgcGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGc                                                               <  1:2687166/148‑1 (MQ=255)
                                                                                       cgcccgcGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGcc                                                              <  2:559088/148‑1 (MQ=255)
                                                                                            gcgccgcgTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGAc                                                         >  2:2172011/1‑148 (MQ=255)
                                                                                              gccgcgTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACgg                                                       <  1:3298778/148‑1 (MQ=255)
                                                                                               ccgcgTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACggg                                                      <  2:3012606/148‑1 (MQ=255)
                                                                                                  cgTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATg                                                   >  1:1598744/1‑148 (MQ=255)
                                                                                                   gTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGc                                                  <  1:2497774/148‑1 (MQ=255)
                                                                                                      gCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATc                                               <  1:3221402/148‑1 (MQ=255)
                                                                                                      gCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATc                                               >  1:2366844/1‑148 (MQ=255)
                                                                                                       cAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCa                                              >  1:2441962/1‑148 (MQ=255)
                                                                                                       cAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCa                                              <  2:2256327/148‑1 (MQ=255)
                                                                                                       cAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCa                                              >  1:1601772/1‑148 (MQ=255)
                                                                                                         gAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAAt                                            <  1:3105045/148‑1 (MQ=255)
                                                                                                          aaTGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                           <  1:2054141/148‑1 (MQ=255)
                                                                                                          aaTGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                           <  2:1720927/148‑1 (MQ=255)
                                                                                                          aaTGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                           <  2:882882/148‑1 (MQ=255)
                                                                                                          aaTGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                           <  1:2521342/148‑1 (MQ=255)
                                                                                                          aaTGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                           >  2:1140965/1‑148 (MQ=255)
                                                                                                               tgttTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTcac                                      >  2:2786608/1‑148 (MQ=255)
                                                                                                                gttTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTcaca                                     <  1:1402908/148‑1 (MQ=255)
                                                                                                                gttTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGACAGAAAATGAGATTCATAACAGCACCAGGGAGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATCCCCCGACGGGATGCATCAATCATcaca                                     >  2:46098/1‑148 (MQ=255)
                                                                                                                       caAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCt                              <  1:2239388/148‑1 (MQ=255)
                                                                                                                       caAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCt                              >  1:2272752/1‑148 (MQ=255)
                                                                                                                         aGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTc                            <  2:350561/148‑1 (MQ=255)
                                                                                                                         aGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTc                            <  2:2379623/148‑1 (MQ=255)
                                                                                                                         aGCATTAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATAACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTc                            >  2:2026032/1‑148 (MQ=255)
                                                                                                                          gCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTcc                           >  2:1859578/1‑148 (MQ=255)
                                                                                                                          gCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTcc                           >  2:3101100/1‑148 (MQ=255)
                                                                                                                          gCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTcc                           >  1:666773/1‑148 (MQ=255)
                                                                                                                           cATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCca                          >  1:215831/1‑148 (MQ=255)
                                                                                                                           cATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCca                          >  1:1496876/1‑148 (MQ=255)
                                                                                                                             tttaTACAGGCGTCGCAGCGGGAGCAATGGCTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGATTTGCGCCAGAGAACAGTGAGTATGCCCCGACTGGATGCATCAATCTTCACATTGACCTTCcatc                        <  2:290245/145‑1 (MQ=255)
                                                                                                                             ttAATACAGGCGTCGCAGCGGGTGAATTGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatc                        <  2:2287629/148‑1 (MQ=255)
                                                                                                                              tAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatca                       >  2:2954930/1‑148 (MQ=255)
                                                                                                                              tAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatca                       >  2:1653595/1‑148 (MQ=255)
                                                                                                                              tAATACAGGCGTCGCAGCGGGAGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatca                       >  1:798575/1‑148 (MQ=255)
                                                                                                                                aTACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCAtt                     >  2:2547097/1‑148 (MQ=255)
                                                                                                                                 tACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTa                    <  2:791125/148‑1 (MQ=255)
                                                                                                                                 tACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTa                    >  1:2244485/1‑148 (MQ=255)
                                                                                                                                 tACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTa                    <  2:3268485/148‑1 (MQ=255)
                                                                                                                                  aCAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTaa                   >  1:2649729/1‑148 (MQ=255)
                                                                                                                                  aCAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTaa                   <  1:2083548/148‑1 (MQ=255)
                                                                                                                                    aGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAAcg                 >  1:2393561/1‑148 (MQ=255)
                                                                                                                                    aGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAAcg                 <  1:1780001/148‑1 (MQ=255)
                                                                                                                                     ggCGTCGCATCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAAcgc                >  1:2049345/1‑148 (MQ=255)
                                                                                                                                      gCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAAcgcg               >  1:2697283/1‑148 (MQ=255)
                                                                                                                                         tCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGctc            >  1:1555004/1‑148 (MQ=255)
                                                                                                                                         tCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGctc            <  1:1210412/148‑1 (MQ=255)
                                                                                                                                          cGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGctct           >  2:479478/1‑148 (MQ=255)
                                                                                                                                            cAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGa         >  2:2404222/1‑148 (MQ=255)
                                                                                                                                                gggTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAAt     <  1:1266006/148‑1 (MQ=255)
                                                                                                                                                gggTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAAt     <  1:3242155/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGA  >  CP009273/2296718‑2297012

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: