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.7% 393.2 / 21.0 183intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (5/10);  new base T (66/95);  total (77/106)
Fisher's exact test for biased strand distribution p-value = 7.84e-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 186 total aligned reads displayed.

TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGC  >  CP009273/3263264‑3263557
                                                                                                                                                  |                                                                                                                                                    
ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                                  <  2:1204751/150‑1 (MQ=37)
     ggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGccc                                                                                                                                             <  2:4587184/150‑1 (MQ=37)
        ataGTGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATTTTTATGccccga                                                                                                                                          >  1:1227242/1‑148 (MQ=255)
         taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTcc                                                                                                                                         <  1:2230126/150‑1 (MQ=37)
                gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                                  <  2:3119294/150‑1 (MQ=37)
                     cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAc                                                                                                                             <  1:1758424/150‑1 (MQ=37)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                           <  2:550818/150‑1 (MQ=37)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                           <  1:2462744/150‑1 (MQ=37)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGAGGACAGTCAATCATTATACTGCCCCTCCAAAATCAAAAAacc                                                                                                                           >  2:3400851/1‑150 (MQ=11)
                       gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTGGGGGATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCGAAAGCaaatcccc                                                                                                                           >  2:3400945/1‑145 (MQ=11)
                          attcatCTTTCTGCCCCCCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                                        <  2:2044752/150‑1 (MQ=255)
                            tcatCTTTCTGCCCCTCCAAATGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGAGAGTCAGTCATGTTTCTGCGCCTGCAAAAGCACAAACccgccgt                                                                                                                      >  2:5321633/1‑149 (MQ=11)
                            tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                                                                                      >  1:2906558/1‑150 (MQ=37)
                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGGGGTGGGTAATTCATTTTTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                                    >  2:1545627/1‑150 (MQ=11)
                                   tCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                               <  2:1332526/150‑1 (MQ=37)
                                    cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                              <  1:4127263/150‑1 (MQ=37)
                                     tggcacTCCACAAGTAAAAACCCGCCGAATCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                             <  1:3462271/145‑1 (MQ=11)
                                     tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTAGGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGgtt                                                                                                             >  2:4465532/1‑150 (MQ=11)
                                     tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATGTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                             >  2:966592/1‑150 (MQ=21)
                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                            <  2:1804635/150‑1 (MQ=37)
                                        cccACCAAAAGTACAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                          <  2:2091218/150‑1 (MQ=11)
                                         cactcAAAAATAACAAACCCGCCGTAGCGGTTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                          <  2:3095111/146‑1 (MQ=11)
                                            ccAACAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                      <  2:4937632/150‑1 (MQ=255)
                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                      >  1:2052100/1‑150 (MQ=255)
                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTATTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTTCGTa                                                                                                      >  2:1342171/1‑150 (MQ=11)
                                               aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                                   >  1:114847/1‑150 (MQ=38)
                                               aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTGTTACGTAAAc                                                                                                   >  1:4362137/1‑150 (MQ=255)
                                               aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCCTTCATTTTTCTGCCCCTCCGAAAGCAAAATCCCGCCGTCGCGGGTGTGTTCGTaata                                                                                                   >  2:4491954/1‑148 (MQ=11)
                                                aaGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                                  <  2:2006467/150‑1 (MQ=255)
                                                    tataaaCCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                                                                               <  2:637528/147‑1 (MQ=255)
                                                     aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg                                                                                             <  2:2264870/150‑1 (MQ=37)
                                                        aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                            >  2:1059233/1‑133 (MQ=34)
                                                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGAGGAGAGTCCTTAATTTTTCTGCCCCTACAAAAGCAAAAAACCGCCGTAGCGGGTTTTTACGTGAAGCAGGTGaaa                                                                                          >  2:1912398/1‑150 (MQ=11)
                                                          cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc                                                                                         <  2:2391169/150‑1 (MQ=17)
                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAAc                                                                                         >  1:4685616/1‑150 (MQ=17)
                                                           ccgccgTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt                                                                                        <  2:4424331/150‑1 (MQ=11)
                                                             gccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                      <  2:3589938/150‑1 (MQ=255)
                                                               cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAAACAGGTGAAACTGAcc                                                                                    >  1:2444328/1‑150 (MQ=255)
                                                                gtaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCg                                                                                   <  2:4458464/148‑1 (MQ=37)
                                                                gAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCg                                                                                   <  2:2834421/150‑1 (MQ=255)
                                                                gAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATATTTCTGCCCCGCCAAAAGCAAAAAACCGCCGAAGCGGGTATTTACGTACAACAGGTGAAACAGACCg                                                                                   >  2:3653739/1‑150 (MQ=11)
                                                                    cGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTTAACCAGGTGAAACTGTCCGTTaa                                                                               >  1:3717655/1‑150 (MQ=255)
                                                                    cGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATaa                                                                               <  1:2151233/150‑1 (MQ=37)
                                                                     gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGGACGATAAg                                                                              >  1:1743738/1‑150 (MQ=255)
                                                                      ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTTGAAAGTCATTCATGTTTCTGCCCCACCAAAAGAAAAAACACGCCGAAGCGGGTTTTTACGTAAAAAAGGTGAAACAGACCGAGaaac                                                                             >  1:3225662/1‑148 (MQ=11)
                                                                          tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGct                                                                         <  2:500412/150‑1 (MQ=38)
                                                                           ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgggt                                                                        >  1:4790165/1‑147 (MQ=0)
                                                                           ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                        >  1:5209668/1‑150 (MQ=255)
                                                                             aCGTACAACATTTGACACTGAGCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCTTTCTATTTACTGCCACTCCACAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                      <  1:1514987/150‑1 (MQ=11)
                                                                              cGTAAAACAGGTTCAACTGACCGATAAGCCGCTGTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                     <  2:4620214/150‑1 (MQ=255)
                                                                              cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctgtct                                                                     >  1:2085735/1‑150 (MQ=255)
                                                                              cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                     <  1:557859/150‑1 (MQ=255)
                                                                               gTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                                                    <  1:1118949/150‑1 (MQ=255)
                                                                                 aaaaCAGGTGAAACTGACCGCTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                  <  2:5022235/150‑1 (MQ=255)
                                                                                 aaaaCAGGTGAAACTGACCGATCAGCCGCTTTCTTTTGGCTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCGCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtatt                                                                  <  2:253390/150‑1 (MQ=14)
                                                                                 aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                  <  1:3889614/150‑1 (MQ=255)
                                                                                  aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCACTCCTTAATTTTTATGCCCCACCCAAAGCAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGTCCCATAAGCCGCGGTCTTTTg                                                                 >  1:4983138/1‑150 (MQ=11)
                                                                                   aaCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                <  2:1868485/150‑1 (MQ=25)
                                                                                      aGGTGAAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTGTggggga                                                             >  2:4837137/1‑148 (MQ=11)
                                                                                       ggTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                            <  1:1042840/150‑1 (MQ=17)
                                                                                         tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                          <  1:4197941/150‑1 (MQ=25)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc                                                         <  2:2906558/150‑1 (MQ=25)
                                                                                          gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                                                                                                          <  1:1032817/85‑1 (MQ=33)
                                                                                             aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                      <  1:4991505/150‑1 (MQ=37)
                                                                                                gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTg                                                   <  1:2473824/150‑1 (MQ=37)
                                                                                                 aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                                  <  1:840029/150‑1 (MQ=37)
                                                                                                   cGATAAGCCGCTTTCTTTTGGGTATCGTGTCGTGGACAGTCATACATTTTTCGGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                                <  2:3524396/150‑1 (MQ=11)
                                                                                                   cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                                <  1:4528957/150‑1 (MQ=37)
                                                                                                   cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCATTCAATTTTCTGCCCCACCAAAAGCAAAAACCCCCCCAAGCGGGTTTTTACGGAAACCCGGTGAAACAGACCGAAAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAc                                                >  1:2465251/1‑150 (MQ=11)
                                                                                                    gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTTACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGGACa                                               >  1:296037/1‑150 (MQ=255)
                                                                                                     aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCTAGGACAg                                              >  2:3327955/1‑150 (MQ=11)
                                                                                                      tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                                             <  1:664045/150‑1 (MQ=37)
                                                                                                        atccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTATTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                           <  1:4648568/148‑1 (MQ=37)
                                                                                                        agacGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCCTTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAACCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                           <  2:1242972/147‑1 (MQ=255)
                                                                                                         gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                          <  1:772517/150‑1 (MQ=37)
                                                                                                         gCCGCGATCTTTTGGGTATAGTGTCGTGGACAGACTTTCTTTTTTCTGCCACTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                          <  2:1654740/150‑1 (MQ=11)
                                                                                                           cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCCCCACCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGAACGATAAGCCGCTTTCTTTAGGGGACAGAGTCATGGAAAGtaaatc                                        >  2:4079453/1‑145 (MQ=11)
                                                                                                                   tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCt                                <  1:391863/150‑1 (MQ=37)
                                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTg                               <  1:1018954/150‑1 (MQ=37)
                                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCGTTGGGGTATAGCGTCGTGGACAGTCATTGCTCGTTCTg                               >  1:3977915/1‑150 (MQ=11)
                                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                             >  2:5023720/1‑72 (MQ=34)
                                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                             <  1:5023731/72‑1 (MQ=34)
                                                                                                                    ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                             >  2:5023731/1‑72 (MQ=34)
                                                                                                                        gTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGTGGCCAGGCATACATAGTTCTGcccc                           >  2:377580/1‑150 (MQ=11)
                                                                                                                               gtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTGTTACGTAAACCAGGTGAAACTGTCCGATAAGCCGCGGTCTGTTGGGGCTAGCGTCGGGGCCAGGCCTGCATGTTACTGccccggcagtc                    >  2:160547/1‑143 (MQ=11)
                                                                                                                                  gTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTAGAGCGTCGTGGGCAGTCATTGATGGTGCTGCGCCGGCAAAAGCa                 >  1:4968563/1‑150 (MQ=11)
                                                                                                                                  gTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAAATGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCATTCAACTTTCTGCACCTCCAAAAGCa                 >  2:4612499/1‑150 (MQ=11)
                                                                                                                                   tGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCaa                >  1:3348087/1‑150 (MQ=255)
                                                                                                                                   tGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGTGGCCAGGCAAGCAAGTTGAGGCCCCACCAAAAGCaa                >  2:4026234/1‑150 (MQ=11)
                                                                                                                                     gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACCGTCCATCATGGTTATGCCCCGCCAAAAGccata              >  2:914383/1‑146 (MQ=11)
                                                                                                                                        aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCATTCGGCCCCAGCCACAGCAAAAAcc           >  2:3157892/1‑150 (MQ=11)
                                                                                                                                         gTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGCCAGGCATTCAACAGTCTGCCCCGCCCACAGCCCAAAccc          >  1:2701975/1‑150 (MQ=11)
                                                                                                                                          tcattcatTTTTCTGCCCCTCCAAAAGCAAAATCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACccg         >  2:1954290/1‑150 (MQ=255)
                                                                                                                                           cattcattTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc        >  2:2026155/1‑150 (MQ=37)
                                                                                                                                            attcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTGGGGGATAGAGTCTATGACAGTCCATCAAAATTAAGCCCCAACAAAATCAAAAACccgcc       >  1:2585416/1‑150 (MQ=11)
                                                                                                                                             ttcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTATTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTAGGGTAGAGCGACGAGGAAAGTCATTCATATTACAGCACCTGCAAGAACAAAAACccgccg      >  2:3457935/1‑150 (MQ=11)
                                                                                                                                               catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTTCCCCTACAAAAGCAAAAACCCGCCGaa    >  2:2660038/1‑150 (MQ=255)
                                                                                                                                               catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCACCTCCAAAAGCAAAAACccgccgca    >  2:2778763/1‑148 (MQ=255)
                                                                                                                                               catTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGTGGACAGGCAGTCATCGTTCTGCCCCGCCAAAAGCAAAAACccgccgca    >  2:4223173/1‑148 (MQ=11)
                                                                                                                                                 tCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGc  <  2:608781/150‑1 (MQ=37)
                                                                                                                                                  |                                                                                                                                                    
TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGC  >  CP009273/3263264‑3263557

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: