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,4090CT91.2% 266.7 / 19.2 120intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (50/56);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 123 total aligned reads displayed.

TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263261‑3263554
                                                                                                                                                    |                                                                                                                                                 
ttctttTGGGTATAGTGTCGTGCACAGTCATACAAATTTCCGCCCCTCCAAAAGTAAAAACCCGCAGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATc                                                                                                                                                   <  2:2560813/149‑1 (MQ=11)
  ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAtttt                                                                                                                                                 <  1:2967995/149‑1 (MQ=37)
   ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                                <  2:3066388/149‑1 (MQ=37)
     ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCt                                                                                                                                              <  1:1393499/149‑1 (MQ=37)
        ggTATTGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGcc                                                                                                                                           <  2:1797877/149‑1 (MQ=255)
        ggTATAGTGTCGTGGAAAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGcc                                                                                                                                           <  2:1797857/149‑1 (MQ=255)
            taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGTTTTTACGTAAAAAAAGTGAAACTGACCGATAACCCGCTTTCCTTTGGGTATAGTGTCGTGGGCAGTCATTAATTTTTCTGccccgg                                                                                                                                       >  1:1688266/1‑147 (MQ=11)
            taGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                       <  2:2586354/149‑1 (MQ=255)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTcc                                                                                                                                      <  2:324509/149‑1 (MQ=255)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTcc                                                                                                                                      <  2:1475994/149‑1 (MQ=37)
               tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAGAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTGTGGACCTCCaa                                                                                                                                    >  1:3857891/1‑149 (MQ=11)
                gtCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGATTAATTTTTCTGCCCCTCCcaa                                                                                                                                   >  1:2654568/1‑149 (MQ=11)
                 tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTTTTGGGTATAGTGTCGTGGAGAGGCATTCATTGTTGTGCCGCtgcgagt                                                                                                                                  >  2:3557950/1‑143 (MQ=11)
                  cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTTGACAGTCATTAATTTTTCTGCCCCTCCAAAAg                                                                                                                                 <  2:1775792/149‑1 (MQ=11)
                      gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa                                                                                                                             <  1:839525/149‑1 (MQ=37)
                      gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                               >  1:272189/1‑147 (MQ=17)
                      gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCa                                                                                                                               <  2:272189/147‑1 (MQ=17)
                           tcattcatCTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                        <  2:1861836/149‑1 (MQ=255)
                               tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                                                                                                    >  2:3063113/1‑149 (MQ=17)
                                  tCTTTCTGCCCCTCCAAAAGTAAAACCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                                 <  1:724401/149‑1 (MQ=255)
                                      tCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                             <  2:1226095/149‑1 (MQ=37)
                                        tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                           <  1:3157951/149‑1 (MQ=37)
                                        tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGTGTAGAGTGTCGTGGACCGTCGTTGATGTTTCGGCCGCTCCAAATGCACATTCCCTCCCACGCGGGtg                                                                                                           >  2:2972764/1‑148 (MQ=11)
                                         gacccTCCAAAAGAAAAAACCCGCCGTAGCGTGTTTTTACGTAAAAAAGGTGAAACTGACCTATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                          <  1:3392760/147‑1 (MQ=11)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                          <  1:1801956/149‑1 (MQ=37)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                          <  1:3392963/149‑1 (MQ=37)
                                            ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                       <  1:81778/149‑1 (MQ=37)
                                             cTCCAAAAGGACAAACCCGCCGCAGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTCAc                                                                                                      <  2:3956835/149‑1 (MQ=11)
                                                  aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTGCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                 >  2:129715/1‑149 (MQ=255)
                                                    aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                               <  2:3504300/149‑1 (MQ=255)
                                                      tAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                             <  2:382401/149‑1 (MQ=255)
                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGGCCGTCATTGATATTTCTGGCCCTCCAAAAGCAAAAAACCGCCGGAGCGGGGTTTTAGGTAAAACAggg                                                                                           >  2:884368/1‑148 (MQ=11)
                                                             ccgccgAAGCGGGTTTTTACGTAACACAGGTGAAACTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc                                                                                      <  2:3273354/149‑1 (MQ=255)
                                                             ccgccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc                                                                                      <  2:1027901/149‑1 (MQ=255)
                                                              cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                     >  1:2396576/1‑149 (MQ=17)
                                                                 cgcaGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                  <  2:2274179/146‑1 (MQ=255)
                                                                  gacgCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                 <  2:3268452/146‑1 (MQ=255)
                                                                   aaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                >  1:1650347/1‑149 (MQ=255)
                                                                    aGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                                                               <  1:2102970/149‑1 (MQ=37)
                                                                      cGGCTTATTACTTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                             <  2:34115/149‑1 (MQ=11)
                                                                        ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTTAACCAGGTGAAACTGACCGATAAg                                                                           >  2:2890447/1‑149 (MQ=255)
                                                                         tttcttACGTACAACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                          <  2:3122466/145‑1 (MQ=255)
                                                                         gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                          <  2:3551732/149‑1 (MQ=37)
                                                                         gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                          <  1:3474301/149‑1 (MQ=37)
                                                                           ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTAATTTTTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACATACCGATAAGCCg                                                                        >  2:1424090/1‑149 (MQ=11)
                                                                              tACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                     >  1:2000234/1‑149 (MQ=255)
                                                                              tACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATTCCCCTCCAAAAGAAAAAAACCCCCGCAGCGGGTTTTTACGTAAATCAGGGGAAAATGGCCGATAAGCCgggt                                                                     >  1:3558633/1‑146 (MQ=1)
                                                                                cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                   >  2:3320026/1‑149 (MQ=255)
                                                                                cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTAATTTGGGTATAGTGTCGTGGACACTCATTCATTTAACAGCCCCACCAAAACCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGACCGAAAAAccccgtaa                                                                   >  2:2890179/1‑143 (MQ=11)
                                                                                 ttAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCTGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                                                  <  1:129715/148‑1 (MQ=14)
                                                                                   aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                                <  1:2905917/149‑1 (MQ=255)
                                                                                   aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACCTAAACCCGGGGAAACTGCCCGATAAGCCGcgtgagtg                                                                >  2:2083619/1‑142 (MQ=11)
                                                                                      aCAGTTTAAACTGACCGATAAGCCGCTTTCTTTTTTTTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                             <  1:501912/149‑1 (MQ=11)
                                                                                        aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGt                                                           >  2:650635/1‑149 (MQ=255)
                                                                                        aGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCATTCATTTTTGTGCCCCACCAAAACCAAAAACACGCCCAATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGTTTATCCGCTTGATTTTgggg                                                           >  1:2274943/1‑148 (MQ=11)
                                                                                          gtaaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                        >  2:2797014/4‑102 (MQ=1)
                                                                                           tcaaaCTGAGATATAAGCCGCTTACTGTTGGCTATAGTGACGTTCCCAGTCATTCATTTTTATCCCACTCCACAAGCAAAAACCCGCCGCATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                        <  2:2014152/147‑1 (MQ=11)
                                                                                           tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                        <  1:4059491/149‑1 (MQ=25)
                                                                                            gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                       <  1:140420/149‑1 (MQ=25)
                                                                                              aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTggggacagtc                                                     >  2:550178/1‑142 (MQ=1)
                                                                                              aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTGTGGGTATAGcg                                                     >  2:3198021/1‑149 (MQ=9)
                                                                                               aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                    <  2:1097275/149‑1 (MQ=37)
                                                                                               aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCAGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTTGGGATAGcgt                                                    >  2:1613318/1‑149 (MQ=11)
                                                                                                  gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                     >  2:944130/1‑97 (MQ=18)
                                                                                                     ctacaaGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                              <  2:2654568/145‑1 (MQ=255)
                                                                                                     cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                              <  2:2654595/149‑1 (MQ=37)
                                                                                                      gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCAATCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                             >  2:918537/1‑149 (MQ=255)
                                                                                                      gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                             <  1:618537/149‑1 (MQ=37)
                                                                                                       aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTAGTGGACa                                            >  2:197499/1‑149 (MQ=255)
                                                                                                       aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAATGACCGATAAGCCGCGTTCTTTAGGGTAAAGCGTAGGGGaaa                                            >  1:2010371/1‑147 (MQ=11)
                                                                                                        tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCAATAATTTTTATGCCCCTCCAAAATCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAAACAGGTGAAACAGACCGCTAAGCCGCTTTCTTTTGGGGATAGCGACGTggcccg                                           >  1:1647449/1‑145 (MQ=11)
                                                                                                         aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCCTCGTGGACAgg                                          >  2:1290933/1‑148 (MQ=255)
                                                                                                          aGCCGGTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                         <  1:700055/149‑1 (MQ=255)
                                                                                                          aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGCGTCGTGGAAAGtc                                         >  1:3775993/1‑149 (MQ=255)
                                                                                                          aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTAAGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCCATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                         >  2:2689679/1‑149 (MQ=255)
                                                                                                          aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                         <  1:700068/149‑1 (MQ=37)
                                                                                                          aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                         <  2:2550509/149‑1 (MQ=255)
                                                                                                          aGCCGCTTTCTTTTGGGTATAGTGACGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                         <  1:700144/149‑1 (MQ=12)
                                                                                                             cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTATTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGCGACGTGGACAgacctt                                      >  2:2598563/1‑144 (MQ=11)
                                                                                                             cGCTATCTTTTGGGTATTGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcatt                                      <  2:3857891/149‑1 (MQ=255)
                                                                                                              gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcaatc                                     >  2:1503670/1‑149 (MQ=255)
                                                                                                               ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                       >  2:412910/1‑82 (MQ=33)
                                                                                                                tttctttTGTGTATAGTGTCGTGGACAGTCATTCAATTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattcat                                   <  2:940231/149‑1 (MQ=255)
                                                                                                                tttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattcat                                   <  2:940937/149‑1 (MQ=255)
                                                                                                                tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCCATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtccataat                                   >  1:972032/1‑143 (MQ=11)
                                                                                                                 ttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTGGTGGACAGtcaatcata                                  >  2:1031399/1‑148 (MQ=255)
                                                                                                                    ttttGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                               <  1:4058014/149‑1 (MQ=255)
                                                                                                                     ttagggTTTCGTGTCGTCGACAGTCTATCTTTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTc                              <  2:2177627/146‑1 (MQ=11)
                                                                                                                          gTATCGTGTCGCGGACAGTCATTCATTTTTCGGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                         <  2:1616040/149‑1 (MQ=11)
                                                                                                                          gTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGGGGGCAGTGCTTCATCGTTGTGccc                         >  1:1920985/1‑149 (MQ=11)
                                                                                                                           tataGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTATTcccc                        >  2:1019467/1‑149 (MQ=11)
                                                                                                                               gtgtCGTGGACAGTCATTCATTTTTCAGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCCATAAGCCGCGTTCTTTTGGGGATAGCGTCGTGGAAAGGCAATCAACTTTATGCCCCGCCa                    >  1:1302279/1‑149 (MQ=11)
                                                                                                                                         cAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaaa          >  1:2399937/1‑149 (MQ=37)
                                                                                                                                           gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGCAAAAAcc        >  2:3115446/1‑149 (MQ=255)
                                                                                                                                           gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc        >  1:2070754/1‑149 (MQ=37)
                                                                                                                                            tcattcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc       <  1:853760/149‑1 (MQ=37)
                                                                                                                                              attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc     <  1:1439795/149‑1 (MQ=37)
                                                                                                                                              attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc     >  2:519539/1‑149 (MQ=255)
                                                                                                                                              attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTCGGGATAGCGGCGTGGACAGGCAATCATATTACTTCCCCACCATAAGAAAAAACccga     >  2:2838569/1‑148 (MQ=11)
                                                                                                                                                 catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa  >  1:2826521/1‑149 (MQ=37)
                                                                                                                                                    |                                                                                                                                                 
TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGA  >  CP009273/3263261‑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: