Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,146,030 T→A 9.1% S194C (AGT→TGT)  yegL ← VMA domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,146,0300TA 9.1% 1045.6 / 73.4 331S194C (AGT→TGT) yegLVMA domain protein
Reads supporting (aligned to +/- strand):  ref base T (130/171);  new base A (11/19);  total (141/190)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 336 total aligned reads displayed.

ATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACT  >  CP009273/2145882‑2146178
                                                                                                                                                    |                                                                                                                                                    
aTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACt                                                                                                                                                      <  1:675469/149‑1 (MQ=255)
 tCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTg                                                                                                                                                     <  2:606415/149‑1 (MQ=255)
 tCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTg                                                                                                                                                     >  2:2216484/1‑149 (MQ=255)
   cTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGa                                                                                                                                                   >  1:2300072/1‑149 (MQ=255)
    tGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGAt                                                                                                                                                  >  1:405606/1‑149 (MQ=255)
    tGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGAt                                                                                                                                                  >  1:1545428/1‑149 (MQ=255)
      aCAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATaa                                                                                                                                                >  2:882159/1‑149 (MQ=255)
       cAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAAc                                                                                                                                               >  1:1128902/1‑149 (MQ=255)
        aaGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAAcc                                                                                                                                              >  2:1207562/1‑149 (MQ=255)
        aaGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAAcc                                                                                                                                              <  2:1847326/149‑1 (MQ=255)
           gCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGc                                                                                                                                           <  2:3026405/149‑1 (MQ=255)
           gCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGc                                                                                                                                           <  1:3117490/149‑1 (MQ=255)
            cAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCt                                                                                                                                          >  1:884120/1‑149 (MQ=255)
             aaGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTa                                                                                                                                         >  1:1980775/1‑149 (MQ=255)
             aaGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTa                                                                                                                                         >  1:4481295/1‑149 (MQ=255)
              aGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTaa                                                                                                                                        <  1:4785803/149‑1 (MQ=255)
               gTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTaaa                                                                                                                                       >  2:1039380/1‑149 (MQ=255)
                tCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAAc                                                                                                                                      <  1:124565/149‑1 (MQ=255)
                tCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAAc                                                                                                                                      <  1:1292195/149‑1 (MQ=255)
                 ccGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACa                                                                                                                                     >  1:2918116/1‑149 (MQ=255)
                    cgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGtt                                                                                                                                  <  2:4626264/149‑1 (MQ=255)
                           gtgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGaa                                                                                                                           <  2:4307097/149‑1 (MQ=255)
                            tgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAAt                                                                                                                          <  1:4188856/149‑1 (MQ=255)
                            tgCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAAt                                                                                                                          <  1:3852902/149‑1 (MQ=255)
                             gCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGACtt                                                                                                                         >  2:4379373/1‑149 (MQ=255)
                             gCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAAtt                                                                                                                         <  1:3322877/149‑1 (MQ=255)
                                aCGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTa                                                                                                                      >  2:4229507/1‑149 (MQ=255)
                                   ttCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATc                                                                                                                   >  2:3846275/1‑149 (MQ=255)
                                   ttCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATc                                                                                                                   <  2:2955228/149‑1 (MQ=255)
                                   ttCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATc                                                                                                                   >  1:2498417/1‑149 (MQ=255)
                                           gTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCaga                                                                                                           <  2:3010119/149‑1 (MQ=255)
                                            ttGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCagag                                                                                                          <  2:2679764/149‑1 (MQ=255)
                                              gATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGc                                                                                                        <  1:1966847/149‑1 (MQ=255)
                                               aTGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCa                                                                                                       <  1:3134827/149‑1 (MQ=255)
                                                     aGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCt                                                                                                 >  1:4343093/1‑149 (MQ=255)
                                                     aGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCt                                                                                                 >  1:1651388/1‑149 (MQ=255)
                                                      gACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTg                                                                                                <  1:2996981/149‑1 (MQ=255)
                                                       aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGa                                                                                               >  2:813170/1‑149 (MQ=255)
                                                       aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGa                                                                                               <  2:4534777/149‑1 (MQ=255)
                                                       aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGa                                                                                               >  1:1394815/1‑149 (MQ=255)
                                                            aCGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAAc                                                                                          <  2:3722406/149‑1 (MQ=255)
                                                            aCGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAAc                                                                                          <  2:2356223/149‑1 (MQ=255)
                                                               ccAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCt                                                                                       <  1:684038/149‑1 (MQ=255)
                                                               ccAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCt                                                                                       >  1:3282668/1‑149 (MQ=255)
                                                               ccAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCt                                                                                       >  1:4581829/1‑149 (MQ=255)
                                                                cAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTa                                                                                      >  2:1589856/1‑149 (MQ=255)
                                                                 aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTaa                                                                                     <  1:4667051/149‑1 (MQ=255)
                                                                 aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTaa                                                                                     <  1:1612193/149‑1 (MQ=255)
                                                                 aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTaa                                                                                     >  1:2183948/1‑149 (MQ=255)
                                                                 aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACGTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCGTGCAGAGGCAAAGGCTGACGAACGCTaa                                                                                     <  2:1209355/149‑1 (MQ=255)
                                                                  gCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAt                                                                                    >  2:3407625/1‑149 (MQ=255)
                                                                  gCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAt                                                                                    <  2:1126420/149‑1 (MQ=255)
                                                                    tCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAttt                                                                                  <  2:1137343/149‑1 (MQ=255)
                                                                     cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTACTTtg                                                                                 >  2:2990680/1‑149 (MQ=255)
                                                                     cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtg                                                                                 <  1:575123/149‑1 (MQ=255)
                                                                     cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtg                                                                                 <  1:4248873/149‑1 (MQ=255)
                                                                      acacTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                >  1:1748074/1‑149 (MQ=255)
                                                                        acTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACGTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGACCGCTAATTTGTGc                                                                              <  2:4137370/149‑1 (MQ=255)
                                                                         cTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGcc                                                                             >  1:4555592/1‑149 (MQ=255)
                                                                          tGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCa                                                                            <  2:3920183/149‑1 (MQ=255)
                                                                          tGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCa                                                                            >  2:3788859/1‑149 (MQ=255)
                                                                           gACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAg                                                                           <  2:3451007/149‑1 (MQ=255)
                                                                            aCGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGc                                                                          <  1:2443291/149‑1 (MQ=255)
                                                                               tCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTc                                                                       <  2:2295542/149‑1 (MQ=255)
                                                                               tCCAGCCTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCt                                                                      <  1:4130729/149‑1 (MQ=255)
                                                                                   gCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCa                                                                   >  2:1485748/1‑149 (MQ=255)
                                                                                       tttGGCGCTTCCAGCACAACTTCCGGTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATc                                                               <  2:4255380/149‑1 (MQ=255)
                                                                                           gcgcTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGcc                                                           <  2:135397/149‑1 (MQ=255)
                                                                                            cgcTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGccc                                                          <  1:755988/149‑1 (MQ=255)
                                                                                             gcTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGcccc                                                         <  2:67746/149‑1 (MQ=255)
                                                                                              cTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCACAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGccccc                                                        >  2:2523775/1‑149 (MQ=255)
                                                                                              cTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGccccc                                                        <  1:1009862/149‑1 (MQ=255)
                                                                                               ttCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCt                                                       >  1:2663164/1‑149 (MQ=255)
                                                                                               ttCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCt                                                       <  1:3940264/149‑1 (MQ=255)
                                                                                                 ccAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGa                                                     >  2:2039565/1‑149 (MQ=255)
                                                                                                    gCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAaca                                                  >  1:2698832/1‑149 (MQ=255)
                                                                                                     cacaACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAacac                                                 >  1:460893/1‑149 (MQ=255)
                                                                                                      acaACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACAcc                                                <  2:3188395/149‑1 (MQ=255)
                                                                                                       caACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCg                                               >  2:358756/1‑149 (MQ=255)
                                                                                                          cTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATg                                            <  2:2174323/149‑1 (MQ=255)
                                                                                                            tCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGa                                          >  1:1911205/1‑149 (MQ=255)
                                                                                                                ttCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaag                                      <  2:369129/149‑1 (MQ=255)
                                                                                                                ttCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaag                                      >  2:4390614/1‑149 (MQ=255)
                                                                                                                  cccGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaa                                    >  2:2629827/1‑149 (MQ=255)
                                                                                                                     ggCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGc                                 <  2:750391/149‑1 (MQ=255)
                                                                                                                        gTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCaaa                              >  1:440983/1‑149 (MQ=255)
                                                                                                                            aGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCtt                          <  2:4676253/149‑1 (MQ=255)
                                                                                                                              cGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCtttt                        <  1:4751631/149‑1 (MQ=255)
                                                                                                                               gggAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTa                       <  1:2694909/149‑1 (MQ=255)
                                                                                                                                ggAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAt                      >  1:2748736/1‑149 (MQ=255)
                                                                                                                                 gAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtc                     >  1:234166/1‑149 (MQ=255)
                                                                                                                                 gAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtc                     <  2:3157895/149‑1 (MQ=255)
                                                                                                                                  aaaCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtct                    <  2:4083130/149‑1 (MQ=255)
                                                                                                                                    aCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttc                  >  2:4241108/1‑149 (MQ=255)
                                                                                                                                     ccGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttct                 <  1:358756/149‑1 (MQ=255)
                                                                                                                                              gCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGa        >  2:51966/1‑149 (MQ=255)
                                                                                                                                                gAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGaaa      >  1:3821843/1‑149 (MQ=255)
                                                                                                                                                 aaCTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAac     >  1:1402655/1‑149 (MQ=255)
                                                                                                                                                   cTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAacac   >  2:1300596/1‑149 (MQ=255)
                                                                                                                                                    tGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACt  <  1:1845597/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACT  >  CP009273/2145882‑2146178

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: