Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 32,659 C→A 93.5% R615S (CGC→AGC)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927332,6590CA93.5% 484.1 / 24.9 170R615S (CGC→AGC) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base A (71/88);  total (78/92)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 170 total aligned reads displayed.

GGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCG  >  CP009273/32515‑32803
                                                                                                                                                |                                                                                                                                                
ggcggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGcc                                                                                                                                               <  2:3334943/148‑1 (MQ=255)
ggcggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGcc                                                                                                                                               <  2:3334891/148‑1 (MQ=255)
 gcggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCct                                                                                                                                              <  2:1558851/148‑1 (MQ=255)
 gcggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCct                                                                                                                                              <  2:1559978/148‑1 (MQ=255)
 gcggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCct                                                                                                                                              <  1:1788299/148‑1 (MQ=255)
  cggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctc                                                                                                                                             >  2:887946/1‑148 (MQ=255)
  cggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctc                                                                                                                                             <  1:3008509/148‑1 (MQ=255)
  cggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctc                                                                                                                                             >  2:2333835/1‑148 (MQ=255)
   ggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctct                                                                                                                                            <  1:3426275/148‑1 (MQ=255)
    gcCCGAACCGTATCGGTCAGGGTCTCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTa                                                                                                                                           >  1:2655752/1‑148 (MQ=255)
      ccGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACt                                                                                                                                         <  1:3251241/148‑1 (MQ=255)
       cGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACtt                                                                                                                                        <  1:977666/148‑1 (MQ=255)
         aaCCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCg                                                                                                                                      <  2:2339895/148‑1 (MQ=255)
           ccGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAg                                                                                                                                    <  2:105918/148‑1 (MQ=255)
           ccGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGACTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAg                                                                                                                                    >  2:498063/1‑148 (MQ=255)
                tCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGt                                                                                                                               <  2:2689833/148‑1 (MQ=255)
                tCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGt                                                                                                                               >  2:1060067/1‑148 (MQ=255)
                  ggTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTaa                                                                                                                             <  1:2358945/148‑1 (MQ=255)
                   gTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAAc                                                                                                                            <  1:1143110/148‑1 (MQ=255)
                   gTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAAc                                                                                                                            >  2:839994/1‑148 (MQ=255)
                    tCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAAct                                                                                                                           <  2:2073575/148‑1 (MQ=255)
                    tCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAAct                                                                                                                           >  1:1419101/1‑148 (MQ=255)
                         gTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGa                                                                                                                      <  2:3254847/148‑1 (MQ=255)
                          gaTCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTACGACAGCCTCTACTTCGAGACGGTAACTCTGGaa                                                                                                                     >  2:1910318/2‑148 (MQ=255)
                            tCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGa                                                                                                                   <  1:2739844/148‑1 (MQ=255)
                             cGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGAt                                                                                                                  >  1:3225250/1‑148 (MQ=255)
                              gAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGAtg                                                                                                                 >  2:940373/1‑148 (MQ=255)
                                 ttCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGc                                                                                                              <  1:1319600/148‑1 (MQ=255)
                                  tCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCt                                                                                                             >  2:2996323/1‑148 (MQ=255)
                                  tCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCt                                                                                                             >  2:2995816/1‑148 (MQ=255)
                                      ctactGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGACGATGTGCTGGaa                                                                                                         >  2:2340961/1‑148 (MQ=255)
                                      ctactGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGaa                                                                                                         <  2:2586507/148‑1 (MQ=255)
                                          tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCg                                                                                                     >  2:643993/1‑148 (MQ=255)
                                          tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCg                                                                                                     <  1:1878420/148‑1 (MQ=255)
                                           gTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGt                                                                                                    >  2:1221963/1‑148 (MQ=255)
                                            ttGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTg                                                                                                   <  1:2832669/148‑1 (MQ=255)
                                            gtGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACGGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTg                                                                                                   <  2:1825199/147‑1 (MQ=255)
                                               cGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt                                                                                                <  1:1755836/148‑1 (MQ=255)
                                               cGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt                                                                                                <  1:728291/148‑1 (MQ=255)
                                                 tACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCGGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt                                                                                              >  1:2884330/1‑148 (MQ=255)
                                                 tACACGCCTCGCTGGCGCGGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt                                                                                              <  2:3052411/148‑1 (MQ=255)
                                                  acacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATc                                                                                             >  1:1026815/1‑148 (MQ=255)
                                                  acacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATc                                                                                             >  1:1853039/1‑148 (MQ=255)
                                                  acacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATc                                                                                             >  1:3119510/1‑148 (MQ=255)
                                                   cacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                            >  1:1072292/1‑148 (MQ=255)
                                                     cGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCgag                                                                                          <  2:2710882/148‑1 (MQ=255)
                                                        cTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAg                                                                                       <  1:1869479/148‑1 (MQ=255)
                                                         tCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                      >  1:2304023/1‑148 (MQ=255)
                                                         tCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                      <  1:1970590/148‑1 (MQ=255)
                                                         tCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                      <  1:1970086/148‑1 (MQ=255)
                                                              ggCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaaa                                                                                 <  1:3037415/148‑1 (MQ=255)
                                                               gcgcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAg                                                                                >  2:3159924/1‑148 (MQ=255)
                                                               gcgcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAg                                                                                <  1:519529/148‑1 (MQ=255)
                                                                cgcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAgg                                                                               <  2:1240420/148‑1 (MQ=255)
                                                                 gcTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGc                                                                              >  1:1068351/1‑148 (MQ=255)
                                                                    gcgcGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGtt                                                                           <  2:572672/148‑1 (MQ=255)
                                                                    gcgcGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGtt                                                                           <  1:794547/148‑1 (MQ=255)
                                                                    gcgcGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGtt                                                                           <  2:2562953/148‑1 (MQ=255)
                                                                     cgcgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTa                                                                          <  2:1260376/148‑1 (MQ=255)
                                                                           gACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTcc                                                                    >  1:1193677/1‑148 (MQ=255)
                                                                           gACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTcc                                                                    >  1:406226/1‑148 (MQ=255)
                                                                           gACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTcc                                                                    >  2:2224493/1‑148 (MQ=255)
                                                                                ttACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAc                                                               >  1:353124/1‑148 (MQ=255)
                                                                                  aCGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcgg                                                             <  2:2408959/148‑1 (MQ=255)
                                                                                     aaaCCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcgg                                                          <  1:621812/148‑1 (MQ=255)
                                                                                         cATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAg                                                      >  1:831996/1‑148 (MQ=255)
                                                                                           ttATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTCTCGTCCAGTACGGCGGTCAGAc                                                    >  2:2099022/1‑148 (MQ=255)
                                                                                               ggTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCg                                                <  2:3522215/148‑1 (MQ=255)
                                                                                               ggTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCg                                                >  1:2697428/1‑148 (MQ=255)
                                                                                               ggTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCg                                                <  2:2038166/148‑1 (MQ=255)
                                                                                                    aCTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                           <  2:1410870/148‑1 (MQ=255)
                                                                                                     cTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaaa                                          <  1:691453/148‑1 (MQ=255)
                                                                                                      tGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAAc                                         >  1:500707/1‑148 (MQ=255)
                                                                                                      tGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAAc                                         <  1:1451670/148‑1 (MQ=255)
                                                                                                       gTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACt                                        >  1:2633209/1‑148 (MQ=255)
                                                                                                        tAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTg                                       <  1:62327/148‑1 (MQ=255)
                                                                                                             cGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcg                                  <  2:3241720/148‑1 (MQ=255)
                                                                                                             cGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcg                                  <  2:1697694/148‑1 (MQ=255)
                                                                                                             cGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcg                                  <  2:195860/148‑1 (MQ=255)
                                                                                                               gAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcgcg                                >  1:2695656/1‑148 (MQ=255)
                                                                                                                 aaCCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcgcgcg                              >  2:2964282/1‑148 (MQ=255)
                                                                                                                 aaCCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcgcgcg                              >  2:692632/1‑148 (MQ=255)
                                                                                                                   ccGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCt                            >  2:798923/1‑148 (MQ=255)
                                                                                                                   ccGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCt                            <  2:2884330/148‑1 (MQ=255)
                                                                                                                     gTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTgg                          >  2:2030570/1‑148 (MQ=255)
                                                                                                                     gTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTgg                          <  1:2224493/148‑1 (MQ=255)
                                                                                                                      tctcCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGa                         >  2:828696/1‑148 (MQ=255)
                                                                                                                        tcCACCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAg                       <  1:2076927/148‑1 (MQ=255)
                                                                                                                           aCCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctg                    <  2:691064/148‑1 (MQ=255)
                                                                                                                           aCCGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctg                    >  2:2505052/1‑148 (MQ=255)
                                                                                                                            ccGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctgc                   <  1:1221963/148‑1 (MQ=255)
                                                                                                                             cGACTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctgct                  <  1:2189877/148‑1 (MQ=255)
                                                                                                                                cTACGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGc               >  1:90002/1‑148 (MQ=255)
                                                                                                                                  aCGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGt             >  1:1611865/1‑148 (MQ=255)
                                                                                                                                   cGACACTTCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTa            <  1:1364825/148‑1 (MQ=255)
                                                                                                                                       acgtCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCgg        <  1:498063/145‑1 (MQ=255)
                                                                                                                                         ttCCGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGtt      <  1:3534785/148‑1 (MQ=255)
                                                                                                                                           ccGACAGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTAt    <  2:272083/148‑1 (MQ=255)
                                                                                                                                             gACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGACACTGGCGCGCGCGCTGGAAGCTGCAGGCGTCCCGGTCATCg  >  2:1607227/1‑148 (MQ=255)
                                                                                                                                             gACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCg  >  2:3369334/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCG  >  CP009273/32515‑32803

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: