Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,809 G→T 91.2% 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,8090GT91.2% 368.2 / 28.5 154E352* (GAA→TAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (10/4);  new base T (74/66);  total (84/70)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 154 total aligned reads displayed.

TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214664‑4214953
                                                                                                                                                 |                                                                                                                                                
tctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTaa                                                                                                                                                <  2:404561/148‑1 (MQ=255)
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAg                                                                                                                                               <  2:2383473/148‑1 (MQ=255)
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAg                                                                                                                                               <  2:2133413/148‑1 (MQ=255)
  tcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGa                                                                                                                                              >  1:1877373/1‑148 (MQ=255)
  tcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGa                                                                                                                                              <  2:1814829/148‑1 (MQ=255)
   cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                             >  2:742395/1‑148 (MQ=255)
     atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAg                                                                                                                                           >  2:2478265/1‑148 (MQ=255)
      tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATatc                                                                                                                                          >  2:179232/1‑146 (MQ=255)
      tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGc                                                                                                                                          <  1:2592201/148‑1 (MQ=255)
       atCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGcc                                                                                                                                         >  2:2751235/1‑148 (MQ=255)
         cgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTg                                                                                                                                       >  2:2440472/1‑148 (MQ=255)
         cgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTg                                                                                                                                       <  1:1398288/148‑1 (MQ=255)
          gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGt                                                                                                                                      <  2:612726/148‑1 (MQ=255)
          gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGt                                                                                                                                      <  2:612732/148‑1 (MQ=255)
            cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGttt                                                                                                                                    >  1:2088713/1‑148 (MQ=255)
              gcggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTtgt                                                                                                                                  >  2:1801458/1‑148 (MQ=255)
                ggcTGCTGTGGCACCACGCCACAACATACTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGa                                                                                                                                <  1:2322365/148‑1 (MQ=255)
                   tgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACg                                                                                                                             >  2:1771302/1‑148 (MQ=255)
                    gctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGt                                                                                                                            <  1:789316/148‑1 (MQ=255)
                       gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCGGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                         <  2:127278/148‑1 (MQ=255)
                       gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                         <  2:56286/148‑1 (MQ=255)
                       gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                         <  2:56592/148‑1 (MQ=255)
                        tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                                        <  2:1365242/148‑1 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGa                                                                                                                      >  1:1004405/1‑148 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                      <  1:1120754/148‑1 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                      <  2:2682210/148‑1 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                      <  2:2378961/148‑1 (MQ=255)
                          gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCGGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGa                                                                                                                      >  1:1004257/1‑148 (MQ=255)
                           caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGaa                                                                                                                     <  1:2853814/148‑1 (MQ=255)
                           caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGaa                                                                                                                     <  1:3128634/148‑1 (MQ=255)
                             ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACg                                                                                                                   >  1:2678910/1‑148 (MQ=255)
                              cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  <  1:3306916/148‑1 (MQ=255)
                              cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  <  1:1467304/148‑1 (MQ=255)
                               acGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCa                                                                                                                 >  2:400061/1‑148 (MQ=255)
                                cGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCAc                                                                                                                >  1:2131651/1‑148 (MQ=255)
                                   cacaACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCaa                                                                                                             <  1:2080911/148‑1 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTTCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATATCCGGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  1:1287507/1‑148 (MQ=255)
                                      aacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGt                                                                                                          <  1:2252023/148‑1 (MQ=255)
                                       acaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTc                                                                                                         >  2:2783750/1‑148 (MQ=255)
                                        caTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCa                                                                                                        >  1:134269/1‑148 (MQ=255)
                                        caTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGGTAAGATAGCCGGTTTGTGAACGTGGGTGAACGCACCAACGTCa                                                                                                        >  1:134293/1‑148 (MQ=255)
                                             tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   <  2:1395302/148‑1 (MQ=255)
                                             tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   <  2:3041474/148‑1 (MQ=255)
                                                  cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                                                              >  1:2852409/1‑148 (MQ=255)
                                                       aGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCGAAGATAGCCTGTTTGTTAACGTGGGTGAACGCACCAAACTCACCGGTTGCGCTAAGt                                                                                         >  1:2825879/1‑148 (MQ=255)
                                                       aGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGCACATTTGCGCAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGt                                                                                         >  1:2825936/1‑148 (MQ=255)
                                                          cGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCa                                                                                      <  2:1877373/148‑1 (MQ=255)
                                                             gCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgc                                                                                   >  2:2654502/1‑148 (MQ=255)
                                                             gCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgc                                                                                   >  2:2654419/1‑148 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACCTTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:946768/1‑148 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGGTAATATAGCCTGTTTGTGAACGTGGGTGAACGCACCCACCTCACCCGTTCCGCGAAGTTCAAGCGccg                                                                               >  1:878087/1‑147 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGTTAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  1:431693/1‑148 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:946758/1‑148 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               <  1:608899/148‑1 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCt                                                                               >  1:2315612/1‑148 (MQ=255)
                                                                   gAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGa                                                                             >  2:1194695/1‑148 (MQ=255)
                                                                   gAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGa                                                                             >  2:1551144/1‑148 (MQ=255)
                                                                    aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAt                                                                            <  1:1901819/148‑1 (MQ=255)
                                                                     aGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGCGAAGATATGCTGTTTTTGTACGTGGGGGTACGCACCAACCGCACCCGGTCCGCGAACTTAAAGCGCCGGATc                                                                           >  1:652541/1‑148 (MQ=255)
                                                                       gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGATCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaa                                                                         <  2:2356164/148‑1 (MQ=255)
                                                                        aTTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaaa                                                                        >  2:1530898/1‑148 (MQ=255)
                                                                         ttAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaag                                                                       <  1:3088281/148‑1 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTCAGATAGCCTGTTTGTGAACGTGGGGGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      >  2:2561348/1‑148 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGGCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTCAGTTCAAGCGCCTGATCAaaga                                                                      >  2:2561366/1‑148 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      >  2:2046851/1‑148 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      >  2:660271/1‑148 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      <  1:2199627/148‑1 (MQ=255)
                                                                            gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaag                                                                    <  1:1549113/148‑1 (MQ=255)
                                                                                 gcgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaaa                                                               <  2:3121048/148‑1 (MQ=255)
                                                                                 gcgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaaa                                                               <  2:212036/148‑1 (MQ=255)
                                                                                 gcgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaaa                                                               <  2:1480572/148‑1 (MQ=255)
                                                                                       aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGc                                                         >  1:1118156/1‑148 (MQ=255)
                                                                                       aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGc                                                         >  2:950969/1‑148 (MQ=255)
                                                                                        cTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCg                                                        >  2:1082197/1‑148 (MQ=255)
                                                                                              gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  <  2:2889851/148‑1 (MQ=255)
                                                                                                 aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg                                               >  2:2527133/1‑148 (MQ=255)
                                                                                                  ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGa                                              >  2:123719/1‑148 (MQ=255)
                                                                                                   tCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTTGCTAAGATAGCCTGTTTGTGGACGTGGGTGAACGCACCAACCTCACCGGTTCCGCTAAGTTCAAGCGCCGGATCAAAGAAGAGAAATACAGCGAGGCGCGGGAt                                             >  1:111111/1‑148 (MQ=255)
                                                                                                    cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  1:1275339/148‑1 (MQ=255)
                                                                                                       gTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcg                                         <  2:245566/148‑1 (MQ=255)
                                                                                                              gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCa                                  >  2:2630802/1‑148 (MQ=255)
                                                                                                              gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCa                                  >  2:2408909/1‑148 (MQ=255)
                                                                                                              gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCa                                  >  1:2996081/1‑148 (MQ=255)
                                                                                                                cGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAAc                                <  1:1221924/148‑1 (MQ=255)
                                                                                                                   ttGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAgg                             >  1:1186974/1‑148 (MQ=255)
                                                                                                                      tCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTgg                          >  2:1026600/1‑148 (MQ=255)
                                                                                                                       ccGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGa                         >  1:2515542/1‑148 (MQ=255)
                                                                                                                           ccTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAAc                     >  1:780642/1‑148 (MQ=255)
                                                                                                                              ggAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgc                  >  1:239012/1‑148 (MQ=255)
                                                                                                                               gAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcg                 >  1:358668/1‑148 (MQ=255)
                                                                                                                                  ccGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcgc              <  2:476731/148‑1 (MQ=255)
                                                                                                                                  ccGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcgc              >  2:3118691/1‑148 (MQ=255)
                                                                                                                                     cTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGAGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGa           >  2:477989/1‑148 (MQ=255)
                                                                                                                                     cTGAACAATGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTACAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATCTCGCGCGTCAACAGGAGGAAAACGGCGCGCAGa           >  1:3059706/1‑148 (MQ=255)
                                                                                                                                       gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         >  1:1029242/1‑148 (MQ=255)
                                                                                                                                       gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         <  2:2678910/148‑1 (MQ=255)
                                                                                                                                        aaCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTa        <  2:2030404/148‑1 (MQ=255)
                                                                                                                                         aCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTAt       <  1:742395/148‑1 (MQ=255)
                                                                                                                                              ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata  <  2:3174920/148‑1 (MQ=255)
                                                                                                                                              ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata  <  2:1535293/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214664‑4214953

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: