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,4090CT85.0% 411.7 / 65.0 194intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (17/11);  new base T (63/99);  total (83/111)
Fisher's exact test for biased strand distribution p-value = 3.84e-02
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 194 total aligned reads displayed.

CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263263‑3263554
                                                                                                                                                  |                                                                                                                                                 
         ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCt                                                                                                                                        <  2:1523912/149‑1 (MQ=37)
            gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCa                                                                                                                                     <  2:839287/149‑1 (MQ=37)
               tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCaaaa                                                                                                                                  <  1:941083/149‑1 (MQ=37)
                 gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATTTTTGTGGCGCTCACAAAGc                                                                                                                                >  2:3198325/1‑149 (MQ=11)
                 gTGGAAAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGc                                                                                                                                <  1:4425190/149‑1 (MQ=255)
                    gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaa                                                                                                                             >  1:323106/1‑149 (MQ=37)
                    gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaa                                                                                                                             <  1:3014641/149‑1 (MQ=37)
                      cagtcagtCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaaa                                                                                                                           <  2:2983144/149‑1 (MQ=255)
                        gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTATTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                         >  2:2749229/1‑149 (MQ=255)
                         tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                        >  1:4308941/1‑149 (MQ=37)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccg                                                                                                                       <  1:2958324/149‑1 (MQ=37)
                           attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATTCCCCGCCAAAAGCAAAAACccgc                                                                                                                      >  1:3761874/1‑149 (MQ=255)
                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                                                                                   >  2:2219027/1‑149 (MQ=21)
                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTCTAGAGTCGTGGACAGTCCTTCATCTTTCTGCCCCGCCAAAAGCACAAACCCGCCGa                                                                                                                   >  2:2377949/1‑149 (MQ=11)
                               taatttCTGCCACTCAAAAATTAAAAACCCGCCGTAGCGTTTTTTTACGTAAAACATGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCAAAAAGCAAAAACCCGCCGaa                                                                                                                  <  1:418855/146‑1 (MQ=11)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                  <  1:128304/149‑1 (MQ=37)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                  >  2:964418/1‑149 (MQ=37)
                                   ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                              <  1:154602/149‑1 (MQ=37)
                                   ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTGATTTTTCTTGCCCTCCAAAAGCAAAATACCGACGACGCgg                                                                                                              >  2:4513188/1‑149 (MQ=11)
                                   ttCGGCCCCTCCACACTTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGTTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCgg                                                                                                              <  2:2266205/149‑1 (MQ=11)
                                       gCCCCTCAAAAATAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                          <  2:527659/149‑1 (MQ=11)
                                        ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATATTGTCGTGGACCGGCATTGATTTTTATGCCCCTCCCAAAGCAATAACCCGCCGAAGCGGGtttt                                                                                                         >  1:2050501/1‑149 (MQ=11)
                                         cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                        <  1:2677728/149‑1 (MQ=37)
                                         cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGTGTCGTGGACAGTCATTCCTGTTGCGGCCCCTCCAAAAGCAAAAACCCGCCGCAGCgggggttt                                                                                                        >  2:3694434/1‑144 (MQ=11)
                                          ccTCCACAAGTAAACACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                       <  2:3721962/149‑1 (MQ=255)
                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCAATCAATTTTATTCCCCTCCAAAAGAAAAAAACCGCCGAAGCGGGTTTTTAAGTaaa                                                                                                 >  1:1645561/1‑149 (MQ=11)
                                                  aGTAAAAACCCGCCGCAGCGGTTTTTTACGTACAACAGGTGAAACTGACCGAGAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                               <  2:1618201/149‑1 (MQ=11)
                                                     aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                                                                            <  1:4015061/149‑1 (MQ=255)
                                                      aaatCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                           <  2:3621770/149‑1 (MQ=37)
                                                      aaacCCCGCCCTAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                           <  2:1982524/149‑1 (MQ=255)
                                                      aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                           >  1:65825/1‑149 (MQ=17)
                                                       aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg                                                                                          <  2:4285980/149‑1 (MQ=37)
                                                          cccGCCGCACCGGGTTCTGACGTACAACCGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                       <  2:4032822/149‑1 (MQ=11)
                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                       <  2:4385346/149‑1 (MQ=37)
                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                       <  1:2794522/149‑1 (MQ=37)
                                                             gccgAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                                                    <  2:2508220/149‑1 (MQ=255)
                                                             gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                                                                    <  2:4121260/149‑1 (MQ=37)
                                                              ccgCAGCGGGTTTTTACGTACAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                   <  2:1800899/149‑1 (MQ=255)
                                                                gAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                 <  1:2093959/149‑1 (MQ=37)
                                                                 ttgCGGTTTTTTACGTATATCTGGTTAATGTGACCGATTAGCCTCTTTCTTTTGTGTATAGTGTCGTGGACAGTAATTCATTTTTCTGCCCCTCCAAAATAAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCg                                                                                <  1:2433330/147‑1 (MQ=11)
                                                                  aGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGCAAAAAACCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGa                                                                               >  1:3319373/1‑149 (MQ=11)
                                                                   gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAt                                                                              <  2:920678/149‑1 (MQ=255)
                                                                    cGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                             <  1:1999150/149‑1 (MQ=255)
                                                                    cGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                             <  1:2192136/149‑1 (MQ=37)
                                                                       gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                          <  2:1845381/149‑1 (MQ=37)
                                                                        tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                                         <  1:3310552/149‑1 (MQ=37)
                                                                        tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAACCCGGTGAAACTGACCGAAAAGcc                                                                         >  2:3332057/1‑149 (MQ=255)
                                                                         ttttACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                        <  2:2618663/149‑1 (MQ=255)
                                                                         ttttACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                        <  2:2618689/149‑1 (MQ=255)
                                                                          tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAAACCCCCGCAGCCGGTTTTTACGTAAACCAGGTTAAAATGACCGATAAACCGc                                                                       >  1:3599219/1‑149 (MQ=11)
                                                                          tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                       <  1:4184132/149‑1 (MQ=38)
                                                                          tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCTCTCCAAAAGCAAAAACCCGCCGACGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGAAAAGCCGc                                                                       >  2:4080018/1‑149 (MQ=255)
                                                                            tACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                     <  2:242676/149‑1 (MQ=255)
                                                                             aCGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcttt                                                                    <  2:1250159/149‑1 (MQ=255)
                                                                              cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                   <  2:4329333/149‑1 (MQ=255)
                                                                               gTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                  <  1:2351213/149‑1 (MQ=255)
                                                                                 acaaCAGGTCAAACTGCCCGATACGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACTGTCCTTCAATTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                                <  2:555190/147‑1 (MQ=11)
                                                                                 aaaaCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                                <  1:3831615/149‑1 (MQ=255)
                                                                                 aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                                >  2:339294/1‑149 (MQ=255)
                                                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                               <  1:1319675/149‑1 (MQ=255)
                                                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTTAGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTAtttt                                                               >  1:72589/1‑149 (MQ=14)
                                                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                               >  1:771760/1‑149 (MQ=255)
                                                                                   atcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                              <  2:4409917/147‑1 (MQ=2)
                                                                                   aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                              <  1:4580239/149‑1 (MQ=38)
                                                                                   aaCAGGTGAAACTGACAGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGTTTCTTTTg                                                              >  2:2290499/1‑149 (MQ=255)
                                                                                     cAGGTGAAACGGACCGATAAGCCGCTTTCTTTTGGGTCTAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                            <  2:1106535/149‑1 (MQ=255)
                                                                                      aGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                                           <  2:4674268/149‑1 (MQ=21)
                                                                                      aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTgggg                                                           >  2:1819699/1‑148 (MQ=255)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATCGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                       <  1:3828983/149‑1 (MQ=25)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                       <  1:223649/149‑1 (MQ=25)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                       <  2:2314779/149‑1 (MQ=25)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                       <  1:223653/149‑1 (MQ=25)
                                                                                            aaCTGACCGATAAGCCGCTTTCTTTTGGGTCTAGTGTCGACGAAAGACTATCTAATTTCTGCCACTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                     <  2:678371/149‑1 (MQ=11)
                                                                                              cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGGAAAGcgtc                                                   >  2:3168058/1‑149 (MQ=255)
                                                                                              cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                   >  2:4540432/1‑149 (MQ=255)
                                                                                               tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcg                                                  >  2:1840239/1‑149 (MQ=37)
                                                                                                gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                 <  1:20918/149‑1 (MQ=37)
                                                                                                    gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGCGTCGAGGAc                                             >  2:4348523/1‑149 (MQ=255)
                                                                                                     aTACGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCACTGTTCTGCCCCTCCAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTACACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                            <  2:4206496/149‑1 (MQ=11)
                                                                                                     aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGGACa                                            >  2:1375874/1‑149 (MQ=255)
                                                                                                     aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAATGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGaga                                            >  2:4776566/1‑147 (MQ=19)
                                                                                                       aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTATTGGACAGt                                          >  1:1121860/1‑149 (MQ=11)
                                                                                                          ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAAATTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAACGATAATCCGCTTTCTTTTGGGTATAGCGTCGTGGTCAGtcca                                       >  1:3504524/1‑147 (MQ=11)
                                                                                                          ccGCTGTCTTTTGGGTCTCGTGTCGTGGACAGTCATTCCTTGTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTACACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                       <  2:3980035/149‑1 (MQ=11)
                                                                                                            gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtccatc                                     >  1:2073713/1‑145 (MQ=255)
                                                                                                             ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGCCAggccatca                                    >  2:871860/1‑142 (MQ=11)
                                                                                                               ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATc                                  <  1:3914190/149‑1 (MQ=37)
                                                                                                                  ttttTGGAATACGGTCTTGGTAAGTCAGTCCTTTTTCGGGCCCTCCACAATCTAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCgtt                               <  2:1227594/149‑1 (MQ=11)
                                                                                                                             gtgtCGTGGACAGTCCTTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAGTAATATTTCAGCCccacca                    >  1:1097910/1‑149 (MQ=11)
                                                                                                                              tgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                   <  1:1767539/149‑1 (MQ=37)
                                                                                                                              tgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCaa                   >  1:2697399/1‑149 (MQ=255)
                                                                                                                                  gTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCTTGGACAGTCATTCATCTTTCTGCACCTCCAAAAGc               >  2:3052251/1‑149 (MQ=255)
                                                                                                                                  gTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaatc               >  1:772227/1‑147 (MQ=37)
                                                                                                                                    ggACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCaa             >  1:103071/1‑149 (MQ=255)
                                                                                                                                     gACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGCaaa            >  1:2615989/1‑149 (MQ=37)
                                                                                                                                      aCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAATCaaaa           >  1:4720891/1‑149 (MQ=255)
                                                                                                                                        aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGTGGAAAGTCATTCATAATTAAGCCCCCACAAAAGCAAAAAc         >  1:4717376/1‑149 (MQ=11)
                                                                                                                                         gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc        >  1:1044942/1‑149 (MQ=37)
                                                                                                                                            tttcatTTTTCTGCCCCTCAAAAAGGTAAAACCAGCCGACGCGGGTTTTTACGTTAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc     <  1:1511538/148‑1 (MQ=11)
                                                                                                                                              tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACCGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg   >  2:2657858/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263263‑3263554

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: