Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,373,129 C→A 21.7% D332Y (GAT→TAT)  menF ← isochorismate synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,373,1290CA21.7% 858.9 / 177.1 320D332Y (GAT→TAT) menFisochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base C (105/145);  new base A (37/33);  total (142/178)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

TTCGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCT  >  CP009273/2372980‑2373278
                                                                                                                                                     |                                                                                                                                                     
ttCGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCAcatc                                                                                                                                                       <  1:497241/150‑1 (MQ=255)
 tCGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCAcatca                                                                                                                                                      >  2:3842206/1‑150 (MQ=255)
  cGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCAcatcat                                                                                                                                                     <  1:4936801/150‑1 (MQ=255)
  cGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCAcatcat                                                                                                                                                     >  1:929122/1‑150 (MQ=255)
    cTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATcc                                                                                                                                                   >  2:1287660/1‑150 (MQ=255)
       tgttgtAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCt                                                                                                                                                >  2:716417/1‑150 (MQ=255)
          tgtAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTttg                                                                                                                                             >  2:4078027/1‑150 (MQ=255)
            tAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTttgtt                                                                                                                                           >  2:4634876/1‑150 (MQ=255)
              aTGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGa                                                                                                                                         <  2:1956189/150‑1 (MQ=255)
                    gaTAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAg                                                                                                                                   <  2:3127854/150‑1 (MQ=255)
                      tAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTc                                                                                                                                 <  2:1615738/150‑1 (MQ=255)
                           cGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGAt                                                                                                                            >  1:4024712/1‑150 (MQ=255)
                            gCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATa                                                                                                                           <  2:3849336/150‑1 (MQ=255)
                             cAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATAc                                                                                                                          >  1:1314818/1‑150 (MQ=255)
                             cAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATAc                                                                                                                          >  2:2243437/1‑150 (MQ=255)
                              aGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACa                                                                                                                         >  2:2297916/1‑150 (MQ=255)
                              aGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACa                                                                                                                         <  2:3548474/150‑1 (MQ=255)
                               gAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAg                                                                                                                        <  1:1814253/150‑1 (MQ=255)
                                 aCCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAgcg                                                                                                                      <  2:4608985/150‑1 (MQ=255)
                                 aCCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAgcg                                                                                                                      <  2:3338052/150‑1 (MQ=255)
                                 aCCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAgcg                                                                                                                      <  1:5018962/150‑1 (MQ=255)
                                   cGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAgcggc                                                                                                                    <  1:218304/150‑1 (MQ=255)
                                   cGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAgcggc                                                                                                                    >  1:1687071/1‑150 (MQ=255)
                                    ggCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGACGTCCAGATACAgcggcg                                                                                                                   >  1:3146224/1‑150 (MQ=255)
                                    ggCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAgcggcg                                                                                                                   >  1:2287163/1‑150 (MQ=255)
                                     gCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGa                                                                                                                  <  2:4409090/150‑1 (MQ=255)
                                       gTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAg                                                                                                                <  2:2908087/150‑1 (MQ=255)
                                          ccATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGAtg                                                                                                             <  1:534496/150‑1 (MQ=255)
                                                tcGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAc                                                                                                       <  1:2707246/150‑1 (MQ=255)
                                                 cGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACt                                                                                                      >  1:3645422/1‑150 (MQ=255)
                                                      gAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacg                                                                                                 >  2:3402886/1‑150 (MQ=255)
                                                        aCGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgca                                                                                               <  1:3414024/150‑1 (MQ=255)
                                                         cGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcag                                                                                              >  1:795459/1‑150 (MQ=255)
                                                            ttCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTGGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacg                                                                                           <  1:2544421/150‑1 (MQ=255)
                                                              cGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgca                                                                                         <  1:2002339/150‑1 (MQ=255)
                                                                 gACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                                                      >  1:1386013/1‑150 (MQ=255)
                                                                 gACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                                                      >  1:1385885/1‑150 (MQ=255)
                                                                 gACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                                                      >  2:4691776/1‑150 (MQ=255)
                                                                   cGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTAcc                                                                                    <  2:1990886/150‑1 (MQ=255)
                                                                   cGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTAcc                                                                                    <  1:2415103/150‑1 (MQ=255)
                                                                    gggCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCt                                                                                   <  1:4740856/150‑1 (MQ=255)
                                                                        gATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTgcgg                                                                               <  2:2263740/150‑1 (MQ=255)
                                                                             aCTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTa                                                                          >  2:1634394/1‑150 (MQ=255)
                                                                               tGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaaa                                                                        <  2:4030792/150‑1 (MQ=255)
                                                                                     gcCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATcc                                                                  >  1:4939316/1‑150 (MQ=255)
                                                                                     gcCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATcc                                                                  >  1:4134428/1‑150 (MQ=255)
                                                                                     gcCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATcc                                                                  <  1:2477155/150‑1 (MQ=255)
                                                                                        aGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGc                                                               >  2:4218022/1‑150 (MQ=255)
                                                                                         gATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCg                                                              <  2:2063563/150‑1 (MQ=255)
                                                                                         gATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCg                                                              >  2:2331167/1‑150 (MQ=255)
                                                                                           tCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTc                                                            <  2:847513/150‑1 (MQ=255)
                                                                                            cgcgcgGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCt                                                           >  1:4353249/1‑150 (MQ=255)
                                                                                                cgGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGt                                                       <  2:153106/150‑1 (MQ=255)
                                                                                                 ggTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTa                                                      <  2:380874/150‑1 (MQ=255)
                                                                                                 ggTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTa                                                      <  2:690643/150‑1 (MQ=255)
                                                                                                   tAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATc                                                    <  1:58770/150‑1 (MQ=255)
                                                                                                    aaGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCg                                                   <  2:4937350/150‑1 (MQ=255)
                                                                                                     aGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCgg                                                  <  2:1701053/150‑1 (MQ=255)
                                                                                                     aGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCgg                                                  <  2:1505301/150‑1 (MQ=255)
                                                                                                     aGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCgg                                                  >  1:3948597/1‑150 (MQ=255)
                                                                                                       ccAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGcc                                                >  1:1536068/1‑150 (MQ=255)
                                                                                                        cAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCt                                               <  2:3719123/150‑1 (MQ=255)
                                                                                                         aGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTg                                              <  1:1281728/150‑1 (MQ=255)
                                                                                                          gCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGt                                             <  1:867553/150‑1 (MQ=255)
                                                                                                           cAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATACGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTa                                            >  2:1780627/1‑150 (MQ=255)
                                                                                                            aaCTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCGTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGtca                                           >  1:3331801/1‑148 (MQ=255)
                                                                                                              ctgctgCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATc                                         <  1:4634876/150‑1 (MQ=255)
                                                                                                              ctgctgCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATc                                         >  2:2157335/1‑150 (MQ=255)
                                                                                                                gctgcCGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGt                                       <  1:1605660/150‑1 (MQ=255)
                                                                                                                    ccGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGAc                                   <  2:1190477/150‑1 (MQ=255)
                                                                                                                     cGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACa                                  <  2:4380947/150‑1 (MQ=255)
                                                                                                                     cGTTGGCTGCAACTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACa                                  <  1:3440072/150‑1 (MQ=255)
                                                                                                                      gTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAg                                 <  1:3754376/150‑1 (MQ=255)
                                                                                                                         ggCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGata                              <  2:4025493/150‑1 (MQ=255)
                                                                                                                         ggCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGata                              >  2:318287/1‑150 (MQ=255)
                                                                                                                         ggCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGata                              <  2:3870652/150‑1 (MQ=255)
                                                                                                                          gCTGCAACTGCTGTAAACAGATCACCTCATCCGCTTTGTTGAGTGAAGTCCAGAGACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGatat                             <  2:4857610/150‑1 (MQ=255)
                                                                                                                            tGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCt                           <  2:2812345/150‑1 (MQ=255)
                                                                                                                               aaCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTcc                        >  2:3272668/1‑150 (MQ=255)
                                                                                                                               aaCTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTcc                        >  1:3353036/1‑150 (MQ=255)
                                                                                                                                aCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTcca                       <  2:2990952/150‑1 (MQ=255)
                                                                                                                                aCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTcca                       <  1:2560468/150‑1 (MQ=255)
                                                                                                                                 cTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccac                      >  1:1087063/1‑150 (MQ=255)
                                                                                                                                 cTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccac                      <  1:2823831/150‑1 (MQ=255)
                                                                                                                                 cTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccac                      <  2:4530370/150‑1 (MQ=255)
                                                                                                                                 cTGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccac                      <  2:609315/150‑1 (MQ=255)
                                                                                                                                  tGATGTAAACAGATCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccacc                     >  1:2596609/1‑150 (MQ=255)
                                                                                                                                   gATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccacca                    >  1:3875310/1‑150 (MQ=255)
                                                                                                                                    aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccaccac                   <  2:1766404/150‑1 (MQ=255)
                                                                                                                                     tGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTccaccacc                  >  2:262470/1‑150 (MQ=255)
                                                                                                                                      gTAAACAGATCACATCATCCGCGGTGTTGAGTGAAGTCCAGATACCGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTCTCGGCCTGTAATCGTTGACAGATATCTTccaccccca                 >  2:2661158/1‑150 (MQ=255)
                                                                                                                                       tAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAg                <  1:5074531/150‑1 (MQ=255)
                                                                                                                                           cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATg            >  2:2459919/1‑150 (MQ=255)
                                                                                                                                             gATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGtt          <  1:72245/150‑1 (MQ=255)
                                                                                                                                              aTCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTtc         >  2:92462/1‑150 (MQ=255)
                                                                                                                                              aTCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTtc         >  1:2926109/1‑150 (MQ=255)
                                                                                                                                               tCACATAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTtct        >  2:1414354/1‑150 (MQ=255)
                                                                                                                                                 acaTAATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcg      >  1:372842/1‑150 (MQ=255)
                                                                                                                                                  cataatCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcgc     >  1:4194221/1‑150 (MQ=255)
                                                                                                                                                   atcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcgcg    <  1:4519785/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCGCTTTGTTGTAATGAGAGATAGCCCGCAGAACCGGCGTACCATTCTCGGGTGAACGGTTCGTGACGGGCGATAAACTGTCGCGCCAGATCGCGCGGTAAGCCAGCAACTGCTGCCGTTGGCTGCAACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCT  >  CP009273/2372980‑2373278

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: