Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,640,268 A→C 6.2% K39Q (AAG→CAG)  ydfV → Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,640,2680AC 6.2% 855.5 / 43.7 260K39Q (AAG→CAG) ydfVQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (124/120);  new base C (4/12);  total (128/132)
Fisher's exact test for biased strand distribution p-value = 6.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 263 total aligned reads displayed.

TGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTT  >  CP009273/1640120‑1640416
                                                                                                                                                    |                                                                                                                                                    
    aCCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTc                                                                                                                                                  <  2:1176114/149‑1 (MQ=255)
     cccATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTcc                                                                                                                                                 >  1:4892754/1‑149 (MQ=255)
        aTGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATc                                                                                                                                              <  1:2104869/149‑1 (MQ=255)
             ttACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTa                                                                                                                                         <  2:2107311/149‑1 (MQ=255)
             ttACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTa                                                                                                                                         >  2:5021505/1‑149 (MQ=255)
              tACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTAt                                                                                                                                        >  2:2589543/1‑149 (MQ=255)
               acacCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTAtt                                                                                                                                       <  2:3656025/149‑1 (MQ=255)
                 acCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATtgt                                                                                                                                     >  1:1260215/1‑149 (MQ=255)
                 acCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATtgt                                                                                                                                     <  2:736940/149‑1 (MQ=255)
                       tGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGt                                                                                                                               <  1:3917903/149‑1 (MQ=255)
                        gTAATTACAAATGTCATTACAAGTATCGCACTACGACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTc                                                                                                                              >  2:2082997/1‑149 (MQ=255)
                          aaTTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAg                                                                                                                            >  2:3872899/1‑149 (MQ=255)
                           aTTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGc                                                                                                                           <  1:4231328/149‑1 (MQ=255)
                           aTTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGc                                                                                                                           <  2:2949825/149‑1 (MQ=255)
                            ttACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCt                                                                                                                          <  1:2445821/149‑1 (MQ=255)
                             tACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTc                                                                                                                         >  2:304190/1‑149 (MQ=255)
                             tACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTc                                                                                                                         >  1:1958836/1‑149 (MQ=255)
                              aCAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCg                                                                                                                        >  1:4406470/1‑149 (MQ=255)
                               cAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGt                                                                                                                       <  1:640183/149‑1 (MQ=255)
                               cAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGt                                                                                                                       <  1:3422350/149‑1 (MQ=255)
                                  aTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGtttt                                                                                                                    >  2:3333104/1‑149 (MQ=255)
                                    gTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTa                                                                                                                  <  1:2536750/149‑1 (MQ=255)
                                    gTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTa                                                                                                                  <  1:2258758/149‑1 (MQ=255)
                                     tCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTaa                                                                                                                 >  2:615950/1‑149 (MQ=255)
                                         tACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAagag                                                                                                             <  1:3443838/149‑1 (MQ=255)
                                          aCAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                            >  2:1976080/1‑149 (MQ=255)
                                          aCAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                            >  2:5245519/1‑149 (MQ=255)
                                          aCAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                            <  2:2399891/149‑1 (MQ=255)
                                           cAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTc                                                                                                           <  2:2635659/149‑1 (MQ=255)
                                            aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                          >  1:1170487/1‑149 (MQ=255)
                                            aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                          >  1:328556/1‑149 (MQ=255)
                                                  cGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcat                                                                                                    >  2:3417728/1‑149 (MQ=255)
                                                  cGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcat                                                                                                    >  2:4687964/1‑149 (MQ=255)
                                                    cACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatca                                                                                                  >  1:4000744/1‑149 (MQ=255)
                                                       tACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCt                                                                                               <  1:305033/149‑1 (MQ=255)
                                                        acaacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCtt                                                                                              >  1:328777/1‑149 (MQ=255)
                                                        acaacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCtt                                                                                              >  1:3149984/1‑149 (MQ=255)
                                                        acaacaTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCtt                                                                                              >  1:16243/1‑149 (MQ=255)
                                                              tGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGt                                                                                        <  1:2666318/149‑1 (MQ=255)
                                                              tGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGt                                                                                        <  1:4341540/149‑1 (MQ=255)
                                                               gCTTATGGAAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTa                                                                                       <  1:5133467/149‑1 (MQ=255)
                                                               gCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTa                                                                                       <  1:181235/149‑1 (MQ=255)
                                                                   aGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGa                                                                                   <  1:5076998/149‑1 (MQ=255)
                                                                        aaGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCa                                                                              >  1:4584406/1‑149 (MQ=255)
                                                                         aGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCat                                                                             <  2:828158/149‑1 (MQ=255)
                                                                         aGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCat                                                                             >  2:2321177/1‑149 (MQ=255)
                                                                            cACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCAtata                                                                          <  1:3804947/149‑1 (MQ=255)
                                                                             aCGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATAc                                                                         <  1:766174/149‑1 (MQ=255)
                                                                             aCGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATAc                                                                         >  1:3767819/1‑149 (MQ=255)
                                                                             aCGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATAc                                                                         >  1:1850251/1‑149 (MQ=255)
                                                                               gaaggaagTCAGAAAGTAGTCGTGAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACtt                                                                       <  1:478882/149‑1 (MQ=255)
                                                                                 aggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTtg                                                                     <  1:349348/149‑1 (MQ=255)
                                                                                  ggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTtgt                                                                    <  2:2754017/149‑1 (MQ=255)
                                                                                  ggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTtgt                                                                    >  2:4087250/1‑149 (MQ=255)
                                                                                     agTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGAc                                                                 <  1:1422647/149‑1 (MQ=255)
                                                                                     agTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGAc                                                                 >  1:3531741/1‑149 (MQ=255)
                                                                                     agTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGAc                                                                 >  1:3019171/1‑149 (MQ=255)
                                                                                         aGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTc                                                             <  1:865320/149‑1 (MQ=255)
                                                                                          gAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCa                                                            >  2:5224005/1‑149 (MQ=255)
                                                                                           aaaGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAg                                                           >  1:5087179/1‑149 (MQ=255)
                                                                                            aaGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTTTGACGGTCAgg                                                          >  2:4405694/1‑149 (MQ=255)
                                                                                             agtagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGc                                                         >  2:768265/1‑149 (MQ=255)
                                                                                              gtagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCt                                                        <  2:495605/149‑1 (MQ=255)
                                                                                               tagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTa                                                       >  2:3571524/1‑149 (MQ=255)
                                                                                               tagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTa                                                       >  1:4138713/1‑149 (MQ=255)
                                                                                                 gtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACa                                                     >  1:4455616/1‑149 (MQ=255)
                                                                                                 gtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACa                                                     <  2:1022249/149‑1 (MQ=255)
                                                                                                 gtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACa                                                     >  2:1415173/1‑149 (MQ=255)
                                                                                                   cGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACata                                                   <  2:5087179/149‑1 (MQ=255)
                                                                                                   cGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACata                                                   <  2:2698140/149‑1 (MQ=255)
                                                                                                    gTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACatat                                                  >  2:2201088/1‑149 (MQ=255)
                                                                                                     tAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATg                                                 >  2:2471274/1‑149 (MQ=255)
                                                                                                          aCGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCATGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACAGGTCCAGACACCATATACTTTGTGACGGTCAGGCTACATATGAAcaa                                            >  2:2278960/1‑149 (MQ=255)
                                                                                                           cGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaac                                           >  2:1136251/1‑149 (MQ=255)
                                                                                                               gATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactcaa                                       >  1:1540482/1‑149 (MQ=255)
                                                                                                                 tCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactcaact                                     <  2:1340966/149‑1 (MQ=255)
                                                                                                                 tCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactcaact                                     >  2:1615422/1‑149 (MQ=255)
                                                                                                                   aCTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTa                                   >  2:2660080/1‑149 (MQ=255)
                                                                                                                     tGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAtt                                 <  1:1016692/149‑1 (MQ=255)
                                                                                                                     tGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAtt                                 <  1:1689735/149‑1 (MQ=255)
                                                                                                                     tGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAtt                                 <  1:4054281/149‑1 (MQ=255)
                                                                                                                             ttGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtt                         <  1:2586980/149‑1 (MQ=255)
                                                                                                                                tgCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAttttt                      <  1:415532/149‑1 (MQ=255)
                                                                                                                                tgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAttttt                      <  1:4719887/149‑1 (MQ=255)
                                                                                                                                 gCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttttt                     <  1:2898488/149‑1 (MQ=255)
                                                                                                                                  ccAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTg                    >  1:540797/1‑149 (MQ=255)
                                                                                                                                    aGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCt                  >  2:3788643/1‑149 (MQ=255)
                                                                                                                                     ggAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCtt                 >  2:451660/1‑149 (MQ=255)
                                                                                                                                      gAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCttt                <  1:2289602/149‑1 (MQ=255)
                                                                                                                                       agcagcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTa               <  1:2231417/149‑1 (MQ=255)
                                                                                                                                           gcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt           >  1:3845877/1‑149 (MQ=255)
                                                                                                                                            ccATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATg          >  1:5002225/1‑149 (MQ=255)
                                                                                                                                             cATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGt         <  1:332167/149‑1 (MQ=255)
                                                                                                                                             cATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGt         <  2:701048/149‑1 (MQ=255)
                                                                                                                                             cATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGt         >  2:1225977/1‑149 (MQ=255)
                                                                                                                                               ttGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCa       <  1:1413332/149‑1 (MQ=255)
                                                                                                                                               ttGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCa       >  2:4811806/1‑149 (MQ=255)
                                                                                                                                                 gCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgt     <  1:2695841/149‑1 (MQ=255)
                                                                                                                                                   tAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgtgt   <  2:3355525/149‑1 (MQ=255)
                                                                                                                                                    aaGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGAATGTCAGTGtt  >  2:4047606/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTT  >  CP009273/1640120‑1640416

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: