Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 31,153 G→T 7.3% V113F (GTC→TTC)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927331,1530GT 7.3% 1078.2 / 37.1 286V113F (GTC→TTC) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (127/138);  new base T (13/8);  total (140/146)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 292 total aligned reads displayed.

AACCTACATCGAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGT  >  CP009273/31005‑31301
                                                                                                                                                    |                                                                                                                                                    
 aCCTACATCGAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGgtgt                                                                                                                                                     >  1:4535608/1‑149 (MQ=255)
  ccTACATCGAGCCGATTCACTGGGAAGTTGTACGAAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTc                                                                                                                                                    <  1:3686073/149‑1 (MQ=255)
    tACATCGAGCCGATTCACTGGTAAGTTGTACGAAATATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTTTCAc                                                                                                                                                  <  1:3834717/149‑1 (MQ=255)
      cATCGAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCa                                                                                                                                                >  2:16250/1‑149 (MQ=255)
        tCGAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATg                                                                                                                                              <  1:3445451/149‑1 (MQ=255)
          gAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGAt                                                                                                                                            <  2:3956601/149‑1 (MQ=255)
             ccGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTgg                                                                                                                                         >  1:4279382/1‑149 (MQ=255)
              cGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGt                                                                                                                                        <  2:693000/149‑1 (MQ=255)
              cGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGt                                                                                                                                        >  1:3196362/1‑149 (MQ=255)
                aTTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGc                                                                                                                                      >  1:2953708/1‑149 (MQ=255)
                   cACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCAc                                                                                                                                   <  2:3204577/149‑1 (MQ=255)
                   cACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCAc                                                                                                                                   >  1:3350025/1‑149 (MQ=255)
                   cACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCAc                                                                                                                                   >  1:2646470/1‑149 (MQ=255)
                         gAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGa                                                                                                                             >  1:1869500/1‑149 (MQ=255)
                           aGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATg                                                                                                                           >  1:741501/1‑149 (MQ=255)
                            gTTGTCCGCAAGAGTAGTGAAACAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGCACGTCAGGGCGTGTTGTAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGc                                                                                                                          <  2:4472169/149‑1 (MQ=255)
                            gTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGc                                                                                                                          <  2:2422629/149‑1 (MQ=255)
                            gTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGc                                                                                                                          <  2:354794/149‑1 (MQ=255)
                            gTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGGGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGc                                                                                                                          <  2:29934/149‑1 (MQ=255)
                             ttGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCg                                                                                                                         >  1:2121308/1‑149 (MQ=255)
                                tACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGAtt                                                                                                                      <  1:374248/149‑1 (MQ=255)
                                 aCGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTg                                                                                                                     <  2:2808580/149‑1 (MQ=255)
                                    cAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATa                                                                                                                  <  2:40502/149‑1 (MQ=255)
                                    cAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATa                                                                                                                  <  2:2271543/149‑1 (MQ=255)
                                      aGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATaaa                                                                                                                <  1:1272205/149‑1 (MQ=255)
                                      aGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATaaa                                                                                                                <  1:2493490/149‑1 (MQ=255)
                                      aGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATaaa                                                                                                                <  2:1948440/149‑1 (MQ=255)
                                       gATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAg                                                                                                               >  2:1837992/1‑149 (MQ=255)
                                           attGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCaga                                                                                                           >  1:2837610/1‑149 (MQ=255)
                                           attGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCaga                                                                                                           <  2:101339/149‑1 (MQ=255)
                                            ttGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCagaa                                                                                                          >  1:4462195/1‑149 (MQ=255)
                                              gAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCagaaga                                                                                                        <  2:4284901/149‑1 (MQ=255)
                                              gAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCagaaga                                                                                                        <  2:903517/149‑1 (MQ=255)
                                               aaaaaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAc                                                                                                       <  2:1478564/149‑1 (MQ=255)
                                               aaaaaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAc                                                                                                       <  1:4086938/149‑1 (MQ=255)
                                               aaaaaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAc                                                                                                       >  2:4361444/1‑149 (MQ=255)
                                                aaaaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAcc                                                                                                      >  2:4443042/1‑149 (MQ=255)
                                                aaaaGAGCGCCCGGACGCGGTGCTGCCAACGAGGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTTTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAcc                                                                                                      <  2:2634492/149‑1 (MQ=255)
                                                 aaaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccg                                                                                                     >  1:2168003/1‑149 (MQ=255)
                                                  aaGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgc                                                                                                    <  2:4598747/149‑1 (MQ=255)
                                                   agagCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgcc                                                                                                   >  2:1238925/1‑149 (MQ=255)
                                                   agagCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgcc                                                                                                   >  2:921997/1‑149 (MQ=255)
                                                   agagCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgcc                                                                                                   >  1:3591960/1‑149 (MQ=255)
                                                    gagCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgccg                                                                                                  >  2:4720786/1‑149 (MQ=255)
                                                    gagCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGAccgccg                                                                                                  >  2:1972386/1‑149 (MQ=255)
                                                      gcgcCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCcgtc                                                                                                >  1:694575/1‑149 (MQ=255)
                                                           cGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTc                                                                                           <  2:1154972/149‑1 (MQ=255)
                                                           cGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTc                                                                                           <  2:1399516/149‑1 (MQ=255)
                                                             gACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGa                                                                                         <  2:504532/149‑1 (MQ=255)
                                                                gcgGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGt                                                                                      >  1:2983364/1‑149 (MQ=255)
                                                                gcgGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGt                                                                                      >  1:2729824/1‑149 (MQ=255)
                                                                   gTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGc                                                                                   <  1:4325716/149‑1 (MQ=255)
                                                                   gTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGc                                                                                   >  2:3052445/1‑149 (MQ=255)
                                                                    tgctgcCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCg                                                                                  >  2:3593904/1‑149 (MQ=255)
                                                                     gctgcCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGa                                                                                 <  1:644152/149‑1 (MQ=255)
                                                                       tgcCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTAGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATg                                                                               <  2:374986/149‑1 (MQ=255)
                                                                          cAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATgaag                                                                            <  1:248045/149‑1 (MQ=255)
                                                                           aaCGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATgaaga                                                                           <  1:4063474/149‑1 (MQ=255)
                                                                             cGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGaaa                                                                         <  1:2051561/149‑1 (MQ=255)
                                                                                 gggCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTg                                                                     <  2:2177774/149‑1 (MQ=255)
                                                                                     ggTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCt                                                                 <  1:310946/149‑1 (MQ=255)
                                                                                      gTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTGTCGACGTAGCGATGAAGAAAATTGGTCTg                                                                >  1:1346091/1‑149 (MQ=255)
                                                                                        cAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGa                                                              <  2:3740202/149‑1 (MQ=255)
                                                                                          gACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaaa                                                            >  2:1291972/1‑149 (MQ=255)
                                                                                           aCGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAAc                                                           <  1:3961412/149‑1 (MQ=255)
                                                                                                gcTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACcgcgc                                                      <  1:3462001/149‑1 (MQ=255)
                                                                                                gcTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACcgcgc                                                      >  1:1241253/1‑149 (MQ=255)
                                                                                                    aaCTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTc                                                  >  2:4507969/1‑149 (MQ=255)
                                                                                                    aaCTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTc                                                  <  2:3059278/149‑1 (MQ=255)
                                                                                                     aCTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTcc                                                 <  2:77951/149‑1 (MQ=255)
                                                                                                     aCTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTcc                                                 >  1:1094397/1‑149 (MQ=255)
                                                                                                     aCTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTcc                                                 <  2:2037300/149‑1 (MQ=255)
                                                                                                              tggagctggaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGca                                        >  1:2244313/1‑149 (MQ=255)
                                                                                                                 agctggaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGcacac                                     >  1:4176331/1‑149 (MQ=255)
                                                                                                                   ctggaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCacacac                                   <  2:1059574/149‑1 (MQ=255)
                                                                                                                   ctggaACGTCAGGGCGTGTTGGAAGAGGTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCacacac                                   <  2:3447855/149‑1 (MQ=255)
                                                                                                                    tggaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACg                                  >  2:4111486/1‑149 (MQ=255)
                                                                                                                    tggaACGTCAGGGCGGGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACg                                  >  2:232535/1‑149 (MQ=255)
                                                                                                                     ggaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGa                                 >  2:2172925/1‑149 (MQ=255)
                                                                                                                      gaACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGGTCCGGTATCGCACACACGAt                                <  2:266115/149‑1 (MQ=255)
                                                                                                                       aaCGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATg                               >  1:1064379/1‑149 (MQ=255)
                                                                                                                             aGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGgaagaa                         >  1:4640546/1‑149 (MQ=255)
                                                                                                                             aGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGgaagaa                         <  1:943779/149‑1 (MQ=255)
                                                                                                                              gggCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGaagaag                        <  2:1161981/149‑1 (MQ=255)
                                                                                                                              gggCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGaagaag                        <  1:1818020/149‑1 (MQ=255)
                                                                                                                                gCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAgcg                      >  2:2906030/1‑149 (MQ=255)
                                                                                                                                 cGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAgcgc                     <  2:4279382/149‑1 (MQ=255)
                                                                                                                                 cGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAgcgc                     <  1:1303321/149‑1 (MQ=255)
                                                                                                                                    gtTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTgg                  >  1:3287599/1‑149 (MQ=255)
                                                                                                                                     ttGGAAGAGTTCGGTTTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGc                 >  1:4187890/1‑149 (MQ=255)
                                                                                                                                         aaGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGACGCGCTGGCGGtt             >  2:2398234/1‑149 (MQ=255)
                                                                                                                                           gagTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGc           >  1:1184553/1‑149 (MQ=255)
                                                                                                                                             gTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCg         >  2:2449468/1‑149 (MQ=255)
                                                                                                                                              ttCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGc        >  2:2814149/1‑149 (MQ=255)
                                                                                                                                              ttCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGc        >  2:3887408/1‑149 (MQ=255)
                                                                                                                                              ttCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATAGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGc        >  1:1609668/1‑149 (MQ=255)
                                                                                                                                               tCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCt       <  1:1436508/149‑1 (MQ=255)
                                                                                                                                                cGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTg      <  1:2181721/149‑1 (MQ=255)
                                                                                                                                                 ggTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGa     >  1:4781481/1‑149 (MQ=255)
                                                                                                                                                    gtCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGt  <  1:171657/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCTACATCGAGCCGATTCACTGGGAAGTTGTACGCAAGATTATTGAAAAAGAGCGCCCGGACGCGGTGCTGCCAACGATGGGCGGTCAGACGGCGCTGAACTGCGCGCTGGAGCTGGAACGTCAGGGCGTGTTGGAAGAGTTCGGTGTCACCATGATTGGTGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACGT  >  CP009273/31005‑31301

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: ‑