Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,088,734 G→T 5.3% V135F (GTT→TTT)  ycdT → diguanylate cyclase, membrane‑anchored

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,088,7340GT 5.3% 924.3 / 17.1 242V135F (GTT→TTT) ycdTdiguanylate cyclase, membrane‑anchored
Reads supporting (aligned to +/- strand):  ref base G (97/132);  new base T (7/6);  total (104/138)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 243 total aligned reads displayed.

GATTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACC  >  CP009273/1088586‑1088880
                                                                                                                                                    |                                                                                                                                                  
gATTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                                                                                                                    <  2:2051097/149‑1 (MQ=255)
 aTTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                                                                                                                   <  2:4670832/149‑1 (MQ=255)
 aTTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                                                                                                                   <  1:680545/149‑1 (MQ=255)
 aTTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                                                                                                                   <  2:1589110/149‑1 (MQ=255)
  tttACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                                                                                                                  >  2:4432793/1‑149 (MQ=255)
   ttACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                                                                                                                 <  1:614409/149‑1 (MQ=255)
    tACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                                                                                                                <  2:2194940/149‑1 (MQ=255)
           ttATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                                                                                                                         >  2:365092/1‑149 (MQ=255)
              taACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                                                                                                                      <  1:48000/149‑1 (MQ=255)
              taACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                                                                                                                      >  1:4056411/1‑149 (MQ=255)
              taACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                                                                                                                      <  1:314109/149‑1 (MQ=255)
               aaCCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                                                                                                                     <  1:4688129/149‑1 (MQ=255)
                 ccGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                                                                                                                   <  1:1472328/149‑1 (MQ=255)
                       ttATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcgc                                                                                                                             <  1:869137/149‑1 (MQ=255)
                          tCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                                                                                                                          <  2:2406423/149‑1 (MQ=255)
                           cATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                                                                                                                         >  2:708676/1‑149 (MQ=255)
                           cATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                                                                                                                         <  1:2468846/149‑1 (MQ=255)
                             ttCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAattta                                                                                                                       >  1:1783982/1‑149 (MQ=255)
                             ttCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                                                                                                                       >  2:3747620/1‑149 (MQ=255)
                               cAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                                                                                                                     <  1:2093259/149‑1 (MQ=255)
                                aGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttattt                                                                                                                    >  1:4194977/1‑149 (MQ=255)
                                aGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAAAGCGCAatttattt                                                                                                                    >  2:4575198/1‑149 (MQ=255)
                                 gCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                                                                                                                   <  1:3010670/149‑1 (MQ=255)
                                          tGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa                                                                                                          <  2:4246738/149‑1 (MQ=255)
                                           gaggagCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt                                                                                                         <  2:252310/149‑1 (MQ=255)
                                           gaggagCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt                                                                                                         <  2:1867962/149‑1 (MQ=255)
                                              gagCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc                                                                                                      <  1:235937/149‑1 (MQ=255)
                                              gagCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc                                                                                                      <  2:4153060/149‑1 (MQ=255)
                                               agCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt                                                                                                     <  2:1779940/149‑1 (MQ=255)
                                               agCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt                                                                                                     <  1:665082/149‑1 (MQ=255)
                                               agCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt                                                                                                     <  2:4168958/149‑1 (MQ=255)
                                                gCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt                                                                                                    >  2:3871016/1‑149 (MQ=255)
                                                  gTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg                                                                                                  <  1:2507185/149‑1 (MQ=255)
                                                   tttGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCACTTTATTTTCCAAAAAAATGACTTtgt                                                                                                 >  1:522949/1‑149 (MQ=255)
                                                     tGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGc                                                                                               >  1:3979182/1‑149 (MQ=255)
                                                      gACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCa                                                                                              >  2:2620733/1‑149 (MQ=255)
                                                        cAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATa                                                                                            >  2:3797954/1‑149 (MQ=255)
                                                         aaGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATaa                                                                                           <  2:4497945/149‑1 (MQ=255)
                                                         aaGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATaa                                                                                           >  2:172771/1‑149 (MQ=255)
                                                          aGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAg                                                                                          <  2:4566984/149‑1 (MQ=255)
                                                          aGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAg                                                                                          <  2:1406897/149‑1 (MQ=255)
                                                            cAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtt                                                                                        <  1:1373866/149‑1 (MQ=255)
                                                            cAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtt                                                                                        <  2:2698533/149‑1 (MQ=255)
                                                            cAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtt                                                                                        <  1:2590974/149‑1 (MQ=255)
                                                                 ttgtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAgtttgttt                                                                                   <  1:2782705/149‑1 (MQ=255)
                                                                     ccAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttg                                                                               <  2:785051/149‑1 (MQ=255)
                                                                     ccAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttg                                                                               <  2:678516/149‑1 (MQ=255)
                                                                      cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgt                                                                              >  1:4621493/1‑149 (MQ=255)
                                                                      cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgt                                                                              <  1:591219/149‑1 (MQ=255)
                                                                       aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgtt                                                                             >  2:532593/1‑149 (MQ=255)
                                                                       aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgtt                                                                             <  2:1992227/149‑1 (MQ=255)
                                                                        aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgttt                                                                            <  2:4678242/149‑1 (MQ=255)
                                                                        aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgttt                                                                            <  2:2953963/149‑1 (MQ=255)
                                                                        aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGtttgttttttgttt                                                                            >  1:3308051/1‑149 (MQ=255)
                                                                          aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTtgttttttgttttt                                                                          >  2:3127377/1‑149 (MQ=255)
                                                                          aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTtgttttttgttttt                                                                          >  1:197982/1‑149 (MQ=255)
                                                                             aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTggg                                                                       <  1:2052016/149‑1 (MQ=255)
                                                                                     aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGTATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCAAAAAAAATGACTTTGTGAATAAGTTTGTTTTTTGTTTTTgggggggggtt                                                               >  2:4387785/1‑144 (MQ=255)
                                                                                      tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATa                                                              >  2:904932/1‑149 (MQ=255)
                                                                                       cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAg                                                             <  2:3157115/149‑1 (MQ=255)
                                                                                       cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAg                                                             >  1:1952869/1‑149 (MQ=255)
                                                                                        tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGt                                                            >  1:1310606/1‑149 (MQ=255)
                                                                                         attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGtt                                                           >  2:1561745/1‑149 (MQ=255)
                                                                                         attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGtt                                                           <  1:2965880/149‑1 (MQ=255)
                                                                                             tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTc                                                       <  1:928436/149‑1 (MQ=255)
                                                                                             tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTc                                                       <  1:1579310/149‑1 (MQ=255)
                                                                                               tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcac                                                     <  2:1244418/149‑1 (MQ=255)
                                                                                                gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcaca                                                    <  2:613979/149‑1 (MQ=255)
                                                                                                gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcaca                                                    >  2:489995/1‑149 (MQ=255)
                                                                                                gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcaca                                                    >  1:3523859/1‑149 (MQ=255)
                                                                                                 tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACAt                                                   <  2:2031352/149‑1 (MQ=255)
                                                                                                 tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACAt                                                   <  2:233682/149‑1 (MQ=255)
                                                                                                 tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACAt                                                   <  2:3304693/149‑1 (MQ=255)
                                                                                                   tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATAc                                                 >  1:1737202/1‑149 (MQ=255)
                                                                                                     gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACtt                                               >  1:3572827/1‑149 (MQ=255)
                                                                                                       cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTc                                             >  1:3854008/1‑149 (MQ=255)
                                                                                                        aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCt                                            >  2:784458/1‑149 (MQ=255)
                                                                                                         gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTa                                           >  2:3053478/1‑149 (MQ=255)
                                                                                                         gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTa                                           <  2:169761/149‑1 (MQ=255)
                                                                                                             aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCa                                       <  2:1281791/149‑1 (MQ=255)
                                                                                                             aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCa                                       >  1:3705212/1‑149 (MQ=255)
                                                                                                             aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCa                                       >  1:2667238/1‑149 (MQ=255)
                                                                                                               tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAtt                                     <  1:4432793/149‑1 (MQ=255)
                                                                                                                ttGTGCATATTAATATTTCTGGCATTGGTGAATAAATTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTa                                    <  1:2835956/149‑1 (MQ=255)
                                                                                                                    gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgag                                <  1:1068667/149‑1 (MQ=255)
                                                                                                                    gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgag                                <  1:2760906/149‑1 (MQ=255)
                                                                                                                     cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgaga                               <  2:2760735/149‑1 (MQ=255)
                                                                                                                       tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGc                             <  1:1111946/149‑1 (MQ=255)
                                                                                                                       tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGc                             >  1:2421945/1‑149 (MQ=255)
                                                                                                                            atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACaa                        >  2:1530962/1‑149 (MQ=255)
                                                                                                                            atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACaa                        >  2:1051721/1‑149 (MQ=255)
                                                                                                                             tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAAt                       <  2:1646260/149‑1 (MQ=255)
                                                                                                                               tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAAttt                     >  1:4736050/1‑149 (MQ=255)
                                                                                                                                  cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTaca                  <  1:3809080/149‑1 (MQ=255)
                                                                                                                                     gCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAt               <  1:825528/149‑1 (MQ=255)
                                                                                                                                        ttGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGc            <  2:347637/149‑1 (MQ=255)
                                                                                                                                        ttGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGc            >  1:578695/1‑149 (MQ=255)
                                                                                                                                           gTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGa         >  2:871508/1‑149 (MQ=255)
                                                                                                                                             gAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGt       <  1:1095922/149‑1 (MQ=255)
                                                                                                                                              aaTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGtt      >  2:3006245/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATTTACTTTGTTATAACCGTCATTATCATTCAGCAAATTATTGAGGAGCGTTTGACAAGCAGTGTTGTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACC  >  CP009273/1088586‑1088880

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: ‑