Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,640,268 A→C 12.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,2680AC12.2% 903.2 / 123.1 342K39Q (AAG→CAG) ydfVQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (140/160);  new base C (30/12);  total (170/172)
Fisher's exact test for biased strand distribution p-value = 2.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 347 total aligned reads displayed.

CTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTG  >  CP009273/1640122‑1640417
                                                                                                                                                  |                                                                                                                                                     
cTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGt                                                                                                                                                    >  2:228252/1‑150 (MQ=255)
cTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGt                                                                                                                                                    <  1:2942433/150‑1 (MQ=255)
           ttACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTAt                                                                                                                                         <  1:2241971/150‑1 (MQ=255)
             acacCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATtg                                                                                                                                       >  1:5323468/1‑150 (MQ=255)
                    ttGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGt                                                                                                                                >  1:1187585/1‑150 (MQ=255)
                     tGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTc                                                                                                                               <  2:2277846/150‑1 (MQ=255)
                           tACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCGTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGATCg                                                                                                                         >  2:4726981/1‑150 (MQ=255)
                             cAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGtt                                                                                                                       <  1:2732717/150‑1 (MQ=255)
                                aTGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGttttt                                                                                                                    >  2:1018507/1‑150 (MQ=255)
                                 tGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTa                                                                                                                   <  2:3481411/150‑1 (MQ=255)
                                  gTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTaa                                                                                                                  <  2:5420586/150‑1 (MQ=255)
                                  gTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTaa                                                                                                                  <  1:1849446/150‑1 (MQ=255)
                                   tCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTaaa                                                                                                                 <  1:2338859/150‑1 (MQ=255)
                                     aTTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAaga                                                                                                               >  2:3762310/1‑150 (MQ=255)
                                       tACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                             >  1:3314533/1‑150 (MQ=255)
                                       tACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                             <  2:5517991/150‑1 (MQ=255)
                                       tACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGt                                                                                                             <  1:4338667/150‑1 (MQ=255)
                                        aCAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTc                                                                                                            >  1:5252581/1‑150 (MQ=255)
                                        aCAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTc                                                                                                            >  1:2528252/1‑150 (MQ=255)
                                         cAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                           >  1:4267441/1‑150 (MQ=255)
                                         cAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                           >  1:2082357/1‑150 (MQ=255)
                                         cAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTcc                                                                                                           >  1:5078386/1‑150 (MQ=255)
                                          aaGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCg                                                                                                          <  2:2428869/150‑1 (MQ=255)
                                           aGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCgg                                                                                                         >  1:5038439/1‑150 (MQ=255)
                                            gTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGc                                                                                                        >  2:42912/1‑150 (MQ=255)
                                             tATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGcc                                                                                                       <  2:5644396/150‑1 (MQ=255)
                                                cGCACTACAACATGCTTAGGGCAAGTCACGAAGGTAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatc                                                                                                    <  2:5064904/150‑1 (MQ=255)
                                                  cACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatcat                                                                                                  <  2:1905747/150‑1 (MQ=255)
                                                  cACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCcatcat                                                                                                  >  2:1700495/1‑150 (MQ=255)
                                                   aCTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCatcatc                                                                                                 >  1:4704173/1‑150 (MQ=255)
                                                   aCTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCatcatc                                                                                                 >  2:4057484/1‑150 (MQ=255)
                                                   aCTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCatcatc                                                                                                 >  1:2286217/1‑150 (MQ=255)
                                                          caTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTgg                                                                                          <  2:141447/150‑1 (MQ=255)
                                                           aTGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGt                                                                                         >  1:135976/1‑150 (MQ=255)
                                                             gCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTAc                                                                                       >  2:4137758/1‑150 (MQ=255)
                                                             gCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTAc                                                                                       <  2:2622394/150‑1 (MQ=255)
                                                              cTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACa                                                                                      >  1:1567078/1‑150 (MQ=255)
                                                               ttAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAg                                                                                     <  1:2697762/150‑1 (MQ=255)
                                                                     cAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCa                                                                               >  1:3453763/1‑150 (MQ=255)
                                                                     cAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGAAGGTCAGCTAGATTTTAAAGAGTCCGGCCATAATCTTACTGGTACAGACACCa                                                                               >  2:2005650/1‑150 (MQ=255)
                                                                      aaGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCat                                                                              >  2:2677898/1‑150 (MQ=255)
                                                                       aGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCata                                                                             <  2:5379805/150‑1 (MQ=255)
                                                                        gTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGAGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATAGTGCAGGGCAGCGCGATTTTAAAGAGACCGGCCAACATCGTACTGGGACAGccaccagag                                                                            >  2:3044248/1‑147 (MQ=255)
                                                                          cACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATAc                                                                          >  2:5180298/1‑150 (MQ=255)
                                                                             gaaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACttt                                                                       >  2:3171792/1‑150 (MQ=255)
                                                                             gaaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACttt                                                                       <  2:832452/150‑1 (MQ=255)
                                                                              aaggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTtg                                                                      <  2:3059702/150‑1 (MQ=255)
                                                                               aggaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTtgt                                                                     >  1:36138/1‑150 (MQ=255)
                                                                                 gaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGa                                                                   <  1:28509/150‑1 (MQ=255)
                                                                                 gaagTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGa                                                                   <  1:782528/150‑1 (MQ=255)
                                                                                      cAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTGACTGGTACAGACACCATATACTTTGTGACGGTc                                                              >  2:5331196/1‑150 (MQ=255)
                                                                                        gAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAg                                                            <  2:5457164/150‑1 (MQ=255)
                                                                                        gAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAg                                                            <  2:3882178/150‑1 (MQ=255)
                                                                                          aaGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGc                                                          <  1:5506585/150‑1 (MQ=255)
                                                                                             tagtCGTAAGAACGGTGATCTCTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTAc                                                       <  2:4561082/150‑1 (MQ=255)
                                                                                             tagtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTAc                                                       >  2:3636594/1‑150 (MQ=255)
                                                                                               gtCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACat                                                     <  2:1603049/150‑1 (MQ=255)
                                                                                                 cGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACatat                                                   <  2:4689980/150‑1 (MQ=255)
                                                                                                 cGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACatat                                                   <  1:972694/150‑1 (MQ=255)
                                                                                                    aaGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGATCGTTTTAAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGca                                                >  2:2858849/1‑150 (MQ=255)
                                                                                                    aaGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGca                                                <  2:4367632/150‑1 (MQ=255)
                                                                                                             gATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCAcaactcaac                                       <  1:4313011/150‑1 (MQ=255)
                                                                                                               tCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACtt                                     <  1:1057527/150‑1 (MQ=255)
                                                                                                                cACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTa                                    >  2:492492/1‑150 (MQ=255)
                                                                                                                 aCTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAt                                   <  1:4684095/150‑1 (MQ=255)
                                                                                                                  cTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAtt                                  <  1:4796447/150‑1 (MQ=255)
                                                                                                                  cTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAtt                                  <  2:1210600/150‑1 (MQ=255)
                                                                                                                     tcagCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCAt                               <  2:2705940/147‑1 (MQ=255)
                                                                                                                     tCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCAt                               <  1:867594/150‑1 (MQ=255)
                                                                                                                      ccGCTTTGTGCCAGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATc                              >  1:372567/1‑150 (MQ=255)
                                                                                                                      ccGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTACAGAGTCCGGCCATCATCTTACTGGTACAGAAACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATc                              <  2:919058/150‑1 (MQ=255)
                                                                                                                           ttGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAttt                         <  2:3142103/150‑1 (MQ=255)
                                                                                                                            tgtgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttt                        <  2:1160559/150‑1 (MQ=255)
                                                                                                                            tgtgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttt                        <  1:3049177/150‑1 (MQ=255)
                                                                                                                             gtgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAttttt                       <  2:3198397/150‑1 (MQ=255)
                                                                                                                              tgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttttt                      >  1:1030109/1‑150 (MQ=255)
                                                                                                                              tgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttttt                      >  1:4136567/1‑150 (MQ=255)
                                                                                                                              tgCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTAtttttt                      >  1:1394285/1‑150 (MQ=255)
                                                                                                                                ccAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTCCATATGCACAACTCAACTTATTCATCTATTTTTTGc                    >  1:3923421/1‑150 (MQ=255)
                                                                                                                                ccAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACAACATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGc                    <  2:2484738/150‑1 (MQ=255)
                                                                                                                                 cAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCATCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACAACATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCt                   <  2:4244484/150‑1 (MQ=255)
                                                                                                                                 cAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCt                   >  1:1940292/1‑150 (MQ=255)
                                                                                                                                  aGGAGCAGCCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTAGTCATCTATTTTTTgagt                  >  2:4467851/1‑147 (MQ=255)
                                                                                                                                  aGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCtt                  >  2:5207771/1‑150 (MQ=255)
                                                                                                                                  aGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTAAAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCtt                  <  2:2176650/150‑1 (MQ=255)
                                                                                                                                   ggAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCttt                 <  2:5171822/150‑1 (MQ=255)
                                                                                                                                    gAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTa                >  1:5216665/1‑150 (MQ=255)
                                                                                                                                       cagcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCa             <  1:5528464/150‑1 (MQ=255)
                                                                                                                                        agcCATTGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt            >  2:3315745/1‑150 (MQ=255)
                                                                                                                                        agcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt            <  2:2082357/150‑1 (MQ=255)
                                                                                                                                        agcCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCAt            <  1:1337544/150‑1 (MQ=255)
                                                                                                                                           cATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTCTTCATCTATTTTTTGCTTTAGCATGTc         >  2:4400345/1‑150 (MQ=255)
                                                                                                                                           cATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTc         >  2:130296/1‑150 (MQ=255)
                                                                                                                                             ttGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  1:2656751/1‑150 (MQ=255)
                                                                                                                                             ttGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAg       >  1:4618297/1‑150 (MQ=255)
                                                                                                                                              tGCTCAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgt      >  2:390900/1‑150 (MQ=255)
                                                                                                                                               gCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAgtg     <  2:1714319/150‑1 (MQ=255)
                                                                                                                                                  aaGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCTGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTg  <  1:4048877/150‑1 (MQ=255)
                                                                                                                                                  aaGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTg  <  1:4520128/150‑1 (MQ=255)
                                                                                                                                                  aaGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTg  >  1:518595/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTACCCATGTCTTACACCTCTTGTAATTACAAATGTCATTACAAGTATCGCACTACAACATGCTTAGGGCAAGTCACGAAGGAAGTCAGAAAGTAGTCGTAAGAACGGTGATCACTGTCCGCTTTGTGCCAGGAGCAGCCATTGCTAAGTCCATCCTGTATTGTGCAGGTCAGCTCGTTTTTAAAGAGTCCGGCCATCATCTTACTGGTACAGACACCATATACTTTGTGACGGTCAGGCTACATATGCACAACTCAACTTATTCATCTATTTTTTGCTTTAGCATGTCAGTGTTG  >  CP009273/1640122‑1640417

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