Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,463,946 C→A 8.7% V2V (GTC→GTA bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,463,9460CA 8.7% 1274.3 / 69.0 358V2V (GTC→GTAbdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base C (161/166);  new base A (17/14);  total (178/180)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 364 total aligned reads displayed.

CGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGA  >  CP009273/4463797‑4464094
                                                                                                                                                     |                                                                                                                                                    
cgacgTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTc                                                                                                                                                      <  1:2297055/150‑1 (MQ=255)
cgacgTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTc                                                                                                                                                      >  1:3842330/1‑150 (MQ=255)
 gacgTACGATAGCGGCCCCGCTACCACCTCTGCCACCGAGGATGAGAACTGTCTGACCTGTAAAAGCGCCCAGAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCa                                                                                                                                                     <  2:2765240/150‑1 (MQ=255)
 gacgTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCa                                                                                                                                                     <  2:521818/150‑1 (MQ=255)
   cgTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACt                                                                                                                                                   <  1:5172770/150‑1 (MQ=255)
       cGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTTATGGTCACTaaaa                                                                                                                                               >  2:4050836/1‑150 (MQ=255)
          tAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAc                                                                                                                                            >  2:648467/1‑150 (MQ=255)
            gCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACaa                                                                                                                                          <  1:3761050/150‑1 (MQ=255)
             cGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAt                                                                                                                                         <  1:1096464/150‑1 (MQ=255)
             cGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAt                                                                                                                                         <  1:1673089/150‑1 (MQ=255)
              ggCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtc                                                                                                                                        >  1:1623230/1‑150 (MQ=255)
               gCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAAtct                                                                                                                                       <  1:3990732/150‑1 (MQ=255)
                 aCCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTcg                                                                                                                                     <  2:4130307/150‑1 (MQ=255)
                     aTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGtt                                                                                                                                 <  1:247658/150‑1 (MQ=255)
                       aCCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTcc                                                                                                                               <  2:59120/150‑1 (MQ=255)
                            gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGgtc                                                                                                                          <  2:5061216/150‑1 (MQ=255)
                             aCTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGgtcg                                                                                                                         <  1:27309/150‑1 (MQ=255)
                                gCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATGTATTATGTAACATGCATTACAAAACTGTTGTAACTTTCTGTCAACAGGTGTTGTAATGGTCACGAAAAAACAATCGCGCGTTCCAGGTCGTcc                                                                                                                      >  2:3957733/1‑150 (MQ=255)
                                  cACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCa                                                                                                                    <  2:498433/150‑1 (MQ=255)
                                  cACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCa                                                                                                                    <  1:2077246/150‑1 (MQ=255)
                                      gAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACg                                                                                                                <  2:447786/150‑1 (MQ=255)
                                       aGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGt                                                                                                               >  1:3309040/1‑150 (MQ=255)
                                          aTGATAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTc                                                                                                            <  1:928880/150‑1 (MQ=255)
                                           tGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                                           <  2:4776594/150‑1 (MQ=255)
                                           tGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                                           >  1:4873606/1‑150 (MQ=255)
                                           tGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                                           <  2:2177344/150‑1 (MQ=255)
                                            gagaACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGc                                                                                                          >  2:3464149/1‑150 (MQ=255)
                                                aCTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                                      <  2:3352618/150‑1 (MQ=255)
                                                      ttACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAg                                                                                                <  2:3248193/150‑1 (MQ=255)
                                                       tACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAgg                                                                                               <  2:2226305/150‑1 (MQ=255)
                                                          cTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCaa                                                                                            >  1:217062/1‑150 (MQ=255)
                                                          cTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCaa                                                                                            >  1:3883911/1‑150 (MQ=255)
                                                          cTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGATTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCaa                                                                                            >  1:2027877/1‑150 (MQ=255)
                                                              aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                                        >  2:235875/1‑150 (MQ=255)
                                                                 aGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGc                                                                                     <  2:470846/150‑1 (MQ=255)
                                                                  gcgcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCg                                                                                    <  1:561696/150‑1 (MQ=255)
                                                                   cgcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCgg                                                                                   <  2:2484387/150‑1 (MQ=255)
                                                                    gcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                                  <  2:301237/150‑1 (MQ=255)
                                                                    gcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                                  >  1:1982579/1‑150 (MQ=255)
                                                                      ccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                                <  2:4841477/150‑1 (MQ=255)
                                                                       cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaaa                                                                               >  2:492251/1‑150 (MQ=255)
                                                                           aGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTg                                                                           >  2:681371/1‑150 (MQ=255)
                                                                           aGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTg                                                                           <  1:711901/150‑1 (MQ=255)
                                                                           aGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTg                                                                           <  1:775110/150‑1 (MQ=255)
                                                                            ggACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGc                                                                          <  2:2867045/150‑1 (MQ=255)
                                                                            ggACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGc                                                                          >  1:623203/1‑150 (MQ=255)
                                                                             gACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCACAAGTGCt                                                                         >  1:4997117/1‑150 (MQ=255)
                                                                                 ccTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                                     >  1:3599283/1‑150 (MQ=255)
                                                                                   ttGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCa                                                                   <  2:3873242/150‑1 (MQ=255)
                                                                                   ttGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCa                                                                   <  2:1839129/150‑1 (MQ=255)
                                                                                    tGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                                  <  1:1307157/150‑1 (MQ=255)
                                                                                        ttattaTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                              >  2:15574/1‑150 (MQ=255)
                                                                                          attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaaa                                                            >  1:5013536/1‑150 (MQ=255)
                                                                                          attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaaa                                                            >  2:148032/1‑150 (MQ=255)
                                                                                            tatTTAACATGCATTACAAAACTGTTTTAACTTTCTTTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                                          <  1:3087950/150‑1 (MQ=255)
                                                                                            tatGTAACATGCATTACAATTATGTTTTAATTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCTTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                                          <  2:3065154/150‑1 (MQ=255)
                                                                                            tatGTAACATGCATTACAAAACTGTTTTAATTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                                          <  2:3756680/150‑1 (MQ=255)
                                                                                              tGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                                        >  2:712306/1‑150 (MQ=255)
                                                                                               gTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGttt                                                       >  1:2924072/1‑150 (MQ=255)
                                                                                                   cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                                   <  2:3718797/150‑1 (MQ=255)
                                                                                                   cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                                   <  2:3372681/150‑1 (MQ=255)
                                                                                                          tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCa                                            <  1:1765699/150‑1 (MQ=255)
                                                                                                          tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAACAAAAAGGTTTCGATGCTGTCa                                            >  2:3626510/1‑150 (MQ=255)
                                                                                                          tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCAAGACGTTTCGATCCAGAGCAGGCAATCTCTGAGGCAAAAGTGCTTTATCACCAAAAAGGTTTCGATGCTGTCa                                            >  2:1251622/1‑150 (MQ=255)
                                                                                                           aCAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                                           <  1:3191925/150‑1 (MQ=255)
                                                                                                              aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgtgt                                        >  2:1447659/1‑150 (MQ=255)
                                                                                                               aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGtt                                       >  2:3683949/1‑150 (MQ=255)
                                                                                                                 cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAGCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                                     >  1:4579289/1‑150 (MQ=255)
                                                                                                                 cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                                     >  2:1930850/1‑150 (MQ=255)
                                                                                                                  tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                                    >  1:1212797/1‑150 (MQ=255)
                                                                                                                  tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                                    <  2:1623230/150‑1 (MQ=255)
                                                                                                                  tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                                    <  2:2589660/150‑1 (MQ=255)
                                                                                                                    ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGa                                  <  1:4395871/150‑1 (MQ=255)
                                                                                                                    ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGa                                  >  2:1065944/1‑150 (MQ=255)
                                                                                                                     tttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCACAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGaa                                 >  2:3809354/1‑150 (MQ=255)
                                                                                                                        aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt                              <  1:1858345/150‑1 (MQ=255)
                                                                                                                        aaCTTTCTGTCAACAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt                              <  1:1500035/150‑1 (MQ=255)
                                                                                                                         aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                             <  2:841388/150‑1 (MQ=255)
                                                                                                                         aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                             >  1:105333/1‑150 (MQ=255)
                                                                                                                            ttCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTg                          <  2:2150318/150‑1 (MQ=255)
                                                                                                                             tCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGGCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa                         <  2:4176773/150‑1 (MQ=255)
                                                                                                                               tGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAtt                       >  2:1388902/1‑150 (MQ=255)
                                                                                                                                gTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCCATCTCGGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAGTa                      >  2:2707963/1‑150 (MQ=255)
                                                                                                                                 tCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTAt                     >  1:3442321/1‑150 (MQ=255)
                                                                                                                                  cAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCACTCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATc                    >  2:1122713/1‑150 (MQ=255)
                                                                                                                                   aaCAGGTTTTGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt                   >  1:3857077/1‑150 (MQ=255)
                                                                                                                                      aGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATACTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg                >  1:1735612/1‑150 (MQ=255)
                                                                                                                                         ttttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTAt             >  2:5090268/1‑150 (MQ=255)
                                                                                                                                          tttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTAtt            >  1:904915/1‑150 (MQ=255)
                                                                                                                                          tttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTAtt            >  1:3529537/1‑150 (MQ=255)
                                                                                                                                          tttGTAATGGTAACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTAtt            <  1:2643916/150‑1 (MQ=255)
                                                                                                                                           ttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTa           >  2:2037148/1‑150 (MQ=255)
                                                                                                                                           ttGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTa           <  1:685550/150‑1 (MQ=255)
                                                                                                                                             gTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAc         >  1:4028883/1‑150 (MQ=255)
                                                                                                                                             gTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAc         >  1:3453050/1‑150 (MQ=255)
                                                                                                                                              tAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAcc        <  1:1390668/150‑1 (MQ=255)
                                                                                                                                               aaTGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAccc       >  1:1889437/1‑150 (MQ=255)
                                                                                                                                                 tGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAccccc     >  2:3639053/1‑150 (MQ=255)
                                                                                                                                                 tGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAAccccc     <  1:232660/150‑1 (MQ=255)
                                                                                                                                                    tCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGa  >  1:5077789/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTATTAACCCCCCGA  >  CP009273/4463797‑4464094

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