Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,809 G→T 100% E352* (GAA→TAA)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,214,8090GT97.6% 790.9 / 3.7 281E352* (GAA→TAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base T (129/141);  total (138/143)
Fisher's exact test for biased strand distribution p-value = 3.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 281 total aligned reads displayed.

TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATC  >  CP009273/4214664‑4214955
                                                                                                                                                 |                                                                                                                                                  
tctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGa                                                                                                                                                >  2:1185174/1‑150 (MQ=255)
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                               <  1:1505252/150‑1 (MQ=255)
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                               <  1:2216146/150‑1 (MQ=255)
  tcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATa                                                                                                                                              <  1:4814071/150‑1 (MQ=255)
    aaTATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCACACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCGGGAGCCGCTGAACAATGGCTAAGATAGc                                                                                                                                            >  1:3583713/1‑150 (MQ=255)
     atctcgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGcc                                                                                                                                           <  2:4816784/147‑1 (MQ=255)
     atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGcc                                                                                                                                           <  1:5084861/150‑1 (MQ=255)
          gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGttt                                                                                                                                      >  2:3074655/1‑150 (MQ=255)
           tcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTtg                                                                                                                                     <  1:3057834/150‑1 (MQ=255)
             ggcggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGTACATTTGCTAAGATATCCGGTTtgtg                                                                                                                                   >  1:4947225/1‑150 (MQ=255)
             ggcggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTtgtg                                                                                                                                   >  2:2171024/1‑150 (MQ=255)
                ggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGAAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCACGTAGCCTGCCGTGTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGaaa                                                                                                                                >  1:1313183/1‑149 (MQ=255)
                   tgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                             <  2:3025533/150‑1 (MQ=255)
                     ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                           >  1:5294094/1‑150 (MQ=255)
                     ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                           <  2:3819978/150‑1 (MQ=255)
                     ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                           >  2:3488221/1‑150 (MQ=255)
                      tgtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGACAGTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                                          >  2:3994354/1‑150 (MQ=255)
                       gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTg                                                                                                                         <  2:1625683/150‑1 (MQ=255)
                        tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                        >  1:5137959/1‑150 (MQ=255)
                        tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                        <  1:3171258/150‑1 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAAc                                                                                                                      <  1:396023/150‑1 (MQ=255)
                           caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACg                                                                                                                     >  2:1859138/1‑150 (MQ=255)
                            accacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                    >  2:5106933/1‑150 (MQ=255)
                                 gCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCaa                                                                                                               >  1:4581266/1‑150 (MQ=255)
                                 gCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCaa                                                                                                               >  2:3318269/1‑150 (MQ=255)
                                  ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                                                              <  1:767533/150‑1 (MQ=255)
                                  ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                                                              <  1:5206529/150‑1 (MQ=255)
                                  ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                                                              <  1:3196920/150‑1 (MQ=255)
                                  ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                                                              <  1:4128224/150‑1 (MQ=255)
                                   cacaACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                             >  2:1740550/1‑150 (MQ=255)
                                    acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGt                                                                                                            <  1:3000135/150‑1 (MQ=255)
                                    acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGt                                                                                                            >  1:2602798/1‑150 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTc                                                                                                           >  1:3097820/1‑150 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTc                                                                                                           >  1:5069537/1‑150 (MQ=255)
                                        caTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                                        >  2:1056005/1‑150 (MQ=255)
                                         atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                       <  1:372279/150‑1 (MQ=255)
                                         atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                       >  2:1964342/1‑150 (MQ=255)
                                         atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGGAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                       <  1:4285578/150‑1 (MQ=255)
                                                   gATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGTGCCGCTGAACCTTGGCTAAGATATCCTGTTTGTGAACGTGGGTGAACGCCCCAACGTCACCGGTTACGCaaa                                                                                             >  1:1144913/1‑150 (MQ=255)
                                                     tGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGt                                                                                           >  2:3631360/1‑150 (MQ=255)
                                                          cGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAg                                                                                      >  1:2170942/1‑150 (MQ=255)
                                                            tGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcg                                                                                    <  2:3850736/150‑1 (MQ=255)
                                                               aGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATATCCTGTTTTTGAACGTGGGTGAACGCACCAACGTCACCGGTTACGCTAAGTTAAAGCGCCt                                                                                 >  1:2745247/1‑150 (MQ=255)
                                                                  aGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAt                                                                              >  2:1847325/1‑150 (MQ=255)
                                                                    aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                            <  2:4455842/150‑1 (MQ=255)
                                                                    aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                            <  1:12875/150‑1 (MQ=255)
                                                                    aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAAAATTGGGAAGATAGCCGGTTTGTGTACGTGGGGGTACGCCCCAACCTCACCCGGTCCGCGAAGTTAAAGCGccgggtaa                                                                            >  1:4505031/1‑144 (MQ=255)
                                                                     aGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGTACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaa                                                                           >  2:1819648/1‑150 (MQ=255)
                                                                       gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCCACGTCACCGGTTACGCAAAGTTCAAGCGCCTGATCaaat                                                                         >  1:3297827/1‑149 (MQ=255)
                                                                        aTTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                        >  2:4006760/1‑150 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaag                                                                      <  2:136636/150‑1 (MQ=255)
                                                                           aGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGGGAACGCACCAACGTCACCGGTTACGCTAAGTTCAAGCGCCGGATAAAagaaga                                                                     >  1:60034/1‑150 (MQ=255)
                                                                            gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAagag                                                                    <  2:263424/150‑1 (MQ=255)
                                                                            gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAagag                                                                    <  2:4898330/150‑1 (MQ=255)
                                                                              gccgcGCAATCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCAGGCATGGAGCCTCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACTCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaa                                                                  <  2:4998776/150‑1 (MQ=255)
                                                                                cgcgCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                                                                <  2:864585/150‑1 (MQ=255)
                                                                                 gcgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATa                                                               <  2:347947/150‑1 (MQ=255)
                                                                                   gcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACa                                                             <  1:2092663/150‑1 (MQ=255)
                                                                                    cAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAg                                                            <  2:605990/150‑1 (MQ=255)
                                                                                     aaaCTGCCGGAAATTCCCGTAGCCTGCCGTGTGTCCGGCCTGGAGCCGCAGAACAATTGCTAACATAGCCTGTTTGTGTACGTGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCAGATCAAAGAAGAGAAATACAga                                                           >  1:2503454/1‑149 (MQ=255)
                                                                                      aaCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTACGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCATGTTTTTGAACGTGGGTGAAAGCACCAACGTCACCGGTTGCGCAAAGTTCAAGAGCCAGATCAAAGAATAGAAATaaatcg                                                          >  1:4767830/1‑145 (MQ=255)
                                                                                       aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCGAAGATAGCCGGTTTGTGAACGTGGGTGTACGCACCAACGTCACCGGTTACGCTAAGTTCAAGCGCCTGATCAAAGTAGAGAAATACAGCGa                                                         >  1:1619857/1‑150 (MQ=255)
                                                                                        cTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAg                                                        <  1:1061044/150‑1 (MQ=255)
                                                                                               aaaTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg                                                 <  1:89822/150‑1 (MQ=255)
                                                                                               aaaTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGGGAAGATAGCCGGTTTGTGAACGTGGGTGAACGCCCCCACGTCACCCGTTGCGCGAAGTTCAAGCGCCGGATCAAAGAAGAGAAATACAGCGAGGCGCTgg                                                 >  1:1946972/1‑150 (MQ=255)
                                                                                                 aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGAt                                               <  2:2924723/150‑1 (MQ=255)
                                                                                                 aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGTACAATTGGTAAGATAGCCTGTTTGTGTACGTGGGGGAACCCACCCACGTCACCCGGTCCGCGAAGTTAAAGCGCCGGATAAAAGAAGAGAAATAAAGCGAGGCGCggggt                                               >  1:7540/1‑148 (MQ=255)
                                                                                                   tCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGt                                             <  2:3912167/150‑1 (MQ=255)
                                                                                                     ccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATTTcg                                           >  2:2214331/1‑150 (MQ=255)
                                                                                                     ccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcg                                           <  2:4544840/150‑1 (MQ=255)
                                                                                                         agccagccGTGTGTCCGGCCTGGAGCCGCTGACCAGTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcgcg                                       <  2:4231809/150‑1 (MQ=255)
                                                                                                         aGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCCACGTCACCGGTTACGCTAAGTTCAAGCGCCGGATCAAATAAGAGAAATAAAGCGAGGCGCTGGATGTcgcgcg                                       >  1:725108/1‑150 (MQ=255)
                                                                                                           ccTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTc                                     >  2:2239805/1‑150 (MQ=255)
                                                                                                           ccTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCAGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTc                                     >  1:2960184/1‑150 (MQ=255)
                                                                                                             tGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCaa                                   <  2:1598860/150‑1 (MQ=255)
                                                                                                              gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAAc                                  >  1:3366176/1‑150 (MQ=255)
                                                                                                               ccGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGGTAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                                 <  1:29169/150‑1 (MQ=255)
                                                                                                               ccGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                                 <  1:4245686/150‑1 (MQ=255)
                                                                                                                    tGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTgg                            <  1:1185174/150‑1 (MQ=255)
                                                                                                                          gCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACg                      <  2:5246573/150‑1 (MQ=255)
                                                                                                                          gCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACg                      <  1:4740235/150‑1 (MQ=255)
                                                                                                                           ccTGGAGCCGCTGAACATTGGCTAATATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGATCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                     >  2:4951541/1‑150 (MQ=255)
                                                                                                                           ccTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                     >  2:4728166/1‑150 (MQ=255)
                                                                                                                            cTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgc                    <  1:136694/150‑1 (MQ=255)
                                                                                                                             tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcg                   >  2:355766/1‑150 (MQ=255)
                                                                                                                              ggAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgc                  <  2:2858231/150‑1 (MQ=255)
                                                                                                                                 gCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCa               <  2:2268430/150‑1 (MQ=255)
                                                                                                                                  ccGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAg              >  1:4673842/1‑150 (MQ=255)
                                                                                                                                   cGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGa             >  1:1337050/1‑150 (MQ=255)
                                                                                                                                    gCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAt            <  2:2769468/150‑1 (MQ=255)
                                                                                                                                    gCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAt            <  1:811740/150‑1 (MQ=255)
                                                                                                                                    gCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAt            <  1:4403447/150‑1 (MQ=255)
                                                                                                                                    gCTGAACATGGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAt            <  2:1200253/150‑1 (MQ=255)
                                                                                                                                     cTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt           <  1:1859138/150‑1 (MQ=255)
                                                                                                                                       gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTAt         >  2:208970/1‑150 (MQ=255)
                                                                                                                                         aCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGATTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCg       <  1:1056005/150‑1 (MQ=255)
                                                                                                                                         aCATTGGCTAAGATAGCCTGTGTGTGAACGTGGGTGAACGCACCATCGTCACCGGGTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATTCAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCg       <  2:4483504/150‑1 (MQ=255)
                                                                                                                                          cATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGa      >  1:2706126/1‑150 (MQ=255)
                                                                                                                                            ttGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata    >  1:651181/1‑150 (MQ=255)
                                                                                                                                              ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGGTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATc  <  1:618626/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATC  >  CP009273/4214664‑4214955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: