Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 598,618 G→T 11.2% E402* (GAA→TAA)  pheP → phenylalanine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273598,6180GT11.2% 862.1 / 56.0 240E402* (GAA→TAA) phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base G (95/118);  new base T (16/11);  total (111/129)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 241 total aligned reads displayed.

CTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGG  >  CP009273/598469‑598765
                                                                                                                                                     |                                                                                                                                                   
cTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTg                                                                                                                                                     >  1:330699/1‑150 (MQ=255)
 ttCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGa                                                                                                                                                    <  1:2502624/150‑1 (MQ=255)
          ggtgTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGtt                                                                                                                                           <  1:2554421/150‑1 (MQ=255)
              ttAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAg                                                                                                                                       >  2:891473/1‑150 (MQ=255)
              ttAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAg                                                                                                                                       >  1:2086197/1‑150 (MQ=255)
               tAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAgg                                                                                                                                      >  1:2564850/1‑150 (MQ=255)
               tAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACTGTTTAAgg                                                                                                                                      >  1:1450589/1‑150 (MQ=255)
               tAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAgg                                                                                                                                      >  1:3579603/1‑150 (MQ=255)
                   cAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgct                                                                                                                                  <  1:3115180/150‑1 (MQ=255)
                   cAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgct                                                                                                                                  <  2:435370/150‑1 (MQ=255)
                      cTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgctgct                                                                                                                               >  2:1293817/1‑150 (MQ=255)
                      cTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgctgct                                                                                                                               <  1:287511/150‑1 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGctct                                                                                                                             >  1:1690174/1‑150 (MQ=255)
                          ctgctgCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTAt                                                                                                                           <  1:4569124/150‑1 (MQ=255)
                           tgctgcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATc                                                                                                                          <  2:1632893/150‑1 (MQ=255)
                            gntgcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATcc                                                                                                                         >  2:3483089/1‑150 (MQ=255)
                                  gCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCgg                                                                                                                   <  1:2288249/150‑1 (MQ=255)
                                  gCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCgg                                                                                                                   <  2:1231710/150‑1 (MQ=255)
                                   cAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGc                                                                                                                  >  1:930275/1‑150 (MQ=255)
                                    aaaaaGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCa                                                                                                                 >  1:654520/1‑150 (MQ=255)
                                          cGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTAtc                                                                                                           >  1:3253843/1‑150 (MQ=255)
                                           gTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTAtct                                                                                                          >  2:2725702/1‑150 (MQ=255)
                                           gTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTAtct                                                                                                          <  1:558179/150‑1 (MQ=255)
                                           gTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTAtct                                                                                                          <  2:4511489/150‑1 (MQ=255)
                                             ttGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATctct                                                                                                        <  2:3762968/150‑1 (MQ=255)
                                               ggTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGc                                                                                                      >  2:3116054/1‑150 (MQ=255)
                                                 tCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCAt                                                                                                    <  2:682126/150‑1 (MQ=255)
                                                  ctgctgATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCAtt                                                                                                   >  1:3022958/1‑150 (MQ=255)
                                                   tgctgATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTg                                                                                                  >  1:943634/1‑150 (MQ=255)
                                                       gATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCtt                                                                                              >  1:1329385/1‑150 (MQ=255)
                                                       gATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCtt                                                                                              <  1:1333751/150‑1 (MQ=255)
                                                       gATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCtt                                                                                              >  2:4150352/1‑150 (MQ=255)
                                                         tGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCAGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTcc                                                                                            >  2:1201873/1‑150 (MQ=255)
                                                         tGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTcc                                                                                            >  1:3041147/1‑150 (MQ=255)
                                                         tGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTGTAAGTCGCTGCTCAATCCGTTCGGCAACTATCTCTGCATTGCCTTca                                                                                            >  2:700166/1‑149 (MQ=255)
                                                           gcgcTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTc                                                                                          <  2:2857573/150‑1 (MQ=255)
                                                           gcgcTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTc                                                                                          <  2:1416989/150‑1 (MQ=255)
                                                            cgcTGGTGGTAGCAACGCTGCTGTTGAACTGGGTTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCg                                                                                         <  2:3696035/150‑1 (MQ=255)
                                                             gcTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCgg                                                                                        <  1:3807893/150‑1 (MQ=255)
                                                             gcTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCgg                                                                                        >  1:1536197/1‑150 (MQ=255)
                                                              cTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGc                                                                                       >  1:1205299/1‑150 (MQ=255)
                                                              cTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGc                                                                                       >  1:2388862/1‑150 (MQ=255)
                                                                 gtggtAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGAATg                                                                                    >  2:901442/1‑150 (MQ=255)
                                                                 gtggtAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGAATg                                                                                    >  2:2999201/1‑150 (MQ=255)
                                                                   ggtAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAt                                                                                  <  1:1243566/150‑1 (MQ=255)
                                                                      aGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAtttt                                                                               <  1:3902762/150‑1 (MQ=255)
                                                                      aGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAtttt                                                                               >  2:4060097/1‑150 (MQ=255)
                                                                      aGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAtttt                                                                               <  1:1371311/150‑1 (MQ=255)
                                                                            gctgctGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTgctgct                                                                         >  2:3943020/1‑150 (MQ=255)
                                                                            gctgctGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTgctgct                                                                         >  2:1557207/1‑150 (MQ=255)
                                                                              tgctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGa                                                                       >  2:2172386/1‑150 (MQ=255)
                                                                              tgctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGa                                                                       >  1:27436/1‑150 (MQ=255)
                                                                              tgctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGa                                                                       <  1:2862010/150‑1 (MQ=255)
                                                                                ctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGAtg                                                                     <  1:2181064/150‑1 (MQ=255)
                                                                                ctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGAtg                                                                     <  2:3256225/150‑1 (MQ=255)
                                                                                  gTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGAtgtg                                                                   >  2:751086/1‑150 (MQ=255)
                                                                                    tGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCa                                                                 >  1:3165202/1‑150 (MQ=255)
                                                                                      aaCTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACg                                                               >  2:4395253/1‑150 (MQ=255)
                                                                                            aTTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgatggat                                                         >  2:557557/1‑150 (MQ=255)
                                                                                               aTGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGgatgat                                                      <  1:1193454/150‑1 (MQ=255)
                                                                                               aTGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGgatgat                                                      <  2:3933708/150‑1 (MQ=255)
                                                                                               aTGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGgatgat                                                      >  1:3989102/1‑150 (MQ=255)
                                                                                                 gATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGatat                                                    >  2:4340847/1‑150 (MQ=255)
                                                                                                   tctgtctgGCGCATCTGCGTTTTCGTGCAGCGCTGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGatattc                                                  <  1:1236557/150‑3 (MQ=255)
                                                                                                   tctgtctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATgc                                                  <  1:3555498/150‑1 (MQ=255)
                                                                                                     tgtctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATgcgc                                                >  2:4623294/1‑150 (MQ=255)
                                                                                                      gtctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCt                                               >  1:561614/1‑150 (MQ=255)
                                                                                                       tctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTGTGCTGCTGATGTGCACGATGGATGATATGCGCtt                                              >  2:2037530/1‑150 (MQ=255)
                                                                                                          ggCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTc                                           <  2:1311446/150‑1 (MQ=255)
                                                                                                          ggCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTc                                           >  1:707452/1‑150 (MQ=255)
                                                                                                           gcgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCa                                          <  1:451134/150‑1 (MQ=255)
                                                                                                            cgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAg                                         <  1:2021991/150‑1 (MQ=255)
                                                                                                             gcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGc                                        <  1:4630361/150‑1 (MQ=255)
                                                                                                              cATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCg                                       <  1:2231801/150‑1 (MQ=255)
                                                                                                                tCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGAt                                     <  1:1863453/150‑1 (MQ=255)
                                                                                                                    cGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCctg                                 <  2:2045473/150‑1 (MQ=255)
                                                                                                                      ttttCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCctgct                               <  1:170828/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  2:2498888/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACCGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  2:593664/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  2:178273/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  1:3359349/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  2:3089823/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                            <  2:2934976/150‑1 (MQ=255)
                                                                                                                          cGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCg                           >  1:3464136/1‑150 (MQ=255)
                                                                                                                          cGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCg                           <  2:1764358/150‑1 (MQ=255)
                                                                                                                             gCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGgtg                        <  2:2870204/150‑1 (MQ=255)
                                                                                                                               aGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGtgtg                      >  1:4382056/1‑150 (MQ=255)
                                                                                                                               aGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGtgtg                      >  1:2542791/1‑150 (MQ=255)
                                                                                                                                  gATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAt                   >  1:666438/1‑150 (MQ=255)
                                                                                                                                  gATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAt                   >  1:4463458/1‑150 (MQ=255)
                                                                                                                                    tGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTg                 <  1:2982301/150‑1 (MQ=255)
                                                                                                                                     gCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGt                >  1:590888/1‑150 (MQ=255)
                                                                                                                                     gCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGt                >  1:261621/1‑150 (MQ=255)
                                                                                                                                      cGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTa               >  2:3545464/1‑150 (MQ=255)
                                                                                                                                       gACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTAt              <  1:3929972/150‑1 (MQ=255)
                                                                                                                                         cGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTc            >  1:475385/1‑150 (MQ=255)
                                                                                                                                            cAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTg         >  2:3626772/1‑150 (MQ=255)
                                                                                                                                            cAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTg         <  1:1023480/150‑1 (MQ=255)
                                                                                                                                               gggCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATAATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGttt      <  2:240819/150‑1 (MQ=255)
                                                                                                                                                ggCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTa     <  1:3594511/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGG  >  CP009273/598469‑598765

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: