Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,628,207 G→A 89.3% 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,2070GA89.3% 593.2 / 93.7 345G110E (GGA→GAA) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (147/161);  total (164/181)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 348 total aligned reads displayed.

AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGG  >  CP009273/1628060‑1628353
                                                                                                                                                   |                                                                                                                                                  
    tggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAat                                                                                                                                              <  2:65147/150‑1 (MQ=17)
     ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatg                                                                                                                                             <  2:3097643/150‑1 (MQ=17)
     ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatg                                                                                                                                             <  2:590751/150‑1 (MQ=17)
     ggtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatg                                                                                                                                             <  2:3069695/150‑1 (MQ=17)
      gtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatga                                                                                                                                            >  2:3743863/1‑150 (MQ=17)
       tgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgat                                                                                                                                           <  1:4245781/150‑1 (MQ=17)
         aTCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatga                                                                                                                                         >  2:5135471/1‑150 (MQ=17)
          tCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgat                                                                                                                                        <  1:4489149/150‑1 (MQ=17)
           cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgatg                                                                                                                                       <  1:1212371/150‑1 (MQ=17)
            aaTTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAatgatgatga                                                                                                                                      >  1:1058908/1‑149 (MQ=17)
                 attgTGACAAAACTGGATCGCCGTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTa                                                                                                                                 >  1:187151/1‑150 (MQ=17)
                  ttgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGGGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTaa                                                                                                                                <  2:3135478/150‑1 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                              >  1:4653917/1‑150 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                              <  1:4235047/150‑1 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                              >  1:3088042/1‑150 (MQ=17)
                    gtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATAAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATAGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAtt                                                                                                                              >  2:4979583/1‑150 (MQ=12)
                     tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAttt                                                                                                                             >  2:3819750/1‑150 (MQ=17)
                     tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAttt                                                                                                                             <  1:2272038/150‑1 (MQ=17)
                     tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAAttt                                                                                                                             >  1:1002732/1‑150 (MQ=17)
                      gACAAAACTGGATCGCCTTGGTTGTAATGAAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTc                                                                                                                            <  2:4821803/150‑1 (MQ=17)
                        cAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcag                                                                                                                          >  1:4322943/1‑150 (MQ=17)
                         aaaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagc                                                                                                                         >  2:4574954/1‑150 (MQ=17)
                          aaaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTcagca                                                                                                                        >  2:1874251/1‑150 (MQ=17)
                            aCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGt                                                                                                                      <  2:3799486/150‑1 (MQ=17)
                            aCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGt                                                                                                                      <  1:1478386/150‑1 (MQ=17)
                             caggATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTc                                                                                                                     >  2:1059102/3‑150 (MQ=17)
                              tGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                                                    <  1:2102793/150‑1 (MQ=17)
                                gATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                                                                  >  1:917817/1‑150 (MQ=17)
                                 aTCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTg                                                                                                                 >  2:1137879/1‑150 (MQ=17)
                                      cTTGGTTGTAATGCAATGGATCTCAGGAAAACTGTGGAACAACTGACCGACACAGGTATCAGAGTGCAGTGGTTAGCATTGTGGGGCATTGACATGACCAGTCCAACAGAAAAATTGATGATGCAAGTAATTTCAGCAGTCGCTGAAttt                                                                                                            <  2:5207307/150‑1 (MQ=11)
                                      cTTGGTTGTAATGCAATGGATCTCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTCGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAttt                                                                                                            <  2:914202/150‑1 (MQ=12)
                                      cTTGGTTGTAATGCAATGGAGATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAttt                                                                                                            <  2:1660012/150‑1 (MQ=17)
                                       ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTg                                                                                                           <  2:3338803/150‑1 (MQ=17)
                                         ggTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGaa                                                                                                         <  2:2328832/150‑1 (MQ=17)
                                             gTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACgag                                                                                                     <  2:1615446/150‑1 (MQ=17)
                                              tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACgaga                                                                                                    >  1:2450273/1‑150 (MQ=17)
                                              tAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCAGTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACgaga                                                                                                    >  1:84417/1‑150 (MQ=17)
                                                  gCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtt                                                                                                >  1:1129746/1‑150 (MQ=17)
                                                   cAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                                                               >  2:2837475/1‑150 (MQ=17)
                                                   cAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCttt                                                                                               >  2:141857/1‑150 (MQ=17)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                              <  1:4926298/150‑1 (MQ=17)
                                                    aaTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtttt                                                                                              >  2:3852794/1‑150 (MQ=17)
                                                     aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTa                                                                                             <  2:687104/150‑1 (MQ=17)
                                                      tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTAc                                                                                            >  2:2218558/1‑150 (MQ=32)
                                                       ggATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACt                                                                                           <  1:4007645/150‑1 (MQ=17)
                                                          tatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                        >  1:1339833/1‑150 (MQ=17)
                                                          tatCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGa                                                                                        >  2:583426/1‑150 (MQ=17)
                                                           atCAGGAAAACAGTGGAACACCTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATGTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGaa                                                                                       <  2:182222/150‑1 (MQ=12)
                                                           atCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGaa                                                                                       >  1:1799089/1‑150 (MQ=17)
                                                           atCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGaa                                                                                       >  1:2142208/1‑150 (MQ=17)
                                                             cAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                                                                     <  1:3203573/150‑1 (MQ=17)
                                                                 aaaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                                                 <  1:1685371/150‑1 (MQ=17)
                                                                  aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                                >  1:4279662/1‑150 (MQ=32)
                                                                  aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                                <  2:4886819/150‑1 (MQ=17)
                                                                  aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                                >  2:2365980/1‑150 (MQ=17)
                                                                    aCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAt                                                                              <  1:4965303/150‑1 (MQ=17)
                                                                    aCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAt                                                                              <  1:4536287/150‑1 (MQ=17)
                                                                     cAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                                             >  1:1264944/1‑150 (MQ=32)
                                                                     cAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                                             >  1:600581/1‑150 (MQ=17)
                                                                       gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTcc                                                                           <  2:4307167/150‑1 (MQ=17)
                                                                        tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCg                                                                          <  1:761629/150‑1 (MQ=17)
                                                                               acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGt                                                                   <  2:1458983/150‑1 (MQ=17)
                                                                                   aCCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAaga                                                               >  1:2833744/1‑150 (MQ=17)
                                                                                    ccGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAagag                                                              >  2:3008465/1‑150 (MQ=17)
                                                                                      gAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGcc                                                            >  1:5100712/1‑150 (MQ=17)
                                                                                      gAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGcc                                                            <  2:2603736/150‑1 (MQ=17)
                                                                                       aaaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATCGTAAGAGccc                                                           >  2:4761423/1‑150 (MQ=17)
                                                                                         aCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcgc                                                         <  1:5244029/150‑1 (MQ=17)
                                                                                               aTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCa                                                   <  2:2914828/150‑1 (MQ=17)
                                                                                               aTCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCa                                                   <  2:1248814/150‑1 (MQ=17)
                                                                                                 cAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg                                                 >  2:35562/1‑150 (MQ=32)
                                                                                                 cAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAgg                                                 <  2:274340/150‑1 (MQ=32)
                                                                                                  agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAggg                                                <  2:3207294/150‑1 (MQ=17)
                                                                                                    agTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaa                                              >  1:1994478/1‑150 (MQ=17)
                                                                                                     gTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaaa                                             >  1:3891184/1‑150 (MQ=17)
                                                                                                     gTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGaaa                                             >  1:5004866/1‑150 (MQ=17)
                                                                                                      tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                                            <  1:3104171/150‑1 (MQ=17)
                                                                                                       gCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACg                                           >  2:388527/1‑150 (MQ=17)
                                                                                                         aTTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGtt                                         <  2:2536162/150‑1 (MQ=17)
                                                                                                         aTTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGAAAGTAAGAGCCCGCGCCGCAGGGAAACGtt                                         >  1:2668302/1‑150 (MQ=12)
                                                                                                           tGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGtttt                                       <  2:4295987/150‑1 (MQ=17)
                                                                                                            gCTTAGCATTGGGGGGCATTGACCAGACCAGTCCAACAGAAAAAAATATGATGCAAGTAAATTCAGAAGTCGCTGAATTTTAACGAGACCATTTACTTGAACACACACAAAACCGGATAGTAATAGCCAGCGGCGCAGCGAAACCTTTTg                                      >  2:4858738/1‑150 (MQ=255)
                                                                                                              ttAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGt                                    <  1:3722440/150‑1 (MQ=17)
                                                                                                              ttAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGt                                    >  2:5048301/1‑150 (MQ=17)
                                                                                                                aGCATTGGGGGGCATTGACCTGACCAGTCCAACCGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAGTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCg                                  <  2:2141141/150‑1 (MQ=17)
                                                                                                                 gCATTGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGa                                 <  1:684365/150‑1 (MQ=17)
                                                                                                                     tGGGGGGCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac                             >  2:866287/1‑150 (MQ=18)
                                                                                                                      ggggggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccacc                            >  1:4566646/1‑150 (MQ=21)
                                                                                                                      ggggggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccacc                            >  1:4039943/1‑150 (MQ=21)
                                                                                                                       gggggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCt                           >  1:3690161/1‑150 (MQ=21)
                                                                                                                        ggggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCtg                          <  2:2289674/150‑1 (MQ=18)
                                                                                                                          ggCATTGACCTGACCAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCtgtg                        <  1:345470/150‑1 (MQ=17)
                                                                                                                                gacctgaccAGTCCAACAGAAAATATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAAt                  <  2:2275868/150‑1 (MQ=17)
                                                                                                                                  cctgaccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATga                >  1:263559/1‑150 (MQ=17)
                                                                                                                                   ctgaccAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATgaa               >  2:396984/1‑150 (MQ=32)
                                                                                                                                       ccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGCAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAAc           >  2:1494098/1‑150 (MQ=17)
                                                                                                                                       ccAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAAc           >  2:4104136/1‑150 (MQ=17)
                                                                                                                                        cAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGTGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACa          <  1:3505968/150‑1 (MQ=17)
                                                                                                                                        cAGTCCAACAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACa          <  2:4418402/150‑1 (MQ=17)
                                                                                                                                                cAGAAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAgg  <  1:989178/150‑1 (MQ=17)
                                                                                                                                                   |                                                                                                                                                  
AATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGG  >  CP009273/1628060‑1628353

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