Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,232,020 C→A 5.7% E226D (GAG→GAT ycgB ← SpoVR family stationary phase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,232,0200CA 5.7% 784.9 / 15.5 209E226D (GAG→GATycgBSpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base C (105/92);  new base A (8/4);  total (113/96)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 215 total aligned reads displayed.

TTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTTCATCAGCGCGT  >  CP009273/1231873‑1232168
                                                                                                                                                   |                                                                                                                                                    
 tGCGGATAAAAATACTGGCTCACCTTACGCACAATACTCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGTTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTtctc                                                                                                                                                    <  1:4013770/149‑1 (MQ=255)
     gATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTcc                                                                                                                                                <  1:1731304/149‑1 (MQ=255)
     gATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCGTCAGCAACCGTTTTCTCTTcc                                                                                                                                                <  1:2557770/149‑1 (MQ=255)
      aTAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGATCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCt                                                                                                                                               >  1:2090258/1‑149 (MQ=255)
       tAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTc                                                                                                                                              >  2:2739240/1‑149 (MQ=255)
            aTACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCt                                                                                                                                         >  2:79480/1‑149 (MQ=255)
                 ggCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGc                                                                                                                                    <  2:4056139/149‑1 (MQ=255)
                    tCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCTGTGGCGCATTTTTCTCCATAAAAGAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAgg                                                                                                                                 <  1:1543964/149‑1 (MQ=255)
                    tCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAgg                                                                                                                                 <  2:1763689/149‑1 (MQ=255)
                      aCCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGt                                                                                                                               >  2:2973659/1‑149 (MQ=255)
                      aCCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGt                                                                                                                               <  1:2692505/149‑1 (MQ=255)
                         ttACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACg                                                                                                                            >  1:638572/1‑149 (MQ=255)
                           aCGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGcc                                                                                                                          >  1:278433/1‑149 (MQ=255)
                             gCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCag                                                                                                                        >  2:2637339/1‑149 (MQ=255)
                             gCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCag                                                                                                                        >  2:1028400/1‑149 (MQ=255)
                               acaATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCagag                                                                                                                      >  1:97183/1‑149 (MQ=255)
                                  aTACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCat                                                                                                                   <  1:2167512/149‑1 (MQ=255)
                                    aCGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCatat                                                                                                                 <  2:2287832/149‑1 (MQ=255)
                                        aGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGAc                                                                                                             >  1:3230479/1‑149 (MQ=255)
                                        aGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGAc                                                                                                             >  1:3641596/1‑149 (MQ=255)
                                          gATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCt                                                                                                           <  1:1298304/149‑1 (MQ=255)
                                            tttCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGa                                                                                                         >  1:562244/1‑149 (MQ=255)
                                            tttCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGa                                                                                                         >  2:684365/1‑149 (MQ=255)
                                               cACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACtt                                                                                                      <  2:2615464/149‑1 (MQ=255)
                                               cACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACtt                                                                                                      >  2:3922841/1‑149 (MQ=255)
                                                aCGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACttt                                                                                                     <  2:2121316/149‑1 (MQ=255)
                                                aCGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACttt                                                                                                     <  2:43088/149‑1 (MQ=255)
                                                  gCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGc                                                                                                   <  2:3294482/149‑1 (MQ=255)
                                                  gCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGc                                                                                                   <  1:3783731/149‑1 (MQ=255)
                                                  gCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGc                                                                                                   <  1:2494389/149‑1 (MQ=255)
                                                  gCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGc                                                                                                   <  1:422630/149‑1 (MQ=255)
                                                   cTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCa                                                                                                  <  1:3014972/149‑1 (MQ=255)
                                                      ccATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGAt                                                                                               >  2:2104134/1‑149 (MQ=255)
                                                      ccATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGAt                                                                                               >  2:1394340/1‑149 (MQ=255)
                                                      ccATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGAt                                                                                               >  1:1030187/1‑149 (MQ=255)
                                                        aTGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATAc                                                                                             >  2:486964/1‑149 (MQ=255)
                                                          gATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATActc                                                                                           <  1:986569/149‑1 (MQ=255)
                                                            ttCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCtt                                                                                         <  2:3909109/149‑1 (MQ=255)
                                                               cagcagTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTcgc                                                                                      >  1:1229988/1‑149 (MQ=255)
                                                               cagcagTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTCTGGTTCGGACGGATAGCGCCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGCTACGCCAGAGCATATTGACCTGACGGTGCAGATACTCTTcgc                                                                                      >  2:3681215/1‑149 (MQ=255)
                                                                agcagTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTcgcg                                                                                     >  2:2945238/1‑149 (MQ=255)
                                                                 gcagTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTGCTCTCGCTTCGGCAGGGTACGCCATAGCATATTGACCTGACTTTGCAGATACTCTTCGCGa                                                                                    <  1:2593150/149‑1 (MQ=255)
                                                                       gcgcATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTc                                                                              >  2:3514686/1‑149 (MQ=255)
                                                                       gcgcATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTc                                                                              >  2:245484/1‑149 (MQ=255)
                                                                            tttttCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCg                                                                         >  1:2258003/1‑149 (MQ=255)
                                                                            tttttCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCg                                                                         >  2:2506162/1‑149 (MQ=255)
                                                                                tctcCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCt                                                                     >  2:3328137/1‑149 (MQ=255)
                                                                                  tcCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCttt                                                                   >  2:1996123/1‑149 (MQ=255)
                                                                                  tcCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCttt                                                                   >  2:1222127/1‑149 (MQ=255)
                                                                                  tcCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCttt                                                                   >  1:750021/1‑149 (MQ=255)
                                                                                   ccATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTtc                                                                  >  1:1101703/1‑149 (MQ=255)
                                                                                   ccATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTtc                                                                  <  2:3230479/149‑1 (MQ=255)
                                                                                    cATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTtct                                                                 >  1:3238248/1‑149 (MQ=255)
                                                                                    cATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTtct                                                                 >  1:3237942/1‑149 (MQ=255)
                                                                                      tAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTctct                                                               <  1:3461669/149‑1 (MQ=255)
                                                                                       aaaaTAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCtctt                                                              >  2:3623637/1‑149 (MQ=255)
                                                                                        aaaTAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCtcttc                                                             >  1:3666573/1‑149 (MQ=255)
                                                                                         aaTAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCtcttct                                                            <  1:3893311/149‑1 (MQ=255)
                                                                                          aTAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTcttctt                                                           <  2:422287/149‑1 (MQ=255)
                                                                                          aTAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTcttctt                                                           <  2:1605361/149‑1 (MQ=255)
                                                                                            agagCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGc                                                         <  1:61065/149‑1 (MQ=255)
                                                                                            agagCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGc                                                         <  1:1111509/149‑1 (MQ=255)
                                                                                            agagCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGc                                                         <  1:513042/149‑1 (MQ=255)
                                                                                             gagCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCACCCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCa                                                        <  2:1800924/149‑1 (MQ=255)
                                                                                             gagCAGGTTTCCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCa                                                        >  2:62158/1‑149 (MQ=255)
                                                                                               gCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGc                                                      >  1:435792/1‑149 (MQ=255)
                                                                                                 aGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCga                                                    <  1:1578934/149‑1 (MQ=255)
                                                                                                  ggTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCgag                                                   >  1:2098375/1‑149 (MQ=255)
                                                                                                   gTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCgaga                                                  >  2:736206/1‑149 (MQ=255)
                                                                                                     tttCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGAtt                                                <  1:201336/149‑1 (MQ=255)
                                                                                                     tttCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGAtt                                                >  2:1338460/1‑149 (MQ=255)
                                                                                                     tttCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGAtt                                                >  2:1338030/1‑149 (MQ=255)
                                                                                                         ttGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGAtttttt                                            <  2:1932896/149‑1 (MQ=255)
                                                                                                             ggTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCgg                                        <  1:1376706/149‑1 (MQ=255)
                                                                                                              gTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCggg                                       >  1:3248704/1‑149 (MQ=255)
                                                                                                                  ggacggaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGtt                                   <  2:3724776/149‑1 (MQ=255)
                                                                                                                  ggacggaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGtt                                   >  2:2549871/1‑149 (MQ=255)
                                                                                                                     cggaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGt                                <  2:562244/149‑1 (MQ=255)
                                                                                                                      ggaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTa                               >  1:1488007/1‑149 (MQ=255)
                                                                                                                      ggaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTa                               >  1:1786522/1‑149 (MQ=255)
                                                                                                                       gaTAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTAc                              <  2:2061632/149‑1 (MQ=255)
                                                                                                                           gcggcgCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGt                          <  2:471916/149‑1 (MQ=255)
                                                                                                                           gcggcgCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGt                          >  2:744568/1‑149 (MQ=255)
                                                                                                                           gcggcgCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGt                          <  1:3986738/149‑1 (MQ=255)
                                                                                                                              gcgcgcTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGtt                                   >  1:3909280/1‑137 (MQ=255)
                                                                                                                              gcgcgcTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCa                       >  1:2659951/1‑149 (MQ=255)
                                                                                                                              gcgcgcTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCa                       >  2:1231358/1‑149 (MQ=255)
                                                                                                                                  gcTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGcc                   >  2:2638906/1‑149 (MQ=255)
                                                                                                                                       aGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGt              <  1:1616820/149‑1 (MQ=255)
                                                                                                                                        gCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGtt             >  2:3144651/1‑149 (MQ=255)
                                                                                                                                        gCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGtt             >  2:4073957/1‑149 (MQ=255)
                                                                                                                                         cAACCGTTTTATCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtc            <  2:910780/149‑1 (MQ=255)
                                                                                                                                          aaCCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtca           >  1:962286/1‑149 (MQ=255)
                                                                                                                                          aaCCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtca           >  1:3268063/1‑149 (MQ=255)
                                                                                                                                          aaCCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtca           <  1:1937260/149‑1 (MQ=255)
                                                                                                                                          aaCCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCACTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtca           >  1:1810256/1‑149 (MQ=255)
                                                                                                                                             cGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtcatca        <  1:79480/149‑1 (MQ=255)
                                                                                                                                             cGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTtcatca        <  1:1863996/149‑1 (MQ=255)
                                                                                                                                               ttttCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTTCATCAgc      >  1:273891/1‑149 (MQ=255)
                                                                                                                                                  tctcTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTTCATCAGcgcg   >  1:2985784/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCGGATAAAAATACTGGCTCACCTTACGCACAATACGCAGGATTTCACGCTGCCATGATTCCAGCAGTGGCGCATTTTTCTCCATAAAATAGAGCAGGTTTTCTTGTGGTTCGGACGGATAGCGGCGCGCTTCAGCAACCGTTTTCTCTTCCTCGCGCTTCGGCAGGGTACGCCAGAGCATATTGACCTGACTTTGCAGATACTCTTCGCGACTTTTCTGCCGGGCTTTCTCTTCTTGCAGCGAGATTTTTTGCGGGCGTTTGTACCGGTCCACGCCGTAGTTCATCAGCGCGT  >  CP009273/1231873‑1232168

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: