Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,628,207 G→A 19.1% G110E (GGA→GAA)  pinQ → Qin prophage; putative site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,628,2070GA19.1% 538.8 / 126.0 231G110E (GGA→GAA) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (91/95);  new base A (22/22);  total (114/117)
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.52e-01

Only 100 of 235 total aligned reads displayed.

AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGC  >  CP009273/1628060‑1628354
                                                                                                                                                   |                                                                                                                                                   
aaTGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGa                                                                                                                                                    <  1:1033685/149‑1 (MQ=32)
aaTGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGa                                                                                                                                                    >  2:2915139/1‑149 (MQ=32)
 aTGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATGGGGGGGCATTGACCTGACCAGTCCAACAGGaa                                                                                                                                                   >  2:2984540/1‑149 (MQ=17)
    tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGaaaaa                                                                                                                                                <  2:2010930/149‑1 (MQ=32)
     ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAat                                                                                                                                               <  2:2142769/149‑1 (MQ=17)
       tgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatga                                                                                                                                             <  1:233574/149‑1 (MQ=17)
          tCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatga                                                                                                                                          <  2:4455969/149‑1 (MQ=32)
           cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAatgatgat                                                                                                                                         <  2:1846240/149‑1 (MQ=32)
              ttgattgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCa                                                                                                                                      >  2:2220397/1‑149 (MQ=32)
               tgattgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCaa                                                                                                                                     <  1:449484/149‑1 (MQ=32)
                   tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTaa                                                                                                                                 <  1:4441719/149‑1 (MQ=32)
                     tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAAtt                                                                                                                               >  1:2038719/1‑149 (MQ=32)
                      gACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAAttt                                                                                                                              >  2:605803/1‑149 (MQ=32)
                         aaaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTcag                                                                                                                           <  2:4505385/149‑1 (MQ=32)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTCTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTcagc                                                                                                                          >  2:3560139/1‑149 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTcagc                                                                                                                          >  2:3745020/1‑149 (MQ=32)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                     >  1:1006697/1‑149 (MQ=32)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                     >  2:3426453/1‑149 (MQ=32)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                     >  1:3496820/1‑149 (MQ=32)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                     >  2:1743599/1‑149 (MQ=17)
                               ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                     >  2:1743803/1‑149 (MQ=17)
                                gATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                                                                    <  1:2453727/149‑1 (MQ=32)
                                gATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                                                                    <  1:1134910/149‑1 (MQ=32)
                                 aTCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                                                                   <  2:1550208/149‑1 (MQ=32)
                                 aTCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                                                                   >  2:3169842/1‑149 (MQ=32)
                                  tCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTg                                                                                                                  >  2:704584/1‑149 (MQ=32)
                                      cTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAtt                                                                                                              <  1:790426/149‑1 (MQ=32)
                                       ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAttt                                                                                                             <  1:2146201/149‑1 (MQ=32)
                                       ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAAGTTCAGCAGTCGCTGAAttt                                                                                                             <  2:3534242/149‑1 (MQ=17)
                                        tGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTg                                                                                                            <  1:470471/149‑1 (MQ=17)
                                         ggTTGGAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGACACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACCGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGa                                                                                                           <  2:1905051/149‑1 (MQ=255)
                                            tGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                                                        <  2:865527/149‑1 (MQ=32)
                                                aTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAc                                                                                                    >  1:3927570/1‑149 (MQ=17)
                                                  gCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCt                                                                                                  <  2:1050658/149‑1 (MQ=32)
                                                  gCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCt                                                                                                  <  2:3369107/149‑1 (MQ=32)
                                                   cAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtt                                                                                                 >  1:1453603/1‑149 (MQ=32)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                                                                >  2:518515/1‑149 (MQ=32)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                                                                <  1:714440/149‑1 (MQ=32)
                                                       ggATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTAc                                                                                             >  1:475992/1‑149 (MQ=32)
                                                       ggATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTAc                                                                                             <  1:1522854/149‑1 (MQ=17)
                                                        gATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACt                                                                                            >  1:1609504/1‑149 (MQ=32)
                                                         atatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACtt                                                                                           <  1:1487753/149‑1 (MQ=32)
                                                          tatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTg                                                                                          >  1:584702/1‑149 (MQ=32)
                                                           atCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                         >  2:2602448/1‑149 (MQ=32)
                                                             cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                                                       <  1:2984540/149‑1 (MQ=32)
                                                             cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                                                       <  2:3864313/149‑1 (MQ=32)
                                                                 aaaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                                                   <  2:3117227/149‑1 (MQ=32)
                                                                  aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                                                  >  2:3340709/1‑149 (MQ=32)
                                                                   aaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                                 >  1:2552073/1‑149 (MQ=32)
                                                                       gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                             <  1:4300646/149‑1 (MQ=32)
                                                                       gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                             >  2:3334446/1‑149 (MQ=17)
                                                                         ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCg                                                                           >  1:4104639/1‑149 (MQ=32)
                                                                          gAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCgg                                                                          <  1:1122025/149‑1 (MQ=32)
                                                                            acaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                                                        <  1:3169842/149‑1 (MQ=32)
                                                                            acaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                                                        >  2:4042141/1‑149 (MQ=32)
                                                                            acaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTGTACTTGAACGCACTCATTCCGGGa                                                                        >  2:3063372/1‑149 (MQ=17)
                                                                            acaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                                                        <  1:825690/149‑1 (MQ=17)
                                                                             caacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGAt                                                                       <  1:704584/149‑1 (MQ=32)
                                                                                  gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTaa                                                                  >  1:4206606/1‑149 (MQ=32)
                                                                                   aCCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGTACTCATTCCGGGATAGTAag                                                                 >  2:4173100/1‑149 (MQ=17)
                                                                                   aCCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGTACTCATTCCGGGATAGTAag                                                                 >  2:3791079/1‑149 (MQ=17)
                                                                                    ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAaga                                                                <  2:4338182/149‑1 (MQ=32)
                                                                                    ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAaga                                                                <  1:1194746/149‑1 (MQ=32)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGcc                                                             >  1:3333048/1‑149 (MQ=32)
                                                                                         aCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                                           >  2:2566754/1‑149 (MQ=32)
                                                                                            ggTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcgg                                                        >  2:450864/1‑149 (MQ=32)
                                                                                               aTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc                                                     <  1:1370685/149‑1 (MQ=32)
                                                                                               aTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc                                                     >  1:4487521/1‑149 (MQ=17)
                                                                                                tCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCa                                                    >  1:3247770/1‑149 (MQ=32)
                                                                                                  agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg                                                  >  1:3150216/1‑149 (MQ=32)
                                                                                                   gagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTGTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAggg                                                 <  2:4387900/149‑1 (MQ=17)
                                                                                                    agTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa                                                >  2:499559/1‑149 (MQ=32)
                                                                                                    agTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa                                                >  2:1039102/1‑149 (MQ=32)
                                                                                                     gTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaa                                               <  2:600488/149‑1 (MQ=32)
                                                                                                      tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaaa                                              <  2:3184094/149‑1 (MQ=17)
                                                                                                      tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaaa                                              >  1:361474/1‑149 (MQ=17)
                                                                                                       gCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                                             >  2:2792650/1‑149 (MQ=32)
                                                                                                       gCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACg                                            >  2:2434175/1‑149 (MQ=17)
                                                                                                         aTTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTCTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGt                                           <  1:3221493/149‑1 (MQ=17)
                                                                                                          ttGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGtt                                          >  2:896847/1‑149 (MQ=32)
                                                                                                           tGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGttt                                         <  2:4296239/149‑1 (MQ=32)
                                                                                                              ttAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGTTGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg                                      <  2:1235867/149‑1 (MQ=17)
                                                                                                              ttAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg                                      <  1:988449/149‑1 (MQ=32)
                                                                                                               tAGCATTGGTGGGCATTGACCTGACCAGTCTAAGAGGAAAATTGATGATGCAAGTAATTTGAGCTGTCGCTGAATTTGTACTAGACCTTTTACTTGAACTCACTGATTACGGGATAGTAAGTGCCCTCGTCTCAGTGAAACGTTTTgtt                                     <  2:1229234/149‑3 (MQ=255)
                                                                                                               tAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGt                                     <  1:901144/149‑1 (MQ=32)
                                                                                                               tAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTCTTGGt                                     <  2:4474125/149‑1 (MQ=17)
                                                                                                                   aTTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGac                                 >  2:3287313/1‑149 (MQ=32)
                                                                                                                    ttGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGacc                                >  1:4295848/1‑149 (MQ=32)
                                                                                                                       gggggCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccacc                             >  2:1370369/1‑149 (MQ=34)
                                                                                                                       gggggCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccacc                             >  1:1262424/1‑149 (MQ=34)
                                                                                                                             aTTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                       <  1:3426453/149‑1 (MQ=32)
                                                                                                                              ttGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTa                      <  2:3999391/149‑1 (MQ=32)
                                                                                                                              ttGACCTGACCAGTCCAAAAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTTAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTa                      <  1:1804966/149‑1 (MQ=12)
                                                                                                                                 acctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCGGTGTTAAAt                   >  2:3504085/1‑149 (MQ=17)
                                                                                                                                             cAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGaaa       <  2:2836155/149‑1 (MQ=32)
                                                                                                                                              aaCAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAc      <  1:2220397/149‑1 (MQ=32)
                                                                                                                                              aaCAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAc      <  1:14673/149‑1 (MQ=32)
                                                                                                                                              aaCAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAc      <  1:3846519/149‑1 (MQ=17)
                                                                                                                                               aCAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACa     >  1:1709176/1‑149 (MQ=32)
                                                                                                                                               aCAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGGGTTAAATGAAGAACAGAAACa     >  2:2684651/1‑149 (MQ=17)
                                                                                                                                                   |                                                                                                                                                   
AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGC  >  CP009273/1628060‑1628354

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: