Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 776,694 T→A 5.6% intergenic (+16/‑131) lysY → / → lysZ tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273776,6940TA 5.6% 392.9 / 2.6 147intergenic (+16/‑131)lysY/lysZtRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base T (63/66);  new base A (8/2);  total (78/69)
Fisher's exact test for biased strand distribution p-value = 9.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

ACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTT  >  CP009273/776549‑776841
                                                                                                                                                 |                                                                                                                                                   
 ccTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                 <  2:4332615/149‑1 (MQ=255)
 ccTCCCTTACATGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGAATTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                 <  1:823428/149‑1 (MQ=255)
 ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAACCTGCACGACCCACCAGTTTTAAAATCAAAAGCAGa                                                                                                                                                 >  1:1523078/1‑149 (MQ=255)
   tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATg                                                                                                                                               <  2:2619799/149‑1 (MQ=255)
       ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTaa                                                                                                                                           <  2:409209/149‑1 (MQ=255)
       ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTaa                                                                                                                                           <  2:409039/149‑1 (MQ=255)
        tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAg                                                                                                                                          <  2:3087339/149‑1 (MQ=255)
        tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGATCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTAGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATAAAACTAAGAAGTTAAg                                                                                                                                          >  1:192469/1‑149 (MQ=14)
        tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTAGTTATCTCAGTTGGTAGAGCTTTTTACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAg                                                                                                                                          <  1:708198/149‑1 (MQ=14)
         aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                         >  1:1017511/1‑149 (MQ=255)
         aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                         >  2:3055884/1‑149 (MQ=255)
         aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGGTTTAACAGAAAACTCAGATGTTAAGc                                                                                                                                         >  1:41539/1‑149 (MQ=255)
           aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGt                                                                                                                                       >  2:2933891/1‑149 (MQ=255)
              gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                    <  2:1832917/149‑1 (MQ=255)
              gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                    <  2:2301520/149‑1 (MQ=255)
               aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                   <  1:3296574/149‑1 (MQ=255)
               aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                   >  2:3949489/1‑149 (MQ=255)
                 ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGa                                                                                                                                 <  2:3692905/149‑1 (MQ=255)
                    gTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGTGAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                              <  1:1259199/149‑1 (MQ=255)
                    gTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATaa                                                                                                                              >  1:4247786/1‑149 (MQ=255)
                     tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGTCTTTTAATCAATTGGTAGCAGGGTCGAATACTGCACGACCCACCAGTTTTTACAAAAAACACAGCTTTTTAGCGGGAAGGATGAc                                                                                                                             >  1:388284/1‑149 (MQ=14)
                          ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgc                                                                                                                        <  2:339099/149‑1 (MQ=255)
                          ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgc                                                                                                                        <  2:4022453/149‑1 (MQ=255)
                            ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACGTGAAAAGAAGATGTTAACCGGGACGGATGACGTTgcgc                                                                                                                      >  2:1140301/1‑149 (MQ=14)
                             tCGATACCGACAACACCCACCACTCGGGTCGTTATGTGTGTTGGTATTGATGTCTTATTTTACACATTACGATGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                     <  1:2390202/149‑1 (MQ=9)
                              cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCa                                                                                                                    <  1:176897/149‑1 (MQ=255)
                              cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTAGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGGGAAGGATAACGTTGCGCCa                                                                                                                    >  2:3976939/1‑149 (MQ=255)
                               gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACCTTGCGCCAg                                                                                                                   >  1:508666/1‑149 (MQ=255)
                                aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCTTGAAGGATAACGTTGCTCCAGc                                                                                                                  >  2:1773562/1‑149 (MQ=255)
                                aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                  >  1:1357115/1‑149 (MQ=255)
                                    cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAACGTTGCGCCCGCAACg                                                                                                              >  2:939879/1‑149 (MQ=255)
                                     gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGTTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgg                                                                                                             >  1:3497780/1‑149 (MQ=255)
                                        atcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTACTTTGTTACGCGTGAAGGAAAACGGTGCGCCAGCAACGGccc                                                                                                          >  1:727731/1‑149 (MQ=14)
                                           acccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGATAACGTTTCGCCAGCAACCGCCCgga                                                                                                       >  1:3325592/1‑147 (MQ=255)
                                                   cTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCGGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgagc                                                                                               <  2:1016807/149‑1 (MQ=255)
                                                       ggTCGTTAGCTCAGTTTGTAGAGCACTTTACTTTTAATAAAATGGTCGCAGGTTCGAATACAGCACGACCCACCACTTTTAAAATAAAACACAGATGTTAATCGTGAAGGATAACGTTGCGCCACCAAAGGCCCGTAGGGCGagcgaag                                                                                           >  1:2710984/1‑149 (MQ=255)
                                                             ttgctatgttgttAGAGCAGTTTACTTTTAATATATTGGGCGCAGTTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCATTTTTTAAGCTTGAATGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                     <  2:1231212/142‑1 (MQ=14)
                                                                 tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTAAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTATGGCGAGCGACGCGAGTCATcc                                                                                 >  2:3254201/1‑149 (MQ=255)
                                                                    gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTGGGGCGAGCGAAGCGAGTCATCCTGc                                                                              >  2:1819577/1‑149 (MQ=255)
                                                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAATCCCAGAAGAGAAGCGTAACGGAAAACATATCGCCAGAACCTGCCCGTGGGACGAGCGCAGCGAGTCATGCTGAACgtc                                                                         >  2:2463800/1‑147 (MQ=9)
                                                                           gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCTGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCACCGGCCCGGAGGGCGCGCGAAGCGAGTCATCCTGCACGAccc                                                                       >  1:1610052/1‑149 (MQ=14)
                                                                            agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                                                                      >  2:870399/1‑149 (MQ=255)
                                                                                gTTGACTTTTAATCATTGGGTCGCTGGTTCTAAGCCTGCACGACCCACCATTTTAAACATCAAACTCAGATTTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                                  <  2:3256723/149‑1 (MQ=14)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                     >  1:4091385/1‑82 (MQ=255)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                     >  1:4091407/1‑82 (MQ=255)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                     <  2:4091385/82‑1 (MQ=255)
                                                                                  tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATc                                                                <  2:2710984/149‑1 (MQ=255)
                                                                                   ttcTTTTAATATTTTGGTCGCTGTTTGTACTCCTGCACGACCCACCTGTTTTAACATAAAACGACGTTTTTATGCGTGATGGAAATCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                               <  1:3113728/147‑1 (MQ=9)
                                                                                    aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTGAAGGATAACTTTGCGCCAGAAACGGCCCGTAGGGCGAGCGTAGCGAGTCATCCTGCACGACCCACCAATCtt                                                              >  2:3216323/1‑149 (MQ=255)
                                                                                     cTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCATTTTTAACATCAAACTCAGATGTTAAGCTTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTa                                                             <  2:2186167/149‑1 (MQ=255)
                                                                                         taataaTTTGGTCGCTGGTTCGAATCCTGCACGACCCACCATTTTTAACATCAAACTCAGATGTTAAGCGTTATGGATAACGTTGCGCCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGa                                                         <  1:3184274/149‑1 (MQ=17)
                                                                                          aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                        >  1:2953068/1‑149 (MQ=255)
                                                                                             catttGGTGGCTGGTTCGATTCGTGCACGCCCCACATTTTTTAACATAAAACTCATATTTTAAGGGTTATGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                     <  1:1942551/146‑1 (MQ=17)
                                                                                             cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                     >  2:4180659/1‑149 (MQ=255)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                    >  2:823478/1‑149 (MQ=255)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                    >  2:2985141/1‑149 (MQ=255)
                                                                                                ttGGTCGCAGTGTCGACTCCGGCACGACCCACCAGTTTTAACATCATACTCTGATGTTATGCGTGAAGTATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGccc                                                  <  1:3883072/149‑1 (MQ=17)
                                                                                                 ggcccccAGGTTCTAATCCCGCACCCCCCACCATTTTTAAAATCAAACTCAGATGTTAAGCGTTAAGTATAACTTTGCGCCACAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                 <  2:1796777/143‑1 (MQ=17)
                                                                                                  ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAACGAttagatcgg                                                >  1:2309926/1‑142 (MQ=255)
                                                                                                  ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCg                                                >  1:3251919/1‑149 (MQ=255)
                                                                                                   gtcgatggTTCTAATCCTGCACGACCCACAATTTTTAAAATAAAACTATGTTTTTATGGGTGATGGATACCTTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                               <  1:938665/143‑1 (MQ=18)
                                                                                                     cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                             >  2:1303662/1‑149 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAGTCAGa                                                                                                                                                 >  2:3402605/1‑46 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                 <  1:3402605/46‑1 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                              >  2:3439749/1‑65 (MQ=255)
                                                                                                        aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                              <  1:3439749/65‑1 (MQ=255)
                                                                                                         ggTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                         <  1:3216909/53‑1 (MQ=255)
                                                                                                         ggTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                         >  2:3216909/1‑53 (MQ=255)
                                                                                                          gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCATATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaa                                        <  2:3034408/149‑1 (MQ=255)
                                                                                                          gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTaaaaa                                        >  1:1266905/1‑149 (MQ=255)
                                                                                                             cGATTCCTGCACGACCCACCAGTTTTAACATCAACCTCAGTTGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                     <  2:776170/149‑1 (MQ=255)
                                                                                                              gAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                    <  1:3932195/149‑1 (MQ=255)
                                                                                                              gAATCCGGCACGCCCGACGCGTTTGCACTTCCCACTCAGTTGTTATGGGTGTAGGATACCGTTGCGCCAGCACCGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                    <  2:1119417/149‑1 (MQ=18)
                                                                                                               aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                   >  1:1210193/1‑149 (MQ=255)
                                                                                                               aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                   >  1:1074526/1‑149 (MQ=255)
                                                                                                                aTCCTGCACGACCCACCAGTTTTAACCTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                  >  1:1845378/1‑149 (MQ=255)
                                                                                                                aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                  >  1:1845420/1‑149 (MQ=255)
                                                                                                                aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                  >  2:1166249/1‑149 (MQ=255)
                                                                                                                 tCCTGCACGACCCACCAGTTTTAACATCAAACTCATTTGTTATGCGTGATGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                 <  2:1455794/149‑1 (MQ=21)
                                                                                                                    tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCTAGTAAAAATCTTTCAgga                              >  1:2767192/1‑148 (MQ=255)
                                                                                                                    tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                              >  1:2094027/1‑149 (MQ=255)
                                                                                                                    tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                              >  2:410046/1‑149 (MQ=255)
                                                                                                                    ggCACGACCCACCATTTTTAACAACAAACTCAGATTTTAAGGGTTAAGGATAAGGTTGCGCCCGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                              <  1:885398/148‑1 (MQ=21)
                                                                                                                     gaaCGACCCACCAGTTTTAAAATAAAAATCTGTTTTTATGCTTTATGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                             <  1:4145725/147‑1 (MQ=21)
                                                                                                                     gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                             <  2:2305460/149‑1 (MQ=255)
                                                                                                                     gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                             <  2:147579/149‑1 (MQ=255)
                                                                                                                         gACCCACCAGTTTTAACATCAAACTCAGATTTTAAGCGTTAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                         <  1:1044037/149‑1 (MQ=255)
                                                                                                                           cccACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAccc                       >  1:2602435/1‑149 (MQ=255)
                                                                                                                              acctGTTTTAACATCAACCTCAGATGTTATGCGTGATGGATAAGGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                    <  1:2076010/149‑1 (MQ=25)
                                                                                                                                 aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                 <  2:1017511/149‑1 (MQ=255)
                                                                                                                                  gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                >  1:3494398/1‑149 (MQ=255)
                                                                                                                                  gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAACTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                >  1:2365234/1‑149 (MQ=255)
                                                                                                                                   ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt               >  1:2314251/1‑149 (MQ=255)
                                                                                                                                    tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTc              <  2:1355730/149‑1 (MQ=255)
                                                                                                                                       aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGtt           >  1:2103585/1‑149 (MQ=255)
                                                                                                                                            aaaaCACAGTTGTTAAGCGTTAAGGATACCGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc      <  1:2933891/148‑1 (MQ=25)
                                                                                                                                             aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTCTGGGTCGTTAGCTCa     >  1:1183372/1‑149 (MQ=255)
                                                                                                                                             aaaCACAGTTGTTACGCGTGAAGGATACCGTTGCGCCACCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa     <  1:585667/149‑1 (MQ=37)
                                                                                                                                              aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAg    >  1:3624895/1‑149 (MQ=255)
                                                                                                                                                cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGtt  >  1:4317615/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTT  >  CP009273/776549‑776841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: