Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 469,775 C→A 13.2% V281V (GTG→GTT tesB ← acyl‑CoA thioesterase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273469,7750CA13.2% 715.2 / 55.6 224V281V (GTG→GTTtesBacyl‑CoA thioesterase 2
Reads supporting (aligned to +/- strand):  ref base C (107/87);  new base A (10/20);  total (117/107)
Fisher's exact test for biased strand distribution p-value = 3.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

TTGGTCTGCGTGTACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGG  >  CP009273/469629‑469923
                                                                                                                                                  |                                                                                                                                                    
  ggTCTGCGTGTACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCAccc                                                                                                                                                  >  2:381772/1‑149 (MQ=255)
    tCTGCGTGTACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCt                                                                                                                                                >  2:1581451/1‑149 (MQ=255)
    tCTGCGTGTACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCt                                                                                                                                                <  1:708703/149‑1 (MQ=255)
           gtACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGaa                                                                                                                                         <  1:3871475/149‑1 (MQ=255)
           gtACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGaa                                                                                                                                         <  1:179006/149‑1 (MQ=255)
               cGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGt                                                                                                                                     <  1:157293/149‑1 (MQ=255)
                   agagCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAg                                                                                                                                 >  1:2811085/1‑149 (MQ=255)
                   agagCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAg                                                                                                                                 <  1:358857/149‑1 (MQ=255)
                     agCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGc                                                                                                                               <  2:713001/149‑1 (MQ=255)
                        aGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAAc                                                                                                                            <  1:3468872/149‑1 (MQ=255)
                                 aaGGGCTGGATGTCAACAGCCACGGCGAGAAGGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGcc                                                                                                                   <  2:1490482/149‑1 (MQ=255)
                                   gggCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGt                                                                                                                 >  1:1812885/1‑149 (MQ=255)
                                   gggCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGt                                                                                                                 >  1:3675833/1‑149 (MQ=255)
                                    ggCTGGATGTCAACAGCCACGGCGAGAATGCCTCTAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAGCGGTCGAGGCAACCAGTACGCCGTc                                                                                                                >  1:4085133/1‑149 (MQ=255)
                                    ggCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCACCCAGTACGCCggc                                                                                                                >  2:3039477/1‑147 (MQ=255)
                                    ggCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTc                                                                                                                >  1:2246904/1‑149 (MQ=255)
                                        ggATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTgg                                                                                                            >  1:3714133/1‑149 (MQ=255)
                                           tGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGta                                                                                                         <  2:3188294/149‑1 (MQ=255)
                                               aaCAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATaaa                                                                                                     >  2:68933/1‑149 (MQ=255)
                                                aCAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATaaaa                                                                                                    >  2:578281/1‑149 (MQ=255)
                                                aCAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCGTGGGTATaaaa                                                                                                    >  2:2743157/1‑149 (MQ=255)
                                                    ccACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCa                                                                                                >  1:2787725/1‑149 (MQ=255)
                                                     cACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCAc                                                                                               <  1:1245361/149‑1 (MQ=255)
                                                        ggCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACcgc                                                                                            >  2:3975695/1‑149 (MQ=255)
                                                        ggCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACcgc                                                                                            >  1:1505666/1‑149 (MQ=255)
                                                               aTGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACaaa                                                                                     <  1:569411/149‑1 (MQ=255)
                                                                 gCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGc                                                                                   <  2:1262533/149‑1 (MQ=255)
                                                                   cTATAACGCGAAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTAAATATTTAATTCTGATTACGCATAACCCCTTACTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTAAAAAACGCACAGCGCACAAAGCCa                                                                                 >  1:3264168/1‑149 (MQ=255)
                                                                    tataACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCAc                                                                                <  2:1187317/149‑1 (MQ=255)
                                                                    tataACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCAc                                                                                <  2:2853036/149‑1 (MQ=255)
                                                                    tataACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCAc                                                                                <  2:3653360/149‑1 (MQ=255)
                                                                      taACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTGGGGTATAAAACTCACCGCGCACAAAGCCCCGt                                                                              >  2:1227608/1‑149 (MQ=255)
                                                                           cgTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCt                                                                         <  2:1454556/149‑1 (MQ=255)
                                                                           cgTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCt                                                                         <  2:779529/149‑1 (MQ=255)
                                                                           cgTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCGTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCt                                                                         <  2:1060625/149‑1 (MQ=255)
                                                                            gTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTg                                                                        >  2:1721487/1‑149 (MQ=255)
                                                                              aaCAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGa                                                                      >  1:2450394/1‑149 (MQ=255)
                                                                              aaCAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGa                                                                      >  1:2450365/1‑149 (MQ=255)
                                                                                 aaGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGc                                                                   <  1:3775024/149‑1 (MQ=255)
                                                                                  aGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGcc                                                                  >  2:1700901/1‑149 (MQ=255)
                                                                                    cACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGGCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGa                                                                >  2:924366/1‑149 (MQ=255)
                                                                                       tGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGt                                                             >  2:3497084/1‑149 (MQ=255)
                                                                                          aaaaaaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGct                                                          <  2:2627553/149‑1 (MQ=255)
                                                                                             aaaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCtctc                                                       >  2:3542139/1‑149 (MQ=255)
                                                                                              aaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTcc                                                      >  1:2777716/1‑149 (MQ=255)
                                                                                              aaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTcc                                                      <  2:808572/149‑1 (MQ=255)
                                                                                              aaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTcc                                                      <  2:257303/149‑1 (MQ=255)
                                                                                              aaCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTcc                                                      <  1:3949327/149‑1 (MQ=255)
                                                                                               aCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTTTGATTACGCATCACCCCTTACTGAACGGTCGAGGCAACCACTACGCCGTCTTGGGTATAAAACTCACCCCGCCCAAAGCCACCTGCGCTGGCCGCCGAGGTGCGCTCCa                                                     >  1:1403501/1‑149 (MQ=255)
                                                                                               aCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCa                                                     <  2:349281/149‑1 (MQ=255)
                                                                                               aCAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCa                                                     >  1:3061491/1‑149 (MQ=255)
                                                                                                  gCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGc                                                  <  1:678497/149‑1 (MQ=255)
                                                                                                   ccGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCt                                                 <  1:1493230/149‑1 (MQ=255)
                                                                                                    cGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCta                                                <  2:457000/149‑1 (MQ=255)
                                                                                                    cGGACGGGTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCta                                                <  2:1186911/149‑1 (MQ=255)
                                                                                                     ggACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCtat                                               <  2:3561310/149‑1 (MQ=255)
                                                                                                      gACGGTTTTCACCTCCGGCTATTTTTTTAATTTTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCtata                                              >  1:790613/1‑149 (MQ=255)
                                                                                                      gACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCtata                                              >  1:3368179/1‑149 (MQ=255)
                                                                                                      gACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCtata                                              <  1:381772/149‑1 (MQ=255)
                                                                                                        cGGTTGGCACCTCCGGCTCTTTTTTTAATTGTGCTTACGCAGCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTCTACACCTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATAca                                            <  2:2155172/149‑1 (MQ=255)
                                                                                                          gTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATAcagc                                          >  2:1637139/1‑149 (MQ=255)
                                                                                                          gTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATAcagc                                          <  1:468107/149‑1 (MQ=255)
                                                                                                           ttttCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATAcagca                                         >  2:2066638/1‑149 (MQ=255)
                                                                                                             ttCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACagcagc                                       <  2:1979819/149‑1 (MQ=255)
                                                                                                             ttCACCTCCGGCTATTTTTTTAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACagcagc                                       <  1:1316836/149‑1 (MQ=255)
                                                                                                               cACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCca                                     <  1:1943542/149‑1 (MQ=255)
                                                                                                                aCCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcat                                    >  1:597982/1‑149 (MQ=255)
                                                                                                                aCCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcat                                    <  1:2371797/149‑1 (MQ=255)
                                                                                                                 ccTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcatt                                   <  2:2332902/149‑1 (MQ=255)
                                                                                                                    ccGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcattcat                                <  2:1789322/149‑1 (MQ=255)
                                                                                                                    ccGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcattcat                                <  2:2760322/149‑1 (MQ=255)
                                                                                                                    ccGGCTATTTTTTTAATTGTGATTACGCATAACCCCGTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCcattcat                                <  2:387684/149‑1 (MQ=255)
                                                                                                                      ggCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCattcattc                              >  2:2602774/1‑149 (MQ=255)
                                                                                                                        cTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaa                            >  2:185668/1‑149 (MQ=255)
                                                                                                                        cTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaa                            >  1:378665/1‑149 (MQ=255)
                                                                                                                        cTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaa                            >  1:2762743/1‑149 (MQ=255)
                                                                                                                        cTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaa                            >  2:608220/1‑149 (MQ=255)
                                                                                                                         tATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAANACTCACCGCGCACAAAGCCACGTGCACTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaaa                           >  1:2083290/1‑149 (MQ=255)
                                                                                                                         tATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACATCAGCCATTCATTCaaa                           >  1:1628303/1‑149 (MQ=255)
                                                                                                                         tATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaaa                           >  2:3120834/1‑149 (MQ=255)
                                                                                                                         tATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGAAGCCATTCATTCaaa                           <  2:1798760/149‑1 (MQ=255)
                                                                                                                         tATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGCACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCaaa                           <  1:2960495/149‑1 (MQ=255)
                                                                                                                           tttttttAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAAtt                         <  1:1598798/149‑1 (MQ=255)
                                                                                                                           tttttttAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAAtt                         >  2:1783182/1‑149 (MQ=255)
                                                                                                                                 tAATTGTGATTACGCATAACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACgg                   <  2:1288065/149‑1 (MQ=255)
                                                                                                                                  aaTTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACggg                  >  1:112870/1‑149 (MQ=255)
                                                                                                                                  aaTTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACggg                  <  1:923473/149‑1 (MQ=255)
                                                                                                                                    ttGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCg                >  2:2518508/1‑149 (MQ=255)
                                                                                                                                      gtgATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGAt              <  2:2811085/149‑1 (MQ=255)
                                                                                                                                       tgATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATg             <  2:3170145/149‑1 (MQ=255)
                                                                                                                                       tgATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATg             >  1:2548074/1‑149 (MQ=255)
                                                                                                                                       tgATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATg             >  1:3213695/1‑149 (MQ=255)
                                                                                                                                        gATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATgg            >  2:229923/1‑149 (MQ=255)
                                                                                                                                         aTTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGa           >  1:978612/1‑149 (MQ=255)
                                                                                                                                          ttACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGaa          >  2:2222396/1‑149 (MQ=255)
                                                                                                                                          ttACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGaa          >  2:1948782/1‑149 (MQ=255)
                                                                                                                                           tACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAAc         <  2:209483/149‑1 (MQ=255)
                                                                                                                                               cATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACcaca     <  1:1034846/149‑1 (MQ=255)
                                                                                                                                                 tCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAATACACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGACAATCATTAAAATTAAACGGGCGATGTAACCACATg   >  1:625449/1‑149 (MQ=255)
                                                                                                                                                 tCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATg   >  1:2259716/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGGTCTGCGTGTACCGGAAGAGCAGGAGCGAGAAGGGCTGGATGTCAACAGCCACGGCGAGAATGCCTATAACGCGTAACAAGCACTGCAAAAAACAGCCGGACGGTTTTCACCTCCGGCTATTTTTTTAATTGTGATTACGCATCACCCCTTCCTGAACGGTCGAGGCAACCAGTACGCCGTCTTGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGG  >  CP009273/469629‑469923

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: