Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 535,589 G→T 5.9% M329I (ATG→ATT allB → allantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273535,5890GT 5.9% 1292.3 / 38.5 337M329I (ATG→ATTallBallantoinase
Reads supporting (aligned to +/- strand):  ref base G (160/157);  new base T (8/12);  total (168/169)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 342 total aligned reads displayed.

AATCGGTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCAAAAAGGCCG  >  CP009273/535442‑535738
                                                                                                                                                   |                                                                                                                                                     
aaTCGGTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTACCGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGaa                                                                                                                                                     <  1:3351683/150‑1 (MQ=255)
aaTCGGTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGaa                                                                                                                                                     >  2:17326/1‑150 (MQ=255)
   cGGTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGc                                                                                                                                                  >  2:2247150/1‑150 (MQ=255)
     gTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCAt                                                                                                                                                <  1:3520323/150‑1 (MQ=255)
     gTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCAt                                                                                                                                                <  1:675263/150‑1 (MQ=255)
      tACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTCCCTGGTTTCCGACCACTCTCCATGCCCGCAGGAAATGAAAGCCGGTAACATCATGAAAGCATg                                                                                                                                               >  2:523801/1‑150 (MQ=255)
       aCTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATgg                                                                                                                                              >  1:5653672/1‑150 (MQ=255)
          ctGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGc                                                                                                                                           >  2:1422380/1‑150 (MQ=255)
              cGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTa                                                                                                                                       <  2:4577412/150‑1 (MQ=255)
              cGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTa                                                                                                                                       <  1:3134396/150‑1 (MQ=255)
               gAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGACATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTAt                                                                                                                                      >  2:2689424/1‑150 (MQ=255)
               gAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTAt                                                                                                                                      >  2:5432512/1‑150 (MQ=255)
               gAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAAAGGTGAAATCGACTGCCTGGTTTCCGACAACTCTCCATGCCCGACGGAAATGAAAGCCGGTAACATAATGAAAGCATGGGGCGGTAt                                                                                                                                      >  2:2806187/1‑150 (MQ=255)
                aaGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATc                                                                                                                                     <  2:1134354/150‑1 (MQ=255)
                     ttCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg                                                                                                                                >  1:2090186/1‑150 (MQ=255)
                         ccgccgATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTg                                                                                                                            <  1:1073258/150‑1 (MQ=255)
                         ccgccgATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTg                                                                                                                            <  1:4534721/150‑1 (MQ=255)
                          cgccgATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGc                                                                                                                           <  1:2932599/150‑1 (MQ=255)
                           gccgATCCGCGATCCGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCa                                                                                                                          <  1:4683358/150‑1 (MQ=255)
                            ccgATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCaa                                                                                                                         >  2:233010/1‑150 (MQ=255)
                             cgATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCaaa                                                                                                                        <  2:3453251/150‑1 (MQ=255)
                              gATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAg                                                                                                                       <  1:780107/150‑1 (MQ=255)
                               aTCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGc                                                                                                                      >  2:2584754/1‑150 (MQ=255)
                                    cgATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCAt                                                                                                                 >  2:83478/1‑150 (MQ=255)
                                    cgATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCAt                                                                                                                 >  1:5455520/1‑150 (MQ=255)
                                           gAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTg                                                                                                          >  1:2406154/1‑150 (MQ=255)
                                            aaaaCCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGa                                                                                                         <  2:4054583/150‑1 (MQ=255)
                                                 cAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTc                                                                                                    <  2:2031005/150‑1 (MQ=255)
                                                  aGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCg                                                                                                   <  2:799789/150‑1 (MQ=255)
                                                      aGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGa                                                                                               <  1:2138692/150‑1 (MQ=255)
                                                        gCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAg                                                                                             <  1:4552616/150‑1 (MQ=255)
                                                         cATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGc                                                                                            >  1:5175856/1‑150 (MQ=255)
                                                         cATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGc                                                                                            <  2:1042597/150‑1 (MQ=255)
                                                         cATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGc                                                                                            <  2:3262999/150‑1 (MQ=255)
                                                         cATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGc                                                                                            <  2:4424238/150‑1 (MQ=255)
                                                           tgtgGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCgg                                                                                          >  2:5546362/1‑150 (MQ=255)
                                                               ggAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACa                                                                                      >  1:3708115/1‑150 (MQ=255)
                                                                gAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAg                                                                                     <  2:4693660/150‑1 (MQ=255)
                                                                 aaaaaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCTTGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGGGATGTTCGATGAAGAGGTACTGa                                                                                    >  2:3027961/1‑150 (MQ=255)
                                                                 aaaaaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGa                                                                                    >  2:2395243/1‑150 (MQ=255)
                                                                   aaaCTGTTTAACGGTGAAATCGACTGCCTGGTTACCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGaaa                                                                                  >  1:3813368/1‑150 (MQ=255)
                                                                    aaCTGTTTAACGTTGAACTCGTCTGCCTGGTTTCCGACCACTCTCCATGCCCGCCTGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAc                                                                                 <  2:2924985/150‑1 (MQ=255)
                                                                    aaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAc                                                                                 <  1:2317945/150‑1 (MQ=255)
                                                                    aaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAc                                                                                 >  1:1358964/1‑150 (MQ=255)
                                                                     aCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAcg                                                                                >  1:2811131/1‑150 (MQ=255)
                                                                       tGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAcgcg                                                                              <  2:2517647/150‑1 (MQ=255)
                                                                       tGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAcgcg                                                                              <  1:2470221/150‑1 (MQ=255)
                                                                       tGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAAcgcg                                                                              <  2:3309869/150‑1 (MQ=255)
                                                                         tttAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGACACGCGGt                                                                            >  2:4402037/1‑150 (MQ=255)
                                                                            aaCGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATg                                                                         >  2:2171483/1‑150 (MQ=255)
                                                                              cGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGtc                                                                       <  2:2088323/150‑1 (MQ=255)
                                                                                 tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTctct                                                                    <  1:3558447/150‑1 (MQ=255)
                                                                                 tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTctct                                                                    <  2:3236719/150‑1 (MQ=255)
                                                                                 tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTctct                                                                    <  2:3760806/150‑1 (MQ=255)
                                                                                 tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTctct                                                                    <  2:4641013/150‑1 (MQ=255)
                                                                                 tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGAATGGACGTGATGTTCGATGAAGCGGTGCAGAAACGCGGTATGTctct                                                                    >  2:4861079/1‑150 (MQ=255)
                                                                                  gAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTg                                                                   >  1:3213848/1‑150 (MQ=255)
                                                                                   aaaTCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGc                                                                  <  2:1509761/150‑1 (MQ=255)
                                                                                   aaaTCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGc                                                                  >  2:892315/1‑150 (MQ=255)
                                                                                       cGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAAt                                                              >  2:238318/1‑150 (MQ=255)
                                                                                        gACGGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATg                                                             >  2:978237/1‑150 (MQ=255)
                                                                                           tGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTc                                                          >  1:1474430/1‑150 (MQ=255)
                                                                                           tGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTc                                                          <  1:1533396/150‑1 (MQ=255)
                                                                                            gCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCg                                                         <  2:1502119/150‑1 (MQ=255)
                                                                                             ccTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCng                                                        <  1:111734/150‑1 (MQ=255)
                                                                                             ccTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCgg                                                        >  1:482936/1‑150 (MQ=255)
                                                                                             ccTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCgg                                                        <  1:1165479/150‑1 (MQ=255)
                                                                                             ccTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCgg                                                        <  1:25529/150‑1 (MQ=255)
                                                                                              cTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGc                                                       >  1:5379504/1‑150 (MQ=255)
                                                                                              cTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGc                                                       <  1:3021075/150‑1 (MQ=255)
                                                                                               tGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCa                                                      >  1:2480184/1‑150 (MQ=255)
                                                                                               tGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAAGGTTCGGCa                                                      >  2:3069399/1‑150 (MQ=255)
                                                                                                 gTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCaaa                                                    <  2:3227874/150‑1 (MQ=255)
                                                                                                 gTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCaaa                                                    >  1:1865987/1‑150 (MQ=255)
                                                                                                  tttCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAAt                                                   >  2:4524131/1‑150 (MQ=255)
                                                                                                    tCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTa                                                 <  1:1617108/150‑1 (MQ=255)
                                                                                                     ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTaa                                                <  1:4440151/150‑1 (MQ=255)
                                                                                                       gACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATg                                              >  2:1301506/1‑150 (MQ=255)
                                                                                                        aCCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATgg                                             >  2:1154752/1‑150 (MQ=255)
                                                                                                                  aTGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACgc                                   >  1:1899846/1‑150 (MQ=255)
                                                                                                                     cccGCCGGAAATGAAAGCCGGTAACATCCTGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACgcagc                                <  1:3989127/150‑1 (MQ=255)
                                                                                                                     cccGCCGGAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACgcagc                                <  2:4334883/150‑1 (MQ=255)
                                                                                                                      ccgccgGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACgcagca                               <  2:5040737/150‑1 (MQ=255)
                                                                                                                        gccgGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGa                             >  1:5059792/1‑150 (MQ=255)
                                                                                                                          cgGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGata                           <  1:5254637/150‑1 (MQ=255)
                                                                                                                            gAAATGAAAGCCGGTAACATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATAtt                         >  2:1447907/1‑150 (MQ=255)
                                                                                                                                 gAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCgg                    >  1:1376528/1‑150 (MQ=255)
                                                                                                                                  aaaGCCGTTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGt                   <  2:3172565/150‑1 (MQ=255)
                                                                                                                                  aaaGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGt                   <  2:2773932/150‑1 (MQ=255)
                                                                                                                                  aaaGACGTTAACATCTTGAAAGCATGGGGAGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGTCTAACGCAGCAGATATTTTCGGt                   <  2:5581644/150‑1 (MQ=255)
                                                                                                                                   aaGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTc                  >  2:5169556/1‑150 (MQ=255)
                                                                                                                                       cGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTgca              >  2:4912930/1‑150 (MQ=255)
                                                                                                                                       cGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTgca              >  1:495436/1‑150 (MQ=255)
                                                                                                                                       cGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTgca              >  2:3682278/1‑150 (MQ=255)
                                                                                                                                           aaCATCATTAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGATGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCaa          >  1:4279748/1‑150 (MQ=255)
                                                                                                                                             catcatGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCaaaa        <  1:3293036/150‑1 (MQ=255)
                                                                                                                                             catcatGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCaaaa        <  1:556122/150‑1 (MQ=255)
                                                                                                                                              atcatGAAAGCCAGGGGCGGTATCGCCGGTCAGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCATATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCaaaaa       <  2:3888916/150‑1 (MQ=255)
                                                                                                                                                  tGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCAAAAAGGcc   <  2:4265862/150‑1 (MQ=255)
                                                                                                                                                   gAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCAAAAAGGcca  <  2:355724/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATCGGTACTCTGGCGAAGTGTTCACCGCCGATCCGCGATCTGGAAAACCAGAAAGGCATGTGGGAAAAACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAGCTGCATGGACGTGATGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAACGCAGCAGATATTTTCGGTCTGCAGCAAAAAGGCCG  >  CP009273/535442‑535738

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: