Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 552,831 A→T 10.7% F123I (TTC→ATC)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273552,8310AT10.7% 874.7 / 67.3 283F123I (TTC→ATC) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (130/121);  new base T (2/28);  total (133/150)
Fisher's exact test for biased strand distribution p-value = 9.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGG  >  CP009273/552682‑552980
                                                                                                                                                     |                                                                                                                                                     
gcCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGaa                                                                                                                                                       <  1:2464769/150‑1 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATAcc                                                                                                                                                    <  1:4197462/150‑1 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc                                                                                                                                                    >  1:744702/1‑150 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc                                                                                                                                                    >  1:2577529/1‑150 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc                                                                                                                                                    >  2:4356972/1‑150 (MQ=255)
      gATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGACACCGTc                                                                                                                                                 >  2:380190/1‑150 (MQ=255)
       atatTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTcc                                                                                                                                                <  2:2756917/150‑1 (MQ=255)
         atTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCAc                                                                                                                                              <  2:626616/150‑1 (MQ=255)
         atTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGAAGGAAACCGTCCAc                                                                                                                                              >  1:4413185/1‑150 (MQ=255)
            cGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTc                                                                                                                                           <  1:1633688/150‑1 (MQ=255)
             gATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCt                                                                                                                                          <  2:542638/150‑1 (MQ=255)
             gATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCt                                                                                                                                          <  1:2177032/150‑1 (MQ=255)
              aTGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCtt                                                                                                                                         >  2:3166023/1‑150 (MQ=255)
                 cgcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGt                                                                                                                                      <  2:4416715/150‑1 (MQ=255)
                 cgcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGt                                                                                                                                      >  2:2281279/1‑150 (MQ=255)
                 cgcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGt                                                                                                                                      >  2:2974749/1‑150 (MQ=255)
                  gcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTACCGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTc                                                                                                                                     <  1:840254/150‑1 (MQ=255)
                    cAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCg                                                                                                                                   <  1:1057674/150‑1 (MQ=255)
                     aaTCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCgg                                                                                                                                  >  2:2497112/1‑150 (MQ=255)
                      aTCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGa                                                                                                                                 >  2:438876/1‑150 (MQ=255)
                       tCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGAt                                                                                                                                <  2:3971240/150‑1 (MQ=255)
                        caccacGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATg                                                                                                                               >  2:3101556/1‑150 (MQ=255)
                          ccacGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGaa                                                                                                                             <  1:3350726/150‑1 (MQ=255)
                          ccacGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGaa                                                                                                                             <  1:3668863/150‑1 (MQ=255)
                            acGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCTCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATa                                                                                                                           <  2:783725/150‑1 (MQ=255)
                            acGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATa                                                                                                                           <  2:2793121/150‑1 (MQ=255)
                              ggCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACg                                                                                                                         <  1:2402383/150‑1 (MQ=255)
                              ggCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACg                                                                                                                         >  1:4094660/1‑150 (MQ=255)
                              ggCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACg                                                                                                                         >  1:1376451/1‑150 (MQ=255)
                                 gTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTc                                                                                                                      <  2:1670288/150‑1 (MQ=255)
                                 gTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTc                                                                                                                      >  1:3050248/1‑150 (MQ=255)
                                        ccGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACt                                                                                                               <  1:2289095/150‑1 (MQ=255)
                                          gAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACATACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACttt                                                                                                             <  2:1064666/150‑1 (MQ=255)
                                           aaGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTg                                                                                                            <  1:2712948/150‑1 (MQ=255)
                                              gTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACg                                                                                                         >  2:2803575/1‑150 (MQ=255)
                                              gTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACg                                                                                                         >  1:3301375/1‑150 (MQ=255)
                                              gTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACg                                                                                                         >  2:739858/1‑150 (MQ=255)
                                               tATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGt                                                                                                        <  1:2982758/150‑1 (MQ=255)
                                                aTCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGtt                                                                                                       >  2:3652552/1‑150 (MQ=255)
                                                 tCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTa                                                                                                      >  2:2666565/1‑150 (MQ=255)
                                                  cAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGACTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTAt                                                                                                     >  2:370333/1‑150 (MQ=255)
                                                   aaTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATc                                                                                                    <  2:1220184/150‑1 (MQ=255)
                                                    aTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCa                                                                                                   <  2:1732224/150‑1 (MQ=255)
                                                       ttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATa                                                                                                <  1:3714245/150‑1 (MQ=255)
                                                              gCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgccg                                                                                         >  2:934060/1‑150 (MQ=255)
                                                              gCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgccg                                                                                         >  1:47881/1‑150 (MQ=255)
                                                                  aaGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTaa                                                                                     >  2:1976214/1‑150 (MQ=255)
                                                                     cagcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACgg                                                                                  >  1:2773381/1‑150 (MQ=255)
                                                                     cagcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACgg                                                                                  <  1:1901062/150‑1 (MQ=255)
                                                                       gcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCa                                                                                <  2:4082202/150‑1 (MQ=255)
                                                                              gATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAAc                                                                         >  2:4524771/1‑150 (MQ=255)
                                                                                tACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCa                                                                       >  1:592931/1‑150 (MQ=255)
                                                                                 aCCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAg                                                                      <  2:3050248/150‑1 (MQ=255)
                                                                                   cgcgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGACCCAGaa                                                                    >  2:3424046/1‑150 (MQ=255)
                                                                                   cgcgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGaa                                                                    <  1:1117828/150‑1 (MQ=255)
                                                                                   cgcgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGaa                                                                    >  2:152789/1‑150 (MQ=255)
                                                                                     cgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATg                                                                  <  1:3998635/150‑1 (MQ=255)
                                                                                       cgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGcc                                                                <  2:4395114/150‑1 (MQ=255)
                                                                                       cgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGcc                                                                <  2:258352/150‑1 (MQ=255)
                                                                                       cgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGcc                                                                <  1:2943256/150‑1 (MQ=255)
                                                                                         gggTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTTTAACGATGGATACCGTCAACTGCTTAGTCCGGATGAAAACGATCCAGCAATTTCACGTTATCAATACCCTCCGGAAACGGCAGATGAACCAGAATGCCAt                                                              >  1:3745748/1‑150 (MQ=255)
                                                                                         gggTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCAt                                                              <  1:3957769/150‑1 (MQ=255)
                                                                                            tGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGa                                                           >  1:3853865/1‑150 (MQ=255)
                                                                                            tGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGa                                                           >  1:899574/1‑150 (MQ=255)
                                                                                            tGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGAAACGGCAGTTGAACCAGAATGCCAACGa                                                           >  1:3853131/1‑150 (MQ=255)
                                                                                              cAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATg                                                         <  2:1325164/150‑1 (MQ=255)
                                                                                                gggACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgt                                                       >  1:868600/1‑150 (MQ=255)
                                                                                                gggACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgt                                                       <  1:991416/150‑1 (MQ=255)
                                                                                                      cAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                                                 <  1:1014472/150‑1 (MQ=255)
                                                                                                      cAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTAGACTTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGACATCGATGGTGTTGTc                                                 >  1:4460076/1‑150 (MQ=255)
                                                                                                       aGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCg                                                >  1:1855781/1‑150 (MQ=255)
                                                                                                        gACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgg                                               >  2:389693/1‑150 (MQ=255)
                                                                                                           gcggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAt                                            <  1:471186/150‑1 (MQ=255)
                                                                                                           gcggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAt                                            <  2:2568910/150‑1 (MQ=255)
                                                                                                             ggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTc                                          >  2:597154/1‑150 (MQ=255)
                                                                                                              gcgcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                                         <  2:1340004/150‑1 (MQ=255)
                                                                                                                gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                                       <  2:2214287/150‑1 (MQ=255)
                                                                                                                   cgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCTGCATTCAGCGTa                                    <  1:129271/150‑1 (MQ=255)
                                                                                                                   cgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTa                                    <  1:129320/150‑1 (MQ=255)
                                                                                                                   cgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTa                                    <  1:2722369/150‑1 (MQ=255)
                                                                                                                    gcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTAt                                   <  2:720419/150‑1 (MQ=255)
                                                                                                                      tGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCg                                 <  2:3280525/150‑1 (MQ=255)
                                                                                                                        gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTTTTGTCGGCATTCAGCGTATCGAt                               <  1:3900776/150‑1 (MQ=255)
                                                                                                                        gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAt                               <  1:4339785/150‑1 (MQ=255)
                                                                                                                        gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAt                               <  1:2093867/150‑1 (MQ=255)
                                                                                                                        gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAt                               >  2:1638754/1‑150 (MQ=255)
                                                                                                                          acaGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa                             <  1:2530555/150‑1 (MQ=255)
                                                                                                                           caGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAg                            >  2:1451454/1‑150 (MQ=255)
                                                                                                                            aGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGc                           <  1:1778392/150‑1 (MQ=255)
                                                                                                                            aGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGCTGTTGTCGGCATTCAGCGTATCGATAAGc                           >  2:2914751/1‑150 (MQ=255)
                                                                                                                                gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTTTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAATCAGCGTATCGATAAGCTCca                       >  1:621805/1‑150 (MQ=255)
                                                                                                                                     acGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCagcagc                  >  1:1391858/1‑150 (MQ=255)
                                                                                                                                      cGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCt                 <  2:1930374/150‑1 (MQ=255)
                                                                                                                                         tGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCg              <  1:1500475/150‑1 (MQ=255)
                                                                                                                                          gTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGc             >  1:4265400/1‑150 (MQ=255)
                                                                                                                                            aaGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCtt           <  2:360302/150‑1 (MQ=255)
                                                                                                                                             aGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTc          >  2:3945784/1‑150 (MQ=255)
                                                                                                                                                 tGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtg      <  2:3224075/150‑1 (MQ=255)
                                                                                                                                                 tGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtg      <  2:2646674/150‑1 (MQ=255)
                                                                                                                                                     aaCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtggtgg  >  1:2817824/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGG  >  CP009273/552682‑552980

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: