Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,428,220 A→G 8.8% intergenic (‑289/+28) pinR ← / ← ynaE Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,428,2200AG 8.8% 293.1 / 20.7 113intergenic (‑289/+28)pinR/ynaERac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base A (42/61);  new base G (4/6);  total (46/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 118 total aligned reads displayed.

TTTAACTCGGTTTGTTTTTCGAGTTCCGGGCGGACTCAAGGAAGAAGAATAGTGTTGCGTGTTATTTTAACCAGATTTCAAGTTGTTTGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGTTC  >  CP009273/1428094‑1428369
                                                                                                                              |                                                                                                                                                     
tttatgtCGGTTTGTTTTTCGAGTTCCGGGCGGACTCAAGGAAGAAGAATAGTGTTGCGTGTTATTTTAACCAGATTTCAAGTTGTTTGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTttat                                                                                                                                <  1:5155328/144‑1 (MQ=1)
                                                                       cataTTTCAAGTTGTTTGGTCGTGGAAAAGAGGAGCAAAAAGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGATTATTATATTACATTTTAACATATTAATTTTTGTTGTCTTTAAGAAATGATATCTCTTTTTGTAATAATTCTTTg                                                         >  2:2964861/4‑150 (MQ=1)
                                                                                    gtatGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCatg                                            >  1:2089445/4‑150 (MQ=1)
                                                                                       tGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCatgatg                                         <  2:3573182/150‑1 (MQ=11)
                                                                                       tGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCatgatg                                         <  2:3573056/150‑1 (MQ=11)
                                                                                       tGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCatgatg                                         <  1:3106765/150‑1 (MQ=11)
                                                                                       tGGTCGTGGAAAAGTGGAGCAAAATGTTGTNAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCatgatg                                         <  2:2426401/150‑1 (MQ=11)
                                                                                        ggTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGt                                        >  1:3008627/1‑150 (MQ=14)
                                                                                        ggTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGt                                        <  1:3340493/150‑1 (MQ=31)
                                                                                         gTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATACAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtt                                       <  2:3244558/150‑1 (MQ=17)
                                                                                         gTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtt                                       >  1:4289609/1‑150 (MQ=32)
                                                                                         gTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtt                                       <  2:3244635/150‑1 (MQ=32)
                                                                                           cGTGGAACAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTAGATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtttt                                     <  2:3470792/150‑1 (MQ=17)
                                                                                           cGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtttt                                     <  1:4457998/150‑1 (MQ=17)
                                                                                           cGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtttt                                     <  2:4390644/150‑1 (MQ=32)
                                                                                           cGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtttt                                     <  1:1227711/150‑1 (MQ=32)
                                                                                           cGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGtttt                                     <  2:1165994/150‑1 (MQ=32)
                                                                                            gTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTc                                    >  2:3572515/1‑150 (MQ=32)
                                                                                             tGGAAAAGTGGAGCAAAATGTTGTTAAATTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCt                                   <  2:1089586/150‑1 (MQ=17)
                                                                                             tGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCt                                   <  1:3748046/150‑1 (MQ=32)
                                                                                             tGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCt                                   <  1:1206072/150‑1 (MQ=32)
                                                                                             tGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCt                                   >  2:4334064/1‑150 (MQ=32)
                                                                                             tGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCt                                   >  1:2058575/1‑150 (MQ=32)
                                                                                              ggAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCtt                                  <  1:488248/150‑1 (MQ=32)
                                                                                              ggAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCtt                                  <  1:488251/150‑1 (MQ=32)
                                                                                              ggAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCtt                                  <  1:14324/150‑1 (MQ=32)
                                                                                              ggAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCtt                                  >  1:1404795/1‑150 (MQ=32)
                                                                                               gAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAGATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCttt                                 >  1:2395728/1‑150 (MQ=17)
                                                                                               gAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCttt                                 >  1:5042158/1‑150 (MQ=32)
                                                                                               gAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTAATATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTAttt                                 >  1:2834588/1‑150 (MQ=17)
                                                                                                aaaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTa                                >  1:3834435/1‑150 (MQ=32)
                                                                                                aaaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAAAGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTa                                <  1:3158662/150‑1 (MQ=17)
                                                                                                 aaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTAt                               <  2:4911570/150‑1 (MQ=32)
                                                                                                 aaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTAt                               <  1:547174/150‑1 (MQ=32)
                                                                                                 aaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTAt                               <  1:768533/150‑1 (MQ=32)
                                                                                                 aaaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTAt                               <  1:1189207/150‑1 (MQ=32)
                                                                                                  aaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATg                              >  1:3212632/1‑150 (MQ=32)
                                                                                                  aaGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATg                              >  1:2568809/1‑150 (MQ=32)
                                                                                                   aGTGGAGCAAAATGTTGTTAATGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTATAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGt                             <  1:2973280/150‑1 (MQ=17)
                                                                                                   aGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGt                             <  2:4451584/150‑1 (MQ=32)
                                                                                                   aGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGt                             >  2:2264176/1‑150 (MQ=32)
                                                                                                   aGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGt                             <  1:2973314/150‑1 (MQ=32)
                                                                                                   aGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGAGGTTTTCTTTATGt                             >  1:1303261/1‑150 (MQ=17)
                                                                                                    gTGGAGCACAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGtt                            <  2:4683658/150‑1 (MQ=17)
                                                                                                    gTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGtt                            <  2:3277628/150‑1 (MQ=32)
                                                                                                     tGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGttt                           <  1:4375358/150‑1 (MQ=32)
                                                                                                     tGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGGTTATTATATTGCATTTTACCATTTAAATTTTGGTTGTCTTTAAGACCTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGttt                           <  1:5089398/150‑1 (MQ=12)
                                                                                                     tGGAGCAAAATGTTGATAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGttt                           <  1:4375305/150‑1 (MQ=17)
                                                                                                      ggAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTc                          <  1:4472149/150‑1 (MQ=32)
                                                                                                      ggAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTc                          <  1:3505066/150‑1 (MQ=32)
                                                                                                      ggAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCAAGCCATGATGTTTTCTTTATGTTTc                          >  2:4941876/1‑150 (MQ=17)
                                                                                                       gAGCAACATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTcc                         <  2:643662/150‑1 (MQ=17)
                                                                                                       gAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTcc                         <  1:4572436/150‑1 (MQ=17)
                                                                                                       gAGCAAAATGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTcc                         <  2:4692417/150‑1 (MQ=17)
                                                                                                        aGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCt                        >  1:1471278/1‑150 (MQ=32)
                                                                                                        aGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCt                        >  1:1471244/1‑150 (MQ=32)
                                                                                                        aGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCt                        <  2:1359509/150‑1 (MQ=32)
                                                                                                        aGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCAATGATGTTTTCTTTATGTTTCCt                        >  2:2775041/1‑150 (MQ=17)
                                                                                                          cAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTc                      <  2:655256/150‑1 (MQ=32)
                                                                                                          cAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTc                      >  2:3771484/1‑150 (MQ=32)
                                                                                                          cAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTc                      >  2:3771701/1‑150 (MQ=32)
                                                                                                            aaaTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCaa                    >  1:1133733/1‑150 (MQ=32)
                                                                                                             aaTTTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   <  2:759374/150‑1 (MQ=17)
                                                                                                             aaTGTTGTTAAAGTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   >  2:4175112/1‑150 (MQ=17)
                                                                                                             aaTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   <  2:1085712/150‑1 (MQ=32)
                                                                                                             aaTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   >  2:2011962/1‑150 (MQ=32)
                                                                                                             aaTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   <  1:428701/150‑1 (MQ=32)
                                                                                                             aaTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAAt                   <  2:759357/150‑1 (MQ=32)
                                                                                                              aTGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATg                  <  1:4440185/150‑1 (MQ=32)
                                                                                                               tgttgtTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGt                 <  2:4083693/150‑1 (MQ=32)
                                                                                                               tgttgtTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGt                 <  2:4083678/150‑1 (MQ=32)
                                                                                                               tgttgtTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGt                 <  1:2969298/150‑1 (MQ=32)
                                                                                                               tgttgtTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGt                 >  2:3087214/1‑150 (MQ=32)
                                                                                                                gttgttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTa                >  1:1584317/1‑150 (MQ=32)
                                                                                                                gttgttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTa                >  1:3070179/1‑150 (MQ=32)
                                                                                                                gttgttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGGTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTa                >  2:4642927/1‑150 (MQ=17)
                                                                                                                 ttgttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTaa               <  2:4091350/150‑1 (MQ=32)
                                                                                                                   gttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATc             >  1:249605/1‑150 (MQ=32)
                                                                                                                    ttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCa            <  2:3465929/150‑1 (MQ=32)
                                                                                                                    ttAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTGTGTAATAATTCTTTGTTATCCAGCCATGATGGTTTCTTTATGTTTCCGTCTATGTAATCa            >  1:2710736/1‑150 (MQ=255)
                                                                                                                     tAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCaa           >  1:1874489/1‑150 (MQ=32)
                                                                                                                     tAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCaa           >  1:1874314/1‑150 (MQ=32)
                                                                                                                      aaaGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAg          >  1:4963274/1‑150 (MQ=32)
                                                                                                                      aaaGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAg          <  2:2276989/150‑1 (MQ=32)
                                                                                                                       aaGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGc         <  1:1549599/150‑1 (MQ=32)
                                                                                                                       aaGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGc         <  2:715196/150‑1 (MQ=32)
                                                                                                                        aGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCa        <  2:4598861/150‑1 (MQ=32)
                                                                                                                         gTGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAAtcaatcaa       >  1:3153315/1‑150 (MQ=17)
                                                                                                                          tGGAGAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCCTGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAAt      >  2:1366962/1‑150 (MQ=17)
                                                                                                                          tGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAAt      <  2:1599948/150‑1 (MQ=32)
                                                                                                                           ggAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATg     >  1:1149663/1‑150 (MQ=32)
                                                                                                                           ggAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATg     >  2:432847/1‑150 (MQ=32)
                                                                                                                           ggAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATg     >  1:3501361/1‑150 (MQ=32)
                                                                                                                           ggAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTATTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATg     <  1:318070/150‑1 (MQ=17)
                                                                                                                            gAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGt    <  1:3595609/150‑1 (MQ=32)
                                                                                                                            gAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGt    >  1:3268363/1‑150 (MQ=32)
                                                                                                                             aaaaaTGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTTTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGtt   >  2:3898194/1‑150 (MQ=17)
                                                                                                                             aaaaaTGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGtt   >  1:1704504/1‑150 (MQ=32)
                                                                                                                             aaaaaTGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGtt   >  1:4405523/1‑150 (MQ=32)
                                                                                                                              aaaaTGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGTTc  >  2:2341266/1‑150 (MQ=32)
                                                                                                                              |                                                                                                                                                     
TTTAACTCGGTTTGTTTTTCGAGTTCCGGGCGGACTCAAGGAAGAAGAATAGTGTTGCGTGTTATTTTAACCAGATTTCAAGTTGTTTGGTCGTGGAAAAGTGGAGCAAAATGTTGTTAAAGTGGAAAAATGATAAAAAAGTAAGTTTATTATATTACATTTTACCATTTAAATTTTGGTTGTCTTTAAGAACTGATATCGCTGTTTGTAATAATTCTTTGTTATCCAGCCATGATGTTTTCTTTATGTTTCCTTCAATGTAATCAAGCAATGTTC  >  CP009273/1428094‑1428369

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