Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,315 T→C 6.8% intergenic (+354/+260) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,3150TC 6.8% 656.6 / 33.0 254intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (53/172);  new base C (7/9);  total (70/184)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 257 total aligned reads displayed.

GCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  CP009273/3263167‑3263462
                                                                                                                                                    |                                                                                                                                                   
gccgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                                                                                                     <  2:2362732/149‑1 (MQ=255)
   gcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATATTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                  >  1:2724066/1‑149 (MQ=255)
      tCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                               <  1:2035713/149‑1 (MQ=255)
      tCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCAACTTTCTGCCCCTCCCAAAGGAAAAAc                                                                                                                                               >  1:3307571/1‑149 (MQ=255)
       ccGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAcc                                                                                                                                              <  1:2779695/149‑1 (MQ=255)
       ccGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAcc                                                                                                                                              <  1:3853377/149‑1 (MQ=255)
        cGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAccc                                                                                                                                             <  2:2388878/149‑1 (MQ=255)
        cGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAccc                                                                                                                                             <  2:188245/149‑1 (MQ=255)
          gCATGAATAAAGCTCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTTGGTATAGTGTCGTGGCCCCTCATTAATATTTCTGCCCCCACAAAAAAAAAAAAccgc                                                                                                                                           >  1:1227594/1‑149 (MQ=21)
          gCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                                                                                                                                           <  1:3109735/149‑1 (MQ=255)
           cATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgcc                                                                                                                                          <  1:3591136/149‑1 (MQ=255)
             tGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGa                                                                                                                                        >  2:2093892/1‑149 (MQ=255)
               aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                      <  1:4041409/149‑1 (MQ=255)
                aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                     <  1:4678799/149‑1 (MQ=255)
                aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                     <  2:1714206/149‑1 (MQ=255)
                   aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                  <  1:3838659/149‑1 (MQ=255)
                    aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                                                                                                 <  2:1284926/149‑1 (MQ=255)
                        cACTCAACAATCTGAACCCCGCCGGAACGGTTTATTACGTACCTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                             <  2:2119797/149‑1 (MQ=255)
                        cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTGTGGGGGAGAGTGTATAGGACAGTCAGTCATCTTTCTGCCCCACCAAATTAAAATTCCCGGCGTAGCGGgttgtt                                                                                                                             >  2:81124/1‑149 (MQ=18)
                         aCTCAACAATCTGAACCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  2:522279/149‑1 (MQ=255)
                         aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                            <  2:3056844/149‑1 (MQ=255)
                          cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                           <  1:603164/149‑1 (MQ=255)
                           tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                          <  2:4746238/149‑1 (MQ=255)
                             aacacTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                        <  2:4039708/149‑1 (MQ=255)
                                aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                     <  1:3269388/149‑1 (MQ=255)
                                aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                     <  1:4387835/149‑1 (MQ=255)
                                 aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                    >  1:4154193/1‑149 (MQ=255)
                                 aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAAGTacaaa                                                                                                                    >  1:2662868/1‑148 (MQ=255)
                                  tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                   >  1:714013/1‑149 (MQ=255)
                                  tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                   >  2:3828974/1‑149 (MQ=255)
                                     gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATGTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                >  2:1064002/1‑149 (MQ=255)
                                     gAAACCCGCCGGAACGGTTTATTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                                                <  2:1163315/149‑1 (MQ=11)
                                        aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                             >  1:4206496/1‑149 (MQ=255)
                                          ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                           <  1:489222/149‑1 (MQ=255)
                                              cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                       <  1:4147539/149‑1 (MQ=255)
                                                gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGtttt                                                                                                                              <  2:2740843/124‑1 (MQ=255)
                                                 aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                    >  1:4340035/1‑149 (MQ=255)
                                                  aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                   <  2:1994947/149‑1 (MQ=255)
                                                  aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTACAAAATTAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAAAATGTGAAACTGAACGAt                                                                                                   >  1:4290473/1‑149 (MQ=14)
                                                   cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                  <  1:2749276/149‑1 (MQ=255)
                                                   cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATa                                                                                                  <  2:740158/149‑1 (MQ=255)
                                                    ggTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                 <  1:4437826/149‑1 (MQ=255)
                                                     gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                >  2:2516489/1‑149 (MQ=255)
                                                      tttATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                               <  1:1697205/149‑1 (MQ=255)
                                                      tttATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGGGAAACTGACCGATAAGc                                                                                               >  1:1982524/1‑149 (MQ=255)
                                                      tttATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaatc                                                                                                                    >  2:2569953/1‑126 (MQ=255)
                                                       ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                              <  1:2551218/149‑1 (MQ=255)
                                                       ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                                                                                              <  1:1031069/149‑1 (MQ=255)
                                                         attaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                                                            >  2:4658453/1‑148 (MQ=255)
                                                           taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt                                                                                          <  2:3834041/149‑1 (MQ=255)
                                                           taCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                                          <  1:3786294/149‑1 (MQ=17)
                                                             cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc                                                                                        <  2:1309961/149‑1 (MQ=255)
                                                               tACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                      <  1:1970787/149‑1 (MQ=40)
                                                               tACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCCAAAGTAAAAACCCCCCGAAGCGGGTGTTTACGTAAAACAGGTGAAACGGACCGATAAGCCGctttctt                                                                                      >  1:920678/1‑149 (MQ=14)
                                                                aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACCAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                     <  2:3327845/149‑1 (MQ=255)
                                                                 cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                    <  1:3053549/149‑1 (MQ=39)
                                                                   tCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  2:4487828/149‑1 (MQ=38)
                                                                   tCAGGTAAAACTGAACGATATGCCGCTTTCTTTTGCGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCACAACTAAAAACCCGCCGAAGCGGGTTTTTACGTTTAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  1:2611914/149‑1 (MQ=11)
                                                                        tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                             <  1:4375733/149‑1 (MQ=35)
                                                                         aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAg                                                                            >  2:389512/1‑149 (MQ=25)
                                                                         aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                            <  1:3756952/149‑1 (MQ=25)
                                                                          aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt                                                                           <  2:355619/149‑1 (MQ=25)
                                                                            aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                         <  1:2884380/149‑1 (MQ=37)
                                                                             cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc                                                                        <  2:3237591/149‑1 (MQ=37)
                                                                               gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACATTCAATCATCTTTCTGCCCCTCCAAAATGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGtatatggtagt                                                                      >  2:2846741/1‑142 (MQ=11)
                                                                                    ttAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                 <  2:3708514/148‑1 (MQ=37)
                                                                                    aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACAGATAATCCGCTTTCTTTTGGGTAAAGTGTCGAGGACa                                                                 >  2:2630370/1‑149 (MQ=255)
                                                                                      aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGt                                                               <  2:3031933/149‑1 (MQ=37)
                                                                                      aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGt                                                               <  1:1221900/149‑1 (MQ=37)
                                                                                      aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGt                                                               <  2:2782302/149‑1 (MQ=37)
                                                                                         ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGTTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcat                                                            <  1:4408495/149‑1 (MQ=37)
                                                                                          cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTTCCCCAACAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAAAGGTTAAACTGACCGATAAGCCGCATTCTTTTGGGTATAGTGTCGTTGACAAtcaat                                                           >  1:2281863/1‑147 (MQ=11)
                                                                                          cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATATTTATGCCCCTCCAAAATTAAAAAACCCCCCAAACGGGTTTTTAAGTAAAAAAGGTGAAACAGAACGATAATCCCCTTTCTTTTTGGTATAGTGTCGTAGAAAgaaaaa                                                           >  1:1845381/1‑144 (MQ=11)
                                                                                           gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                          <  2:2480172/149‑1 (MQ=37)
                                                                                           gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                          <  1:2110442/149‑1 (MQ=37)
                                                                                             tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcat                                                        <  2:4732413/149‑1 (MQ=37)
                                                                                             tttctttTGGGTATAGTGTCGTGGACACTCAATCATCATTATGCCCCTCCAAAAGTAAAAACCCGCCCAAGCGGGTTTTTAAGTAAAACAGGTGAAAAAGACCGAAAATCCCCTTTATTTTGGGTATAGTTTCGTGGTCAAtcaataaa                                                        >  1:4535046/1‑144 (MQ=11)
                                                                                                ctttTGGGTATCGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCtt                                                     <  1:299826/149‑1 (MQ=37)
                                                                                                ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCtt                                                     <  1:2692789/149‑1 (MQ=37)
                                                                                                         ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCt                                            <  2:1523912/149‑1 (MQ=37)
                                                                                                            gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCa                                         <  2:839287/149‑1 (MQ=37)
                                                                                                               tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCaaaa                                      <  1:941083/149‑1 (MQ=37)
                                                                                                                    gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCGCCTCCAAAAGCaaa                                 >  1:323135/1‑149 (MQ=255)
                                                                                                                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccg                           <  2:4081550/149‑1 (MQ=37)
                                                                                                                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccg                           <  1:2958324/149‑1 (MQ=37)
                                                                                                                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTCTAGAGTCGTGGACAGTCCTTCATCTTTCTGCCCCGCCAAAAGCACAAACCCGCCGa                       >  2:2377949/1‑149 (MQ=11)
                                                                                                                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                      <  1:128304/149‑1 (MQ=37)
                                                                                                                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                      >  2:964418/1‑149 (MQ=37)
                                                                                                                                tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                     <  2:4562170/149‑1 (MQ=37)
                                                                                                                                  tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                   >  1:1071430/1‑149 (MQ=37)
                                                                                                                                   ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                  <  2:4290473/149‑1 (MQ=37)
                                                                                                                                   ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                  <  1:154530/149‑1 (MQ=37)
                                                                                                                                   ttCGGCCCCTCCACACTTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGTTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCgg                  <  2:2266205/149‑1 (MQ=11)
                                                                                                                                       gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt              <  2:527660/149‑1 (MQ=37)
                                                                                                                                        ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATATTGTCGTGGACCGGCATTGATTTTTATGCCCCTCCCAAAGCAATAACCCGCCGAAGCGGGtttt             >  1:2050501/1‑149 (MQ=11)
                                                                                                                                          ccTCCACAAGTAAACACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa           <  2:3721962/149‑1 (MQ=255)
                                                                                                                                           cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCCTTTTTCTGCCCCTGCAAAAGCAAAATCCCGCCGAAGCGGGTTTTGAc          >  1:3511426/1‑149 (MQ=11)
                                                                                                                                              cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa       <  2:1656668/149‑1 (MQ=37)
                                                                                                                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCAATCAATTTTATTCCCCTCCAAAAGAAAAAAACCGCCGAAGCGGGTTTTTAAGTaaa     >  1:1645561/1‑149 (MQ=11)
                                                                                                                                                   gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTGGGGGAGAGTGTCGAGGCCAGTCCTTCAATATTCTGCCCCTCCAAGAGCAAAAACCCGCCGCAGCGGGGTGGTACGGAAAACa  >  2:2603551/1‑149 (MQ=11)
                                                                                                                                                    |                                                                                                                                                   
GCCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  CP009273/3263167‑3263462

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