Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,392 G→T 8.0% G213C (GGC→TGC)  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,3920GT 8.0% 915.4 / 34.5 245G213C (GGC→TGC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (101/123);  new base T (9/11);  total (110/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

ATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCA  >  CP009273/4214243‑4214541
                                                                                                                                                     |                                                                                                                                                     
 tCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCgg                                                                                                                                                      >  2:1560282/1‑150 (MQ=255)
   gagaGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCa                                                                                                                                                    <  1:2989352/150‑1 (MQ=255)
    agagTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCAc                                                                                                                                                   <  2:460398/150‑1 (MQ=255)
        tCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATc                                                                                                                                               <  2:2036663/150‑1 (MQ=255)
          caccaAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGTAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCAc                                                                                                                                             <  1:1880269/150‑1 (MQ=255)
              aaaGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGAc                                                                                                                                         <  1:2842729/150‑1 (MQ=255)
               aaGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACg                                                                                                                                        >  1:749620/1‑150 (MQ=255)
                  cgcTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCt                                                                                                                                     <  1:3874869/150‑1 (MQ=255)
                   gcTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTc                                                                                                                                    <  2:2180496/150‑1 (MQ=255)
                    cTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTcc                                                                                                                                   <  1:2985462/150‑1 (MQ=255)
                     tggtggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCg                                                                                                                                  <  1:1953912/150‑1 (MQ=255)
                       gtggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCggg                                                                                                                                <  2:2357517/150‑1 (MQ=255)
                        tggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGgc                                                                                                                               <  1:496283/150‑1 (MQ=255)
                         ggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGgcg                                                                                                                              >  2:4150239/1‑150 (MQ=255)
                         ggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGgcg                                                                                                                              >  1:1166680/1‑150 (MQ=255)
                            aGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCAc                                                                                                                           <  1:3232078/150‑1 (MQ=255)
                             ggTGGCGCTGATCTGATCCTTATTGAAACAGTTTTCGACACCCTTAACGACAATGAGGCTGTATTTGAGGTGAAAACTGAGTTTGAAGCGCTTGGAGTTGAGCTGACGATTATTATCTCCGGCACCATCACCGACGCCTCCGGGCGCACg                                                                                                                          <  2:167759/150‑1 (MQ=255)
                                ggCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGctc                                                                                                                       <  2:1060088/150‑1 (MQ=255)
                                 gcgcGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGctct                                                                                                                      <  1:3330181/150‑1 (MQ=255)
                                   gcgGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTcc                                                                                                                    <  1:1028779/150‑1 (MQ=255)
                                   gcgGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTcc                                                                                                                    <  2:1538649/150‑1 (MQ=255)
                                    cgGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCg                                                                                                                   <  1:3115247/150‑1 (MQ=255)
                                     ggATCTGATCCTGATTTAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCAACACCGACGCCTCCGGGCGCACGCTCTCCgg                                                                                                                  >  1:601376/1‑150 (MQ=255)
                                     ggATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCgg                                                                                                                  >  1:4066640/1‑150 (MQ=255)
                                      gATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCggg                                                                                                                 <  2:1793714/150‑1 (MQ=255)
                                       aTCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGc                                                                                                                <  2:2835710/150‑1 (MQ=255)
                                         cTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAg                                                                                                              <  2:4061869/150‑1 (MQ=255)
                                           gATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGac                                                                                                            >  1:4401746/1‑150 (MQ=255)
                                             tCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGacca                                                                                                          <  2:2362136/150‑1 (MQ=255)
                                             tCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGacca                                                                                                          >  1:1421973/1‑150 (MQ=255)
                                              ccTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGaccac                                                                                                         <  2:2902088/150‑1 (MQ=255)
                                                tGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCg                                                                                                       <  2:4424503/150‑1 (MQ=255)
                                                 gATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGa                                                                                                      <  1:171135/150‑1 (MQ=255)
                                                    tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                   >  2:1666699/1‑150 (MQ=255)
                                                    tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                   >  2:439623/1‑150 (MQ=255)
                                                     gAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCa                                                                                                  <  2:1366815/150‑1 (MQ=255)
                                                     gAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCa                                                                                                  <  1:4518969/150‑1 (MQ=255)
                                                      aaaCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAt                                                                                                 >  2:1912539/1‑150 (MQ=255)
                                                      aaaCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAt                                                                                                 >  1:1080934/1‑150 (MQ=255)
                                                        aCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAttt                                                                                               <  2:1152061/150‑1 (MQ=255)
                                                        aCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAttt                                                                                               <  1:4208963/150‑1 (MQ=255)
                                                         ccGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                              <  2:2451367/150‑1 (MQ=255)
                                                         ccGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                              >  2:3312726/1‑150 (MQ=255)
                                                         ccGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                              <  1:3487935/150‑1 (MQ=255)
                                                         ccGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCAtttt                                                                                              <  1:3487965/150‑1 (MQ=255)
                                                          cGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTa                                                                                             <  1:633092/150‑1 (MQ=255)
                                                           gTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTAc                                                                                            >  2:4454560/1‑150 (MQ=255)
                                                              ttCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCCTCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGACGCATTTTACAAc                                                                                         <  2:57632/150‑1 (MQ=255)
                                                                cGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTc                                                                                       >  2:4007596/1‑150 (MQ=255)
                                                                cGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTc                                                                                       <  2:2554223/150‑1 (MQ=255)
                                                                 gACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCa                                                                                      >  1:177951/1‑150 (MQ=255)
                                                                 gACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCa                                                                                      <  1:336721/150‑1 (MQ=255)
                                                                  acacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAt                                                                                     >  2:350782/1‑150 (MQ=255)
                                                                  acacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAt                                                                                     <  2:223000/150‑1 (MQ=255)
                                                                   cacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAtt                                                                                    <  2:2437582/150‑1 (MQ=255)
                                                                     cccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgc                                                                                  >  2:729606/1‑150 (MQ=255)
                                                                      ccTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcg                                                                                 <  1:1921467/150‑1 (MQ=255)
                                                                       cTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgcgc                                                                                <  1:949714/150‑1 (MQ=255)
                                                                             gCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGcc                                                                          >  1:2673017/1‑150 (MQ=255)
                                                                             gCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGcc                                                                          <  2:532188/150‑1 (MQ=255)
                                                                              ccAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCg                                                                         <  1:96914/150‑1 (MQ=255)
                                                                               cAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGa                                                                        >  2:1511224/1‑150 (MQ=255)
                                                                                aaaGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGaa                                                                       >  2:1638164/1‑150 (MQ=255)
                                                                                  aGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGc                                                                     <  1:2656503/150‑1 (MQ=255)
                                                                                        ggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGAc                                                               >  2:1487895/1‑150 (MQ=255)
                                                                                         gTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGAcc                                                              <  1:302549/150‑1 (MQ=255)
                                                                                         gTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGAcc                                                              >  2:3905871/1‑150 (MQ=255)
                                                                                           aTTTGCGGTGACAACGGAGGTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCTCCCTCTCCGACGCGTCCGTGCTCACGCTCTCGGGGCAGACCACCGAAGCATTTTACAACTCTTTGCGCCACGCGGACGCTCTGACCtt                                                            <  2:4619083/150‑1 (MQ=255)
                                                                                           aTTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCtt                                                            >  2:3060331/1‑150 (MQ=255)
                                                                                             ttGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTg                                                          <  1:1200325/150‑1 (MQ=255)
                                                                                                 ggTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCt                                                      >  1:4156380/1‑150 (MQ=255)
                                                                                                    gAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGaa                                                   <  1:2016175/150‑1 (MQ=255)
                                                                                                      aaaCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACt                                                 >  1:4391802/1‑150 (MQ=255)
                                                                                                      aaaCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACt                                                 <  2:2551083/150‑1 (MQ=255)
                                                                                                      aaaCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACt                                                 <  2:3442311/150‑1 (MQ=255)
                                                                                                          ggAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTgc                                             <  1:1813367/150‑1 (MQ=255)
                                                                                                          ggAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTgc                                             <  2:2058455/150‑1 (MQ=255)
                                                                                                            aGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTgcgc                                           >  1:3594070/1‑150 (MQ=255)
                                                                                                              tttGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTg                                         >  2:2320404/1‑150 (MQ=255)
                                                                                                                 gAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTgggg                                      >  2:2629141/1‑150 (MQ=255)
                                                                                                                  aaGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGc                                     <  1:2432889/150‑1 (MQ=255)
                                                                                                                  aaGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGAACACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGc                                     >  1:2104607/1‑150 (MQ=255)
                                                                                                                   aGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGcc                                    >  1:3773163/1‑150 (MQ=255)
                                                                                                                     cgcTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCg                                  >  1:341099/1‑150 (MQ=255)
                                                                                                                      gcTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCgc                                 >  1:289315/1‑149 (MQ=255)
                                                                                                                           gCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAAc                            <  2:1338941/150‑1 (MQ=255)
                                                                                                                            cGTTGAGCTGCCGATTATGATCTCCTGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGaaat                           <  2:3027354/150‑3 (MQ=255)
                                                                                                                            cGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACt                           <  1:1519172/150‑1 (MQ=255)
                                                                                                                             gTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTg                          >  2:3548740/1‑150 (MQ=255)
                                                                                                                             gTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTg                          <  2:1012730/150‑1 (MQ=255)
                                                                                                                              ttGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTgc                         <  2:863625/150‑1 (MQ=255)
                                                                                                                                    tGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGt                   <  2:648759/150‑1 (MQ=255)
                                                                                                                                       cGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACg                >  1:2497953/1‑150 (MQ=255)
                                                                                                                                                tctcCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGc       <  2:3302477/150‑1 (MQ=255)
                                                                                                                                                 ctcCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCt      <  2:3560104/150‑1 (MQ=255)
                                                                                                                                                  tcCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTg     >  2:3928775/1‑150 (MQ=255)
                                                                                                                                                   ccGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGt    <  1:845003/150‑1 (MQ=255)
                                                                                                                                                    cGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTc   >  2:2705827/1‑150 (MQ=255)
                                                                                                                                                    cGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTc   >  1:3429232/1‑150 (MQ=255)
                                                                                                                                                     ggCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCa  <  1:2954131/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCA  >  CP009273/4214243‑4214541

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: