Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,474,526 C→A 11.9% D415Y (GAC→TAC)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,474,5260CA11.9% 1360.7 / 100.7 389D415Y (GAC→TAC) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (158/184);  new base A (31/15);  total (189/200)
Fisher's exact test for biased strand distribution p-value = 7.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 395 total aligned reads displayed.

CGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGT  >  CP009273/4474377‑4474675
                                                                                                                                                     |                                                                                                                                                     
cgcCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTc                                                                                                                                                       >  2:976834/1‑150 (MQ=255)
  cccAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCAc                                                                                                                                                     >  2:1644341/1‑150 (MQ=255)
  cccAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAAc                                                                                                                                                     >  2:56394/1‑150 (MQ=255)
           aCTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGtt                                                                                                                                            <  2:2478476/150‑1 (MQ=255)
            cTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTc                                                                                                                                           >  1:638687/1‑150 (MQ=255)
              tACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAg                                                                                                                                         <  1:250341/150‑1 (MQ=255)
              tACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAg                                                                                                                                         <  2:2187512/150‑1 (MQ=255)
                    tAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCgg                                                                                                                                   <  1:3471485/150‑1 (MQ=255)
                     aaaGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGc                                                                                                                                  >  2:862191/1‑150 (MQ=255)
                       aGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCaa                                                                                                                                <  2:4023896/150‑1 (MQ=255)
                        gCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAAt                                                                                                                               <  1:4056370/150‑1 (MQ=255)
                            gACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGt                                                                                                                           <  1:3543447/150‑1 (MQ=255)
                              cAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTc                                                                                                                         >  1:3645817/1‑150 (MQ=255)
                                ggAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCAt                                                                                                                       >  1:4654851/1‑150 (MQ=255)
                                 gAAGCAACCTGCGGTAATCGCCCCACCAGGACTAACGCCAATGTTCGCCTTATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCTTCACCCAGCGGTAAGCTCCATGCGCGGTTACCGGATTGTTCAGACGCGTCAATCAGTTCATg                                                                                                                      <  2:806558/150‑1 (MQ=255)
                                 gAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATg                                                                                                                      <  2:4968912/150‑1 (MQ=255)
                                 gAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATg                                                                                                                      >  1:1511052/1‑150 (MQ=255)
                                    gCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGc                                                                                                                   >  1:1754956/1‑150 (MQ=255)
                                    gCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGc                                                                                                                   >  2:3090282/1‑150 (MQ=255)
                                    gCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGc                                                                                                                   <  1:1069936/150‑1 (MQ=255)
                                      aaCCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCa                                                                                                                 <  2:3326715/150‑1 (MQ=255)
                                      aaCCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCa                                                                                                                 <  2:1896843/150‑1 (MQ=255)
                                          tGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCgg                                                                                                             >  1:2923390/1‑150 (MQ=255)
                                           gCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGa                                                                                                            >  1:4272914/1‑150 (MQ=255)
                                            cGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGAt                                                                                                           >  1:3201651/1‑150 (MQ=255)
                                             ggTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCACTCAGTTCATGGGCCAGCGGAtt                                                                                                          >  1:2315086/1‑150 (MQ=255)
                                             ggTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGAtt                                                                                                          >  2:1261329/1‑150 (MQ=255)
                                             ggTAATCGCCCCACCAGGACGACCGCCAAGGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCACTCAGTTCATGGGCCAGCGGAtt                                                                                                          <  2:2974459/150‑1 (MQ=255)
                                              gTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTa                                                                                                         >  2:2378862/1‑150 (MQ=255)
                                                   cGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGtt                                                                                                    <  2:137465/150‑1 (MQ=255)
                                                   cGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGtt                                                                                                    >  2:4628563/1‑150 (MQ=255)
                                                     ccccACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGAGTATGGTTCg                                                                                                  >  2:3443279/1‑150 (MQ=255)
                                                      cccACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGc                                                                                                 >  2:3724295/1‑150 (MQ=255)
                                                      cccACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGc                                                                                                 >  2:3909475/1‑150 (MQ=255)
                                                         accaGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCAt                                                                                              <  1:1002742/150‑1 (MQ=255)
                                                           caGGACGACCGCCAATGTTCGCCATATAGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCa                                                                                            <  1:1287715/150‑1 (MQ=255)
                                                             ggACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGa                                                                                          <  1:4808028/150‑1 (MQ=255)
                                                               acgacCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGAcc                                                                                        >  2:1504991/1‑150 (MQ=255)
                                                                cgacCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCa                                                                                       <  1:4520591/150‑1 (MQ=255)
                                                                 gacCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAg                                                                                      <  1:4730843/150‑1 (MQ=255)
                                                                  acCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCTCCGGTTTGTTCAGACGCGGCAATCAGTTCATGTGCCAGCGGATTATGGTTCGCCATCAGACCAGt                                                                                     <  2:4123401/150‑1 (MQ=255)
                                                                    cGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTaa                                                                                   <  1:4363371/150‑1 (MQ=255)
                                                                    cGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTaa                                                                                   <  1:661810/150‑1 (MQ=255)
                                                                     gCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAat                                                                                  <  1:1039891/150‑1 (MQ=255)
                                                                      ccAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAata                                                                                 >  1:4903197/1‑150 (MQ=255)
                                                                      ccAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAata                                                                                 >  1:5111911/1‑150 (MQ=255)
                                                                       cAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAatat                                                                                >  1:3269046/1‑150 (MQ=255)
                                                                       cAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAatat                                                                                <  1:1609278/150‑1 (MQ=255)
                                                                        aaTGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATatg                                                                               >  2:1782787/1‑150 (MQ=255)
                                                                              cGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAcc                                                                         <  1:1851860/150‑1 (MQ=255)
                                                                              cGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAcc                                                                         <  2:4876579/150‑1 (MQ=255)
                                                                              cGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAcc                                                                         >  1:3972811/1‑150 (MQ=255)
                                                                               gCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGAccc                                                                        >  2:4413825/1‑150 (MQ=255)
                                                                                ccATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCa                                                                       >  1:2396548/1‑150 (MQ=255)
                                                                                 cATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAg                                                                      <  2:1194445/150‑1 (MQ=255)
                                                                                 cATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAg                                                                      >  1:1948459/1‑150 (MQ=255)
                                                                                  atatCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAgc                                                                     >  2:4211504/1‑150 (MQ=255)
                                                                                      cGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGa                                                                 >  1:1831330/1‑150 (MQ=255)
                                                                                      cGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGa                                                                 >  2:4502584/1‑150 (MQ=255)
                                                                                          aaaaTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCGTTCAGACCAGTAATATGATGACCGTGCGCGAtgtc                                                             <  2:4741561/150‑4 (MQ=255)
                                                                                            aaTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGc                                                           >  2:1159613/1‑150 (MQ=255)
                                                                                            aaTTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGc                                                           >  1:2989707/1‑150 (MQ=255)
                                                                                                    gccAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCg                                                   <  2:1038295/149‑1 (MQ=255)
                                                                                                       aGTTGTTCCTGATACTCGTCACCCAGCTGTAAGCGACATGCGAGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCTATCACGAATTCACCGGta                                                <  2:4346955/150‑2 (MQ=255)
                                                                                                          tgttCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGc                                             <  2:1967051/150‑1 (MQ=255)
                                                                                                           gttCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCg                                            >  2:3792417/1‑150 (MQ=255)
                                                                                                            ttCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGt                                           >  1:3315885/1‑150 (MQ=255)
                                                                                                            ttCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATAATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGt                                           >  2:1846659/1‑150 (MQ=255)
                                                                                                             tCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTc                                          <  2:1544450/150‑1 (MQ=255)
                                                                                                                  aTACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCAc                                     >  1:3111234/1‑150 (MQ=255)
                                                                                                                     cTCGTCACCCAGCGGTAAGCGCCAGGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTc                                  <  1:4675127/150‑1 (MQ=255)
                                                                                                                      tCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCa                                 >  2:2477449/1‑150 (MQ=255)
                                                                                                                        gTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAAt                               >  2:282199/1‑150 (MQ=255)
                                                                                                                         tCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATc                              <  2:1850799/150‑1 (MQ=255)
                                                                                                                           aCCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCAc                            <  2:1834564/150‑1 (MQ=255)
                                                                                                                           aCCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCAc                            <  2:777146/150‑1 (MQ=255)
                                                                                                                            cccAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCAcc                           <  2:2687611/150‑1 (MQ=255)
                                                                                                                            cccAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCAcc                           >  2:4591629/1‑150 (MQ=255)
                                                                                                                             ccAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGAAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCAcag                          >  1:1351858/1‑148 (MQ=255)
                                                                                                                             ccAGCGGTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCg                          >  2:1745955/1‑150 (MQ=255)
                                                                                                                               aGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCt                        >  1:2799885/1‑150 (MQ=255)
                                                                                                                                 cGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTc                      <  1:258490/150‑1 (MQ=255)
                                                                                                                                  ggTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTcc                     <  2:3654885/150‑1 (MQ=255)
                                                                                                                                  ggTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTcc                     <  1:5803/150‑1 (MQ=255)
                                                                                                                                   gTAAGCGCCATGCGCGGTAACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCg                    <  1:4101620/150‑1 (MQ=255)
                                                                                                                                      atcgcCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCt                 >  1:2694397/3‑150 (MQ=255)
                                                                                                                                          ccATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCaaa             <  2:638701/150‑1 (MQ=255)
                                                                                                                                           cATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCaaaa            >  2:3599294/1‑150 (MQ=255)
                                                                                                                                           cATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCaaaa            >  1:3875999/1‑150 (MQ=255)
                                                                                                                                           cATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCaaaa            <  1:3727332/150‑1 (MQ=255)
                                                                                                                                            aTGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAAc           <  2:2654595/150‑1 (MQ=255)
                                                                                                                                             tGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACg          <  1:862191/150‑1 (MQ=255)
                                                                                                                                               cgcgGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCt        >  1:1010224/1‑150 (MQ=255)
                                                                                                                                               cgcgGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCt        <  2:3405832/150‑1 (MQ=255)
                                                                                                                                                gcgGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTc       <  1:2191181/150‑1 (MQ=255)
                                                                                                                                                 cgGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCa      <  1:3378670/150‑1 (MQ=255)
                                                                                                                                                   gTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAAc    >  2:5059482/1‑150 (MQ=255)
                                                                                                                                                   gTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAAc    <  1:957438/150‑1 (MQ=255)
                                                                                                                                                    tCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACg   >  1:2017810/1‑150 (MQ=255)
                                                                                                                                                     cACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGt  >  2:3475142/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATTGGACTCCAGTTGTTCCTGATACTCGTCACCCAGCGGTAAGCGCCATGCGCGGTCACCGGATTGTTCAGACGCGGCAATCAGTTCATGGGCCAGCGGATTATGGTTCGCCATCAGACCAGTAATATGATGACCCAGCGCGATCACGCAGGCACCGGTCAGCGTCGCCACGTCAATCACCGCTTCCGGCTCAAAACGCTCAACGT  >  CP009273/4474377‑4474675

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: