Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,475,510 G→T 12.0% R87S (CGT→AGT)  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,475,5100GT12.0% 1044.5 / 76.7 300R87S (CGT→AGT) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base G (132/132);  new base T (19/17);  total (151/149)
Fisher's exact test for biased strand distribution p-value = 8.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 305 total aligned reads displayed.

CGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTT  >  CP009273/4475361‑4475659
                                                                                                                                                     |                                                                                                                                                     
cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACg                                                                                                                                                       <  1:3287105/150‑1 (MQ=255)
    tctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCa                                                                                                                                                   >  1:4480251/1‑150 (MQ=255)
       tttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATcc                                                                                                                                                >  1:2458177/1‑150 (MQ=255)
        ttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCa                                                                                                                                               <  2:20754/150‑1 (MQ=255)
        ttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCa                                                                                                                                               >  2:4061782/1‑150 (MQ=255)
            cGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTc                                                                                                                                           >  2:2287314/1‑150 (MQ=255)
             gTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCa                                                                                                                                          <  2:3427956/150‑1 (MQ=255)
              tctcGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCAc                                                                                                                                         <  1:4255541/150‑1 (MQ=255)
              tctcGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCAc                                                                                                                                         >  2:4129912/1‑150 (MQ=255)
                tcGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGt                                                                                                                                       >  1:3687238/1‑150 (MQ=255)
                 cGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGtt                                                                                                                                      >  1:3868318/1‑150 (MQ=255)
                  gACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTc                                                                                                                                     >  1:1771898/1‑150 (MQ=255)
                      gCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTg                                                                                                                                 <  2:2791242/150‑1 (MQ=255)
                          gACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGc                                                                                                                             >  1:2312655/1‑150 (MQ=255)
                            cGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCaa                                                                                                                           <  2:1268320/150‑1 (MQ=255)
                               aCTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCa                                                                                                                        >  1:505110/1‑150 (MQ=255)
                                 tttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAAt                                                                                                                      >  1:957369/1‑150 (MQ=255)
                                 tttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAAt                                                                                                                      <  2:3669145/150‑1 (MQ=255)
                                  ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCTGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                                                                                     <  2:1890722/150‑1 (MQ=255)
                                  ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                                                                                     >  1:4126257/1‑150 (MQ=255)
                                  ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                                                                                     >  1:821162/1‑150 (MQ=255)
                                   tCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATaa                                                                                                                    >  2:1658495/1‑150 (MQ=255)
                                       gtagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGgag                                                                                                                <  2:4070156/150‑1 (MQ=255)
                                        tagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGgaga                                                                                                               <  2:3351868/150‑1 (MQ=255)
                                          gttgttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAAt                                                                                                             <  2:1800851/150‑1 (MQ=255)
                                             gttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCg                                                                                                          >  1:2888049/1‑150 (MQ=255)
                                              ttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGc                                                                                                         <  1:4558485/150‑1 (MQ=255)
                                               tACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCt                                                                                                        <  2:3984432/150‑1 (MQ=255)
                                               tACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCt                                                                                                        <  1:1847800/150‑1 (MQ=255)
                                                 cGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTAGCTCg                                                                                                      <  1:4078729/150‑1 (MQ=255)
                                                 cGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCg                                                                                                      >  1:3471902/1‑150 (MQ=255)
                                                  ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                                                                                     >  2:1438827/1‑150 (MQ=255)
                                                  ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                                                                                     <  2:3575672/150‑1 (MQ=255)
                                                  ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                                                                                     <  1:3734979/150‑1 (MQ=255)
                                                     cTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGaaa                                                                                                  >  1:1322721/1‑150 (MQ=255)
                                                      tttAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAg                                                                                                 <  1:3180642/150‑1 (MQ=255)
                                                         aaCGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTAc                                                                                              <  1:1935719/150‑1 (MQ=255)
                                                         aaCGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTAc                                                                                              <  2:3570259/150‑1 (MQ=255)
                                                           cGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACAt                                                                                            <  1:2959198/150‑1 (MQ=255)
                                                             tGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTc                                                                                          <  1:1013571/150‑1 (MQ=255)
                                                                aGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGa                                                                                       <  2:802750/150‑1 (MQ=255)
                                                                  cTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAAc                                                                                     <  1:546329/150‑1 (MQ=255)
                                                                    cagtcagAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACAt                                                                                   <  1:4450037/150‑1 (MQ=255)
                                                                     agtcagAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATg                                                                                  <  2:2705078/150‑1 (MQ=255)
                                                                     agtcagAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATg                                                                                  <  2:4317051/150‑1 (MQ=255)
                                                                      gtcagAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATgg                                                                                 >  2:2018929/1‑150 (MQ=255)
                                                                        cagAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTg                                                                               <  1:1400588/150‑1 (MQ=255)
                                                                          gAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgca                                                                             >  2:3869455/1‑150 (MQ=255)
                                                                          gAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgca                                                                             >  2:1092229/1‑150 (MQ=255)
                                                                           aaaGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcag                                                                            >  1:3564358/1‑150 (MQ=255)
                                                                           aaaGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcag                                                                            >  1:437673/1‑150 (MQ=255)
                                                                             aGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcagca                                                                          >  1:2757354/1‑150 (MQ=255)
                                                                                aGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGcaaca                                                                       <  1:2068857/150‑1 (MQ=255)
                                                                                    cGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGt                                                                   <  2:3188214/150‑1 (MQ=255)
                                                                                       ttCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTTCATTCGGAACATGGTGCAGCAACAATGTCTg                                                                >  2:2556353/1‑150 (MQ=255)
                                                                                       ttCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTg                                                                >  2:3353126/1‑150 (MQ=255)
                                                                                        tCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGc                                                               <  1:79022/150‑1 (MQ=255)
                                                                                         ccATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTTcc                                                              <  2:4039556/150‑1 (MQ=255)
                                                                                           aTTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGcccc                                                            >  2:2637032/1‑150 (MQ=255)
                                                                                            ttGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCg                                                           >  2:4544193/1‑150 (MQ=255)
                                                                                             tGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCgg                                                          >  1:2068521/1‑150 (MQ=255)
                                                                                              gAGCCAGTATCATACAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAAAGTCTGCCCCGGt                                                         >  2:1018918/1‑150 (MQ=255)
                                                                                               aGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGtt                                                        >  2:1314196/1‑150 (MQ=255)
                                                                                               aGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGtt                                                        >  2:866541/1‑150 (MQ=255)
                                                                                                gCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGttt                                                       <  1:4339221/150‑1 (MQ=255)
                                                                                                 ccAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGtttt                                                      >  1:1561071/1‑150 (MQ=255)
                                                                                                  cAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttc                                                     >  1:2477838/1‑150 (MQ=255)
                                                                                                    gTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttcct                                                   <  2:2993800/150‑1 (MQ=255)
                                                                                                    gTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttcct                                                   >  1:618070/1‑150 (MQ=255)
                                                                                                     tATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttcctt                                                  >  1:667494/1‑150 (MQ=255)
                                                                                                     tATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttcctt                                                  >  1:3646224/1‑150 (MQ=255)
                                                                                                      aTCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTttccttc                                                 >  1:1400176/1‑150 (MQ=255)
                                                                                                        cATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCa                                               <  1:2930284/150‑1 (MQ=255)
                                                                                                        cATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCa                                               <  1:735996/150‑1 (MQ=255)
                                                                                                         aTTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAg                                              >  1:1965750/1‑150 (MQ=255)
                                                                                                           tCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGtt                                            <  2:2165333/150‑1 (MQ=255)
                                                                                                           tCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGtt                                            <  2:3757013/150‑1 (MQ=255)
                                                                                                           tCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGtt                                            <  1:2880342/150‑1 (MQ=255)
                                                                                                             aGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCg                                          >  1:3358685/1‑150 (MQ=255)
                                                                                                                 tATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCTGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccc                                      >  2:509644/1‑150 (MQ=255)
                                                                                                                 tATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGcccc                                      <  2:3709440/150‑1 (MQ=255)
                                                                                                                    tAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGAc                                   >  1:282375/1‑150 (MQ=255)
                                                                                                                       tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTa                                <  2:2085796/150‑1 (MQ=255)
                                                                                                                       tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTa                                >  2:1495819/1‑150 (MQ=255)
                                                                                                                       tGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTa                                <  1:2592350/150‑1 (MQ=255)
                                                                                                                          ttttCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCa                             <  2:1167423/150‑1 (MQ=255)
                                                                                                                            ttCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAAGTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAgg                           >  2:837800/1‑150 (MQ=255)
                                                                                                                             tCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAggg                          >  2:603187/1‑150 (MQ=255)
                                                                                                                              cTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGgc                         >  1:4044537/1‑150 (MQ=255)
                                                                                                                              cTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGgc                         >  2:734848/1‑150 (MQ=255)
                                                                                                                              cTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGgc                         <  1:4292896/150‑1 (MQ=255)
                                                                                                                                   tAACCTGCTTGTACTGACTCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGa                    >  2:3056445/1‑150 (MQ=255)
                                                                                                                                       cTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTa                >  1:2421851/1‑150 (MQ=255)
                                                                                                                                       cTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTa                >  1:111189/1‑150 (MQ=255)
                                                                                                                                         gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTAcc              >  1:12767/1‑150 (MQ=255)
                                                                                                                                          cTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTAccc             <  2:2462764/150‑1 (MQ=255)
                                                                                                                                             gTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATc          <  1:379631/150‑1 (MQ=255)
                                                                                                                                                cTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCt       >  1:3830127/1‑150 (MQ=255)
                                                                                                                                                 tGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTg      <  1:188360/150‑1 (MQ=255)
                                                                                                                                                    cGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGAtt   <  1:3925515/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTT  >  CP009273/4475361‑4475659

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: