Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,409 C→T 100% intergenic (+448/+166) 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,4090CT90.7% 229.2 / 16.0 109intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base T (43/53);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 111 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552
                                                                                                                                                |                                                                                                                                               
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                              <  1:749842/148‑1 (MQ=37)
 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                             <  1:1207133/148‑1 (MQ=37)
   gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTg                                                                                                                                           <  1:2139031/148‑1 (MQ=37)
   gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTg                                                                                                                                           <  1:2138977/148‑1 (MQ=37)
    gggataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGc                                                                                                                                          <  2:3452917/145‑1 (MQ=255)
     gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGcc                                                                                                                                         <  2:1526375/148‑1 (MQ=37)
      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGccc                                                                                                                                        <  1:443323/148‑1 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                     <  1:845757/148‑1 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                                     >  2:2044908/1‑148 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                                     >  2:2044993/1‑148 (MQ=37)
               gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAg                                                                                                                               <  2:2074020/148‑1 (MQ=37)
               gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAg                                                                                                                               <  2:2073933/148‑1 (MQ=37)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                             >  2:2337754/1‑148 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                          <  2:468290/148‑1 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                          <  2:2244780/148‑1 (MQ=37)
                     aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCATTTGGGTATAGTGTAGTGGAAAGTCATTAATTTTTCTTCCCCTCcaaaaacaaaaa                                                                                                                         >  2:2102984/1‑148 (MQ=11)
                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAAcc                                                                                                                       >  1:2068598/1‑148 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                      <  2:2805671/148‑1 (MQ=37)
                        cagtcatCTTTCTGCCCCTCCAAACGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                      <  2:1871892/145‑1 (MQ=255)
                            catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATGTTTGTGGCCCTCCAAAAGCAAAAACcagccg                                                                                                                  >  2:668417/1‑148 (MQ=255)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                <  2:285938/148‑1 (MQ=37)
                                    tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                          <  2:825150/148‑1 (MQ=37)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                         <  1:596045/148‑1 (MQ=37)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                         <  1:2384855/148‑1 (MQ=37)
                                      ccccTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  1:3371255/148‑1 (MQ=12)
                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  1:3371262/148‑1 (MQ=37)
                                       cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                       >  1:278796/1‑148 (MQ=37)
                                         cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGGCCCTGCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTa                                                                                                     >  1:1958685/1‑148 (MQ=11)
                                             aaaaGAACAAACCCGCCGCAGCGGGTTATTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                 <  2:1818358/148‑1 (MQ=11)
                                              taaGTAACACCCCGCCGAAGCGGGTTTTTACGTAACACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTCTAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  2:2096905/147‑1 (MQ=11)
                                              aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  1:416275/148‑1 (MQ=37)
                                              aaaCTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  1:416395/148‑1 (MQ=25)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                               >  2:2953017/1‑148 (MQ=255)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                               <  2:2821523/148‑1 (MQ=37)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                          >  2:12479/1‑148 (MQ=12)
                                                     aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                         <  2:1513895/148‑1 (MQ=37)
                                                     aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCTTGGCCAGTCCGTCATTGTTATGCCCCGCCCAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTcaaacgggg                                                                                         >  2:137441/1‑142 (MQ=11)
                                                       aCCCGCCGAAGCGGGTTTTTACGTACAACCGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                                                                       <  1:448094/148‑1 (MQ=255)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCCAGCAATTTGCTCCCCCACCAAAAGCAGAATCCCGACGTAGCGGGGTTTTCAGTAAACAAGGTgt                                                                                       >  2:1408458/1‑147 (MQ=11)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                                                                       <  2:278796/148‑1 (MQ=17)
                                                             cttatcGGTTTTTTAAGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                 <  2:841431/143‑1 (MQ=11)
                                                                 gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGa                                                                             >  2:2109857/1‑148 (MQ=255)
                                                                    ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTATGCCCCGCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATaa                                                                          >  1:1315007/1‑148 (MQ=11)
                                                                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTTAACCAGGTGAAACTGCCCGCTAAg                                                                         >  2:2885713/1‑148 (MQ=255)
                                                                      tttttACGTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                        <  2:151460/148‑1 (MQ=255)
                                                                      tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                        >  1:240304/1‑148 (MQ=12)
                                                                         ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                     <  2:3273678/148‑1 (MQ=38)
                                                                           aCGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                   >  2:1741323/1‑148 (MQ=255)
                                                                            cGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcttt                                                                  <  2:711337/148‑1 (MQ=255)
                                                                              tacatcAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                <  2:3183250/143‑1 (MQ=12)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcgttct                                                                >  1:3179427/1‑148 (MQ=255)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                <  1:641266/148‑1 (MQ=255)
                                                                               aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                                               <  2:2153618/148‑1 (MQ=255)
                                                                                aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                              <  1:355037/148‑1 (MQ=255)
                                                                                aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                              <  1:2147413/148‑1 (MQ=255)
                                                                                 caCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:219659/147‑1 (MQ=255)
                                                                                 aaCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:1744765/148‑1 (MQ=255)
                                                                                 aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:220449/148‑1 (MQ=38)
                                                                                 aaCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:342198/148‑1 (MQ=16)
                                                                                   cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAGCAGGTTAAACTGACCGATAAGCCGGTTTCTTTTgg                                                           >  2:151414/1‑148 (MQ=255)
                                                                                     ggTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTgggg                                                         >  2:708967/1‑147 (MQ=255)
                                                                                       tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                       <  2:926729/148‑1 (MQ=25)
                                                                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                      >  2:3318317/1‑148 (MQ=2)
                                                                                         aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCTGGTTTTTACGTTAACCAGGTGAAACAGACCCATAAGCCGCGTTATTTTGGGtattg                                                     >  2:431780/1‑146 (MQ=11)
                                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                  >  2:3198406/1‑148 (MQ=37)
                                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                  <  1:1044662/148‑1 (MQ=37)
                                                                                                ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTAAATTTTGTTCCCCACCAAAAACAAAAAAACCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAACAGACCCATAATCCGCTTTGTTTTGGGTGTAGCGTCGTg                                              >  1:1235600/1‑148 (MQ=2)
                                                                                                 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                             <  1:703004/148‑1 (MQ=37)
                                                                                                   aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGagtcattgaa                                           >  2:2312499/1‑142 (MQ=11)
                                                                                                    tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTATGCCCCACCAAAATCAAAAAAACGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCCCTAAGCCGCTTTCTTTTGGGTATAGCGTCGTggtac                                          >  2:2703706/1‑145 (MQ=11)
                                                                                                     aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGAACGATAAGCCGCTTTCTTTTTGGTATAGCGTCGTGGCCAg                                         >  1:275841/1‑148 (MQ=11)
                                                                                                       gccgccTACTTTTGTGTACAGTGTGGAGGACAGACATTCAACCTTCAGCCCCCCCAACCGCACTAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                       <  2:761209/148‑1 (MQ=11)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                       <  1:2337754/148‑1 (MQ=37)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTATGCCCCTCCAAAATCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAAAAGCCGCTTTGTTTTGGGTATAGCGTCGTGGGCCGtca                                      >  2:2480091/1‑148 (MQ=11)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                      >  2:3189492/1‑148 (MQ=255)
                                                                                                         cGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAAAAGACCGATAACCCGCTTTCTTTTTGGTATATCGTCGGGGagaggtaa                                     >  1:3332236/1‑144 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAATGAACGATAAGCCGCGTTCTGTTGGGGATAGCGTATGGGACAGgcaat                                    >  2:3243189/1‑146 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTCACCCACCAAAAGCAAAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTGAAACTTACCGATAAGCCGCTTAGTTGTGGGTATAGCGTCGTGGCCAGtcctt                                    >  2:3242980/1‑148 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       <  1:1804563/65‑1 (MQ=35)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       >  2:1804563/1‑65 (MQ=35)
                                                                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTTGAAAGtcaatcaaaa                               >  1:2980481/1‑145 (MQ=11)
                                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                              <  1:2044993/148‑1 (MQ=37)
                                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                              <  1:2044908/148‑1 (MQ=37)
                                                                                                                 tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                             <  1:2109857/148‑1 (MQ=37)
                                                                                                                 tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                             <  1:2109847/148‑1 (MQ=255)
                                                                                                                     ggTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCtttaagc                         >  2:1955826/1‑144 (MQ=255)
                                                                                                                       tataGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAGCGTCGTGGCACGTCAATCATAATTCTGgcc                       >  2:2422766/1‑148 (MQ=11)
                                                                                                                        ataGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAACAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAATCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCAATCAAATTTGTGcccc                      >  1:60830/1‑148 (MQ=11)
                                                                                                                         gcgtgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCt                     <  1:331288/146‑1 (MQ=37)
                                                                                                                          aGTGTCGTGGACAGTCATTCATTTGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                    <  1:2441435/148‑1 (MQ=255)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa           >  1:883087/1‑148 (MQ=37)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATAATTCTGCCCCTCCAAAAGCaa           >  2:1150348/1‑148 (MQ=255)
                                                                                                                                    aCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa          >  1:3432054/1‑148 (MQ=37)
                                                                                                                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGCGTCGTGGAAAGTCATTCATCTTTATGCCCCTCCAAAAGCaaa          >  2:1974627/1‑148 (MQ=255)
                                                                                                                                      aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaac        >  2:2897636/1‑147 (MQ=37)
                                                                                                                                      aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTCATCTTTCTGCCCCGCCAAAAGCaaaaa        >  1:3179184/1‑148 (MQ=255)
                                                                                                                                         cattcattTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCCTCGTGGACAGTAATAAAACTTTATGCCCCGCAAAAAGCAAAAAacc     >  1:1178427/1‑148 (MQ=11)
                                                                                                                                           ttcatTTTTCTGCCCCTCCACACGCAACAACCCGCCGAAGCGGGTTTTTACGTCTTCCCGGTGACACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc   <  2:1672982/148‑1 (MQ=11)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc  >  2:290724/1‑148 (MQ=37)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATAAAATTTCTGCCCCACCAAAATCAAAAACccgcc  >  1:833932/1‑148 (MQ=11)
                                                                                                                                                |                                                                                                                                               
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552

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: