Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,628,207 G→A 86.6% 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,2070GA86.6% 400.2 / 90.6 254G110E (GGA→GAA) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (20/14);  new base A (113/107);  total (133/121)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 257 total aligned reads displayed.

AAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAC  >  CP009273/1628059‑1628350
                                                                                                                                                    |                                                                                                                                               
aaaTGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGa                                                                                                                                                >  1:1636282/1‑150 (MQ=32)
   tgtgGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAaaaa                                                                                                                                             >  2:3284233/1‑150 (MQ=17)
    gtggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGaaaaaa                                                                                                                                            >  1:4200075/1‑150 (MQ=17)
     tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAat                                                                                                                                           <  2:151827/150‑1 (MQ=17)
     tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAat                                                                                                                                           >  2:4429304/1‑150 (MQ=17)
        tgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgat                                                                                                                                        >  1:3393901/1‑150 (MQ=17)
        tgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgat                                                                                                                                        >  1:1712048/1‑150 (MQ=17)
         gATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatg                                                                                                                                       <  1:4412309/150‑1 (MQ=17)
          aTCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatga                                                                                                                                      <  2:3056734/150‑1 (MQ=17)
          aTCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatga                                                                                                                                      <  2:577657/150‑1 (MQ=17)
                   ttgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTaa                                                                                                                             <  2:31650/150‑1 (MQ=17)
                    tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAt                                                                                                                            <  1:2314357/150‑1 (MQ=17)
                     gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                           <  2:2746223/150‑1 (MQ=17)
                      tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAttt                                                                                                                          >  2:657050/1‑150 (MQ=17)
                       gACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTc                                                                                                                         >  2:1008235/1‑150 (MQ=32)
                       gACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTc                                                                                                                         >  1:4033713/1‑150 (MQ=17)
                       gACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTc                                                                                                                         >  2:4029572/1‑150 (MQ=17)
                        aCAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTca                                                                                                                        <  1:4439996/150‑1 (MQ=17)
                        aCAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTca                                                                                                                        <  1:3204222/150‑1 (MQ=17)
                          aaaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagc                                                                                                                      <  2:3835230/150‑1 (MQ=17)
                           aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                     <  2:1577744/150‑1 (MQ=17)
                             aCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGt                                                                                                                   >  2:1695054/1‑150 (MQ=17)
                               tGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                 <  1:1757731/150‑1 (MQ=32)
                               tGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                 <  1:3731662/150‑1 (MQ=32)
                               tGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                 <  1:1925788/150‑1 (MQ=17)
                               tGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                 <  2:789758/150‑1 (MQ=17)
                                ggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGc                                                                                                                >  1:1572819/1‑150 (MQ=17)
                                 gATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                                                               >  2:3301327/1‑150 (MQ=17)
                                   tCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGa                                                                                                             >  1:1674662/1‑150 (MQ=17)
                                    cGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGaa                                                                                                            >  1:3315587/1‑150 (MQ=17)
                                     gCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAt                                                                                                           <  2:2105721/150‑1 (MQ=17)
                                      ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAtt                                                                                                          >  1:4356379/1‑150 (MQ=17)
                                         tGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGa                                                                                                       >  2:2265308/1‑150 (MQ=17)
                                           gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAAc                                                                                                     >  2:2870247/1‑150 (MQ=32)
                                           gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAAc                                                                                                     <  2:2427021/150‑1 (MQ=17)
                                           gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAAc                                                                                                     <  2:3410918/150‑1 (MQ=17)
                                            ttGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                                                    <  1:876935/150‑1 (MQ=17)
                                                aaTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAc                                                                                                <  2:1317401/150‑1 (MQ=17)
                                                aaTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGAc                                                                                                <  1:1154724/150‑1 (MQ=17)
                                                   gCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtt                                                                                             >  2:2753623/1‑150 (MQ=32)
                                                     aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                           >  1:1253996/1‑150 (MQ=32)
                                                     aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                           <  1:2062063/150‑1 (MQ=17)
                                                     aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTAGCTGAATTTGAACGAGACCtttt                                                                                           >  1:2044378/1‑150 (MQ=17)
                                                       tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTAc                                                                                         >  1:3338319/1‑150 (MQ=17)
                                                       tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTAc                                                                                         <  1:1300757/150‑1 (MQ=17)
                                                             tCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                                                   >  2:1028699/1‑150 (MQ=17)
                                                             tCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                                                   <  1:1014386/150‑1 (MQ=17)
                                                              cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                                                                  >  1:323586/1‑150 (MQ=17)
                                                              cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                                                                  >  1:2719358/1‑150 (MQ=17)
                                                               aGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                                                 <  1:1218529/150‑1 (MQ=17)
                                                                 gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                                               <  2:1984579/150‑1 (MQ=17)
                                                                 gAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                                               <  2:3602118/150‑1 (MQ=17)
                                                                   aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                             >  1:4040262/1‑150 (MQ=17)
                                                                    aaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCa                                                                            <  2:1087767/150‑1 (MQ=17)
                                                                      cAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                                          <  2:365894/150‑1 (MQ=17)
                                                                       aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGTATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                         >  1:431891/1‑150 (MQ=17)
                                                                       aGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                         >  1:2906841/1‑150 (MQ=17)
                                                                         tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCg                                                                       >  2:2389639/1‑150 (MQ=17)
                                                                          ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                                         <  2:3445260/147‑1 (MQ=17)
                                                                           gAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCggg                                                                     <  2:2072681/150‑1 (MQ=32)
                                                                            aacaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                                                    >  1:2438900/1‑150 (MQ=17)
                                                                              caacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATa                                                                  >  1:4353894/1‑150 (MQ=32)
                                                                              caacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATa                                                                  >  2:470863/1‑150 (MQ=17)
                                                                                acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGt                                                                >  1:663145/1‑150 (MQ=17)
                                                                                   gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAag                                                             <  1:107628/150‑1 (MQ=17)
                                                                                       gAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGcc                                                         <  1:1452606/150‑1 (MQ=17)
                                                                                       gAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGcc                                                         >  1:3585675/1‑150 (MQ=17)
                                                                                         aaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                                       <  1:721704/150‑1 (MQ=17)
                                                                                          aCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcgc                                                      <  1:3524494/150‑1 (MQ=17)
                                                                                          aCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcgc                                                      >  2:2940315/1‑150 (MQ=17)
                                                                                            aGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcgg                                                    >  2:900846/1‑150 (MQ=17)
                                                                                            aGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcgg                                                    >  1:3029926/1‑150 (MQ=17)
                                                                                                aTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCa                                                >  1:765197/1‑150 (MQ=17)
                                                                                                      gTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaaa                                          <  2:1381741/150‑1 (MQ=17)
                                                                                                        gCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACg                                        <  1:2657703/150‑1 (MQ=32)
                                                                                                           ttGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCCGGGAAACGttt                                     >  1:199625/1‑150 (MQ=17)
                                                                                                            tGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGtttt                                    >  2:4066273/1‑150 (MQ=17)
                                                                                                             gCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTg                                   <  1:1259415/150‑1 (MQ=17)
                                                                                                              cTTAGCATTGGGGGGCATTGACCTTACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTAGCTGAATTTGAACGAGACATTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg                                  >  2:751185/1‑150 (MQ=255)
                                                                                                              cTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAATTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGACACGTTGTgg                                  >  1:4257584/1‑150 (MQ=12)
                                                                                                               ttAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGt                                 <  2:4034002/150‑1 (MQ=17)
                                                                                                                tAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTc                                <  2:3107762/150‑1 (MQ=17)
                                                                                                                tAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTc                                >  1:2396166/1‑150 (MQ=17)
                                                                                                                 aGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCg                               <  1:589848/150‑1 (MQ=17)
                                                                                                                  gCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGa                              >  2:4107340/1‑150 (MQ=17)
                                                                                                                     ttGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGacca                           <  1:4499976/150‑1 (MQ=17)
                                                                                                                             cATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                   >  2:2190929/1‑150 (MQ=32)
                                                                                                                             cATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                   >  2:4337200/1‑150 (MQ=17)
                                                                                                                             cATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                   <  1:4136499/150‑1 (MQ=17)
                                                                                                                             cATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGtt                   >  2:2060169/1‑150 (MQ=17)
                                                                                                                              aTTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTa                  <  1:3130514/150‑1 (MQ=17)
                                                                                                                              aTTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTa                  <  1:3939756/150‑1 (MQ=17)
                                                                                                                                  acctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATg              >  1:3530520/1‑150 (MQ=17)
                                                                                                                                  acctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATg              >  1:3762242/1‑150 (MQ=17)
                                                                                                                                   cctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATga             <  1:1641946/150‑1 (MQ=17)
                                                                                                                                     tgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATgaag           <  2:838198/150‑1 (MQ=17)
                                                                                                                                        ccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAAc        >  2:436611/1‑150 (MQ=17)
                                                                                                                                        ccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAAc        >  2:192548/1‑150 (MQ=17)
                                                                                                                                         cAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACa       <  2:1627998/150‑1 (MQ=17)
                                                                                                                                              cAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAc  >  1:2831389/1‑150 (MQ=17)
                                                                                                                                                    |                                                                                                                                               
AAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAAC  >  CP009273/1628059‑1628350

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: