Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,212,264 G→T 34.2% H58N (CAC→AAC)  aer ← fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,212,2640GT34.2% 556.4 / 291.1 341H58N (CAC→AAC) aerfused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base G (108/116);  new base T (62/55);  total (170/171)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

TTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGG  >  CP009273/3212117‑3212413
                                                                                                                                                   |                                                                                                                                                     
ttCCCTCGCGCACCATCGGTACCGCATGGGCCCGCACCCAATAATGGTAACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTTAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGAc                                                                                                                                                     <  1:1830324/150‑1 (MQ=255)
 tCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACg                                                                                                                                                    <  2:778909/150‑1 (MQ=255)
 tCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACg                                                                                                                                                    >  1:3386262/1‑150 (MQ=255)
    cTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCAc                                                                                                                                                 >  1:4940223/1‑150 (MQ=255)
     tCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCAcc                                                                                                                                                >  2:4612141/1‑150 (MQ=255)
          cACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCAtgtt                                                                                                                                           <  1:4666667/150‑1 (MQ=255)
             cATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTtgtg                                                                                                                                        >  1:4244561/1‑150 (MQ=255)
                 ggTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGc                                                                                                                                    >  1:1120047/1‑150 (MQ=255)
                 ggTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGc                                                                                                                                    >  2:5194666/1‑150 (MQ=255)
                  gTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACTATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCt                                                                                                                                   >  2:5252662/1‑150 (MQ=255)
                  gTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCt                                                                                                                                   <  2:3228357/150‑1 (MQ=255)
                  gTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCt                                                                                                                                   <  2:4774426/150‑1 (MQ=255)
                      cGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGccc                                                                                                                               >  1:2391336/1‑150 (MQ=255)
                      cGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGccc                                                                                                                               <  1:1989865/150‑1 (MQ=255)
                       gCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCACAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCt                                                                                                                              >  1:4783071/1‑150 (MQ=255)
                        cATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCtt                                                                                                                             <  1:2572905/150‑1 (MQ=255)
                         aTTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTg                                                                                                                            <  1:1734355/150‑1 (MQ=255)
                         aTTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCGTg                                                                                                                            >  2:1781319/1‑150 (MQ=255)
                          ttGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGt                                                                                                                           >  1:4170691/1‑150 (MQ=255)
                          ttGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGt                                                                                                                           >  2:4350722/1‑150 (MQ=255)
                          ttGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGt                                                                                                                           <  2:1716032/150‑1 (MQ=255)
                            ggCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGATGCACCATGTTGTGCGGCTGCCCTTGTaa                                                                                                                         <  1:2120637/150‑1 (MQ=255)
                                 gCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCtt                                                                                                                    <  1:3499347/150‑1 (MQ=255)
                                     ccAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCaa                                                                                                                >  2:2318536/1‑150 (MQ=255)
                                     ccAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCaa                                                                                                                >  1:4088694/1‑150 (MQ=255)
                                       aataatGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAgg                                                                                                              >  1:3080467/1‑150 (MQ=255)
                                         taatGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGta                                                                                                            <  1:1775104/150‑1 (MQ=255)
                                          aatGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGtat                                                                                                           >  2:2832235/1‑150 (MQ=255)
                                              gTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGcc                                                                                                       <  2:3566200/150‑1 (MQ=255)
                                                 aCCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCt                                                                                                    <  2:593271/150‑1 (MQ=255)
                                                            ggcgATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACaa                                                                                         <  1:2736329/150‑1 (MQ=255)
                                                                aTTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAgt                                                                                     >  2:2560470/1‑150 (MQ=255)
                                                                 tttttCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCACACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAgtg                                                                                    <  1:2299804/150‑1 (MQ=255)
                                                                 tttttCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAgtg                                                                                    <  2:3770114/150‑1 (MQ=255)
                                                                  ttttCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAgtgt                                                                                   >  2:1928923/1‑150 (MQ=255)
                                                                    ttCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCa                                                                                 >  2:3947889/1‑150 (MQ=255)
                                                                      cACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCAtt                                                                               <  2:4506698/150‑1 (MQ=255)
                                                                       aCGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTa                                                                              >  1:1587832/1‑150 (MQ=255)
                                                                       aCGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTa                                                                              <  2:2993091/150‑1 (MQ=255)
                                                                        cGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAg                                                                             <  1:104998/150‑1 (MQ=255)
                                                                        cGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAg                                                                             <  2:5042365/150‑1 (MQ=255)
                                                                         gATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCAGTAGc                                                                            >  1:648657/1‑150 (MQ=255)
                                                                         gATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGGATAGCCGCTCACCTGCACAAAAGTGTCATTAGc                                                                            >  1:4703223/1‑150 (MQ=255)
                                                                                  ccAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGat                                                                   <  1:431743/150‑1 (MQ=255)
                                                                                  ccAGGGCTCCCCTTTTGTCAGGGTGAACCACATCTCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGat                                                                   <  2:964478/150‑1 (MQ=255)
                                                                                    aGGGCTCCCCTTTTTTCAGGGTTAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCCCCATGGTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGatat                                                                 >  1:1409155/1‑150 (MQ=255)
                                                                                    aGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGatat                                                                 >  2:408286/1‑150 (MQ=255)
                                                                                      ggCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAg                                                               >  2:1364773/1‑150 (MQ=255)
                                                                                      ggCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAg                                                               >  1:5296368/1‑150 (MQ=255)
                                                                                      ggCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAg                                                               <  1:2086884/150‑1 (MQ=255)
                                                                                        cTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCt                                                             >  1:3433621/1‑150 (MQ=255)
                                                                                        cTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCt                                                             <  1:5211829/150‑1 (MQ=255)
                                                                                         tCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCtt                                                            >  2:3537013/1‑150 (MQ=255)
                                                                                          ccccTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCttt                                                           <  2:1108708/150‑1 (MQ=255)
                                                                                          ccccTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCttt                                                           <  2:4587050/150‑1 (MQ=255)
                                                                                           cccTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTg                                                          <  1:2878351/150‑1 (MQ=255)
                                                                                            ccTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGc                                                         <  2:2947606/150‑1 (MQ=255)
                                                                                               tttttCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGa                                                      >  2:4274473/1‑150 (MQ=255)
                                                                                                ttttCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGAt                                                     <  1:962171/150‑1 (MQ=255)
                                                                                                  ttCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGCGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCg                                                   >  1:827483/1‑150 (MQ=255)
                                                                                                   tCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCgg                                                  >  1:4511742/1‑150 (MQ=255)
                                                                                                    cAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGt                                                 >  1:4801294/1‑150 (MQ=255)
                                                                                                    cAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGt                                                 >  2:3711837/1‑150 (MQ=255)
                                                                                                     aGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTa                                                >  2:3033058/1‑150 (MQ=255)
                                                                                                     aGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTa                                                <  1:1502803/150‑1 (MQ=255)
                                                                                                       ggTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGt                                              <  1:4011820/150‑1 (MQ=255)
                                                                                                       ggTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGt                                              <  1:2506852/150‑1 (MQ=255)
                                                                                                         tGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTgg                                            >  2:2630567/1‑150 (MQ=255)
                                                                                                          gAACCACATCTCCGCACACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGa                                           <  2:87979/150‑1 (MQ=255)
                                                                                                           aaCCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGAc                                          >  1:1287477/1‑150 (MQ=255)
                                                                                                            aCCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACa                                         >  2:318421/1‑150 (MQ=255)
                                                                                                              cacaTATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATc                                       <  2:2748370/150‑1 (MQ=255)
                                                                                                               acaTATCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCa                                      >  2:3989640/1‑150 (MQ=255)
                                                                                                                caTATCCGCAAACGCCTCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAAATCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATTAGTGATATAGCTTTGCAGATCGGTAGTGGACATCag                                     <  1:5181348/150‑1 (MQ=255)
                                                                                                                 atatCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCaga                                    >  2:4964245/1‑150 (MQ=255)
                                                                                                                    tCCGCAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTg                                 <  1:1004359/150‑1 (MQ=255)
                                                                                                                       gCAAACGCCGC‑TTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTAt                             <  2:269826/150‑1 (MQ=255)
                                                                                                                        cAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTAt                             >  1:739498/1‑150 (MQ=255)
                                                                                                                        cAAACGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTAt                             <  2:3837858/150‑1 (MQ=255)
                                                                                                                         aaaCGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATc                            >  1:3179070/1‑150 (MQ=255)
                                                                                                                          aaCGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCg                           >  1:3266228/1‑150 (MQ=255)
                                                                                                                           aCGCCGCTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGt                          <  1:4612141/150‑1 (MQ=255)
                                                                                                                              ccgcTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCg                       >  2:936611/1‑150 (MQ=255)
                                                                                                                               cgcTTTTGGCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGc                      <  2:1252043/150‑1 (MQ=255)
                                                                                                                               cgcTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGc                      >  2:3352481/1‑150 (MQ=255)
                                                                                                                               cgcTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGc                      >  1:2005504/1‑150 (MQ=255)
                                                                                                                                gcTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGcc                     >  2:4183642/1‑150 (MQ=255)
                                                                                                                                    ttGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCg                 <  1:3947728/150‑1 (MQ=255)
                                                                                                                                       gCATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCgggg              >  1:1627804/1‑150 (MQ=255)
                                                                                                                                        cATATCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGt             <  1:1467029/150‑1 (MQ=255)
                                                                                                                                         atatCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTa            <  2:3178456/150‑1 (MQ=255)
                                                                                                                                         atatCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTa            >  2:2500918/1‑150 (MQ=255)
                                                                                                                                         atatCCGGGTTACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTa            >  2:630510/1‑150 (MQ=255)
                                                                                                                                          tatCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTAt           <  2:1993274/150‑1 (MQ=255)
                                                                                                                                          tatCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTAt           <  2:1560340/150‑1 (MQ=255)
                                                                                                                                          tatCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCAGTAGCATGAGTGCTATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGtct           >  2:4043856/1‑148 (MQ=255)
                                                                                                                                           atCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTAtt          >  2:4637320/1‑150 (MQ=255)
                                                                                                                                           atCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTAtt          >  2:2168402/1‑150 (MQ=255)
                                                                                                                                             ccGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTAtttt        <  1:2994030/150‑1 (MQ=255)
                                                                                                                                                   gACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTggg  <  2:4603614/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTGACGCACCATGTTGTGCGGCTGCCCTTGTAACTCTTGCAAGGTATAGCCGCTCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGG  >  CP009273/3212117‑3212413

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