Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 552,831 A→T 25.4% 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,8310AT25.4% 671.7 / 265.2 338F123I (TTC→ATC) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (124/128);  new base T (45/41);  total (169/169)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 347 total aligned reads displayed.

GCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGG  >  CP009273/552682‑552980
                                                                                                                                                     |                                                                                                                                                     
gcCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGaa                                                                                                                                                       >  1:85798/1‑150 (MQ=255)
gcCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGaa                                                                                                                                                       >  2:2668317/1‑150 (MQ=255)
 ccAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGaaa                                                                                                                                                      >  2:1256297/1‑150 (MQ=255)
  cAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAc                                                                                                                                                     >  1:5138024/1‑150 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATAcc                                                                                                                                                    >  2:508119/1‑150 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc                                                                                                                                                    >  2:4078900/1‑150 (MQ=255)
   aaCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAcc                                                                                                                                                    >  1:3916107/1‑150 (MQ=255)
      gATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTc                                                                                                                                                 >  2:3541790/1‑150 (MQ=255)
      gATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTc                                                                                                                                                 <  2:1494452/150‑1 (MQ=255)
        tatTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCa                                                                                                                                               >  2:2520625/1‑150 (MQ=255)
        tatTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCa                                                                                                                                               >  1:4511717/1‑150 (MQ=255)
         atTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCAc                                                                                                                                              <  1:2727878/150‑1 (MQ=255)
           tCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGt                                                                                                                                            >  1:4923505/1‑150 (MQ=255)
            cGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTc                                                                                                                                           <  2:2899777/150‑1 (MQ=255)
                 cgcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGt                                                                                                                                      <  2:4841303/150‑1 (MQ=255)
                    cAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCg                                                                                                                                   <  1:3681346/150‑1 (MQ=255)
                     aaTCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCgg                                                                                                                                  <  1:1494208/150‑1 (MQ=255)
                     aaTCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCgg                                                                                                                                  >  2:4594626/1‑150 (MQ=255)
                       tCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGAt                                                                                                                                <  2:1376644/150‑1 (MQ=255)
                        caccacGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATg                                                                                                                               >  1:3139360/1‑150 (MQ=255)
                            acGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATa                                                                                                                           >  1:3208286/1‑150 (MQ=255)
                              ggCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACg                                                                                                                         >  2:169910/1‑150 (MQ=255)
                                cGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGtt                                                                                                                       <  2:1812512/150‑1 (MQ=255)
                                   tGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCa                                                                                                                    >  2:3653230/1‑150 (MQ=255)
                                    gAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAg                                                                                                                   <  1:4930155/150‑1 (MQ=255)
                                         cGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACtt                                                                                                              >  2:665724/1‑150 (MQ=255)
                                         cGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACtt                                                                                                              <  2:3180165/150‑1 (MQ=255)
                                          gAAGGTATCAATGTTGTAACGCTCAGGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACttt                                                                                                             <  1:2291019/150‑1 (MQ=255)
                                          gAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACttt                                                                                                             >  1:1336319/1‑150 (MQ=255)
                                          gAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACttt                                                                                                             <  1:350065/150‑1 (MQ=255)
                                           aaGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTg                                                                                                            <  2:2212080/150‑1 (MQ=255)
                                           aaGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTg                                                                                                            >  1:3577028/1‑150 (MQ=255)
                                           aaGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTg                                                                                                            >  2:4808944/1‑150 (MQ=255)
                                             ggTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGAGGAATACGTTCCAGCACTTTGAc                                                                                                          >  2:3748606/1‑150 (MQ=255)
                                              gTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACg                                                                                                         >  1:621391/1‑150 (MQ=255)
                                                aTCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGtt                                                                                                       >  1:4024080/1‑150 (MQ=255)
                                                 tCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTa                                                                                                      <  2:3032331/150‑1 (MQ=255)
                                                     tgttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCaa                                                                                                  <  1:1402053/150‑1 (MQ=255)
                                                     tgttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCaa                                                                                                  >  1:4776897/1‑150 (MQ=255)
                                                       ttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATa                                                                                                <  2:2671626/150‑1 (MQ=255)
                                                       ttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATa                                                                                                <  1:767729/150‑1 (MQ=255)
                                                          tAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATAccc                                                                                             <  2:4134911/150‑1 (MQ=255)
                                                            aCGCTAAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgc                                                                                           >  1:1092136/1‑150 (MQ=255)
                                                               cTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCgg                                                                                        <  2:1016694/150‑1 (MQ=255)
                                                                tCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGt                                                                                       >  1:4520977/1‑150 (MQ=255)
                                                                tCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGt                                                                                       <  1:717817/150‑1 (MQ=255)
                                                                 cAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTa                                                                                      <  2:2829337/150‑1 (MQ=255)
                                                                  aaGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTaa                                                                                     <  2:3817516/150‑1 (MQ=255)
                                                                  aaGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACACGCCGGTaa                                                                                     >  1:191041/1‑150 (MQ=255)
                                                                    gcagcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACg                                                                                   >  1:3307297/1‑150 (MQ=255)
                                                                    gcagcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACg                                                                                   <  2:3408442/150‑1 (MQ=255)
                                                                      agcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGc                                                                                 >  1:4331066/1‑150 (MQ=255)
                                                                      agcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGc                                                                                 <  1:2189798/150‑1 (MQ=255)
                                                                      agcGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGc                                                                                 <  2:1943694/150‑1 (MQ=255)
                                                                        cGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAg                                                                               >  1:1651831/1‑150 (MQ=255)
                                                                          tGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGtt                                                                             >  2:257968/1‑150 (MQ=255)
                                                                          tGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGtt                                                                             >  2:2448828/1‑150 (MQ=255)
                                                                                tACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCa                                                                       <  2:4910903/150‑1 (MQ=255)
                                                                                   cgcgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGaa                                                                    <  2:4721648/150‑1 (MQ=255)
                                                                                    gcgcgGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAAt                                                                   <  1:302072/150‑1 (MQ=255)
                                                                                        ggggTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCa                                                               >  2:2272509/1‑150 (MQ=255)
                                                                                        ggggTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCa                                                               >  1:4605336/1‑150 (MQ=255)
                                                                                          ggTGCAGGGACGCAGACGCGGCGCGCGCTGGTACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATc                                                             <  1:2223170/150‑1 (MQ=255)
                                                                                          ggTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATc                                                             >  2:4029644/1‑150 (MQ=255)
                                                                                           gTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCg                                                            <  2:636821/150‑1 (MQ=255)
                                                                                              cAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATg                                                         <  2:3253174/150‑1 (MQ=255)
                                                                                               aGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATgg                                                        <  2:1980398/150‑1 (MQ=255)
                                                                                                  gACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtgt                                                     >  2:2097600/1‑150 (MQ=255)
                                                                                                   aCGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgtt                                                    <  1:42852/150‑1 (MQ=255)
                                                                                                    cGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgttg                                                   <  1:2077409/150‑1 (MQ=255)
                                                                                                     gCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgttgt                                                  >  1:3960340/1‑150 (MQ=255)
                                                                                                     gCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgttgt                                                  <  1:904429/150‑1 (MQ=255)
                                                                                                      cAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                                                 >  1:2617959/1‑150 (MQ=255)
                                                                                                          cgcgGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCa                                             >  1:906971/1‑150 (MQ=255)
                                                                                                            cggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAtt                                           >  2:2677205/1‑150 (MQ=255)
                                                                                                             ggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTc                                          <  2:774781/150‑1 (MQ=255)
                                                                                                              gcgcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                                         <  2:3460299/150‑1 (MQ=255)
                                                                                                                gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                                       <  2:2756957/150‑1 (MQ=255)
                                                                                                                gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                                       <  2:4511717/150‑1 (MQ=255)
                                                                                                                 cgcgcTGGCACAGACGACCGACGTTGTAAGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCg                                      <  1:759293/150‑1 (MQ=255)
                                                                                                                    gcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTAt                                   >  2:4997904/1‑150 (MQ=255)
                                                                                                                    gcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTAt                                   <  1:1779089/150‑1 (MQ=255)
                                                                                                                     cTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATc                                  >  2:471405/1‑150 (MQ=255)
                                                                                                                     cTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATc                                  <  1:5125593/150‑1 (MQ=255)
                                                                                                                      tGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCg                                 <  1:1231873/150‑1 (MQ=255)
                                                                                                                          acaGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa                             <  1:807663/150‑1 (MQ=255)
                                                                                                                             gacgacCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCt                          <  1:5037841/150‑1 (MQ=255)
                                                                                                                                gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCca                       <  2:4267291/150‑1 (MQ=255)
                                                                                                                                gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCca                       <  2:1158355/150‑1 (MQ=255)
                                                                                                                                gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCca                       >  1:5111380/1‑150 (MQ=255)
                                                                                                                                gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCca                       <  1:4147900/150‑1 (MQ=255)
                                                                                                                                    gacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagcag                   >  2:1321724/1‑150 (MQ=255)
                                                                                                                                      cGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCt                 <  2:1523113/150‑1 (MQ=255)
                                                                                                                                      cGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCt                 <  2:74647/150‑1 (MQ=255)
                                                                                                                                       gTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTc                >  1:3567858/1‑150 (MQ=255)
                                                                                                                                         tGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTTACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCg              <  2:1169693/150‑1 (MQ=255)
                                                                                                                                           tAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCt            <  2:4391090/150‑1 (MQ=255)
                                                                                                                                            aaGGATGGATACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCtt           >  2:1259896/1‑150 (MQ=255)
                                                                                                                                            aaGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCtt           <  2:3916107/150‑1 (MQ=255)
                                                                                                                                                    aaaCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtggtg   <  2:1623324/150‑1 (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: