Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,008,236 G→T 13.7% P608P (CCG→CCT metE → 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,008,2360GT13.7% 1018.1 / 93.6 313P608P (CCG→CCTmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (141/129);  new base T (19/24);  total (160/153)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 314 total aligned reads displayed.

CTCTGCTGGTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTG  >  CP009273/4008087‑4008383
                                                                                                                                                     |                                                                                                                                                   
  ctgctgGTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCt                                                                                                                                                   <  1:1239704/150‑1 (MQ=255)
   tgctgGTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTg                                                                                                                                                  >  1:2922897/1‑150 (MQ=255)
    gctgGTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGc                                                                                                                                                 <  2:333035/150‑1 (MQ=255)
       ggTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGcgtc                                                                                                                                              <  2:2163600/150‑1 (MQ=255)
        gTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGcgtcg                                                                                                                                             >  1:1003650/1‑150 (MQ=255)
          cGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTa                                                                                                                                           <  1:1747389/150‑1 (MQ=255)
             tCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCg                                                                                                                                        <  2:1603361/150‑1 (MQ=255)
             tCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCg                                                                                                                                        >  2:1617200/1‑150 (MQ=255)
              cccGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGa                                                                                                                                       >  1:344192/1‑150 (MQ=255)
                   gTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTggg                                                                                                                                  >  1:3003048/1‑150 (MQ=255)
                    tGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGa                                                                                                                                 >  1:3366241/1‑150 (MQ=255)
                     gAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGAt                                                                                                                                >  1:1151120/1‑150 (MQ=255)
                      aaGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATg                                                                                                                               >  2:2566461/1‑150 (MQ=255)
                         aTGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCTTCGTAGCGACTGGGATGCGt                                                                                                                            >  2:3150092/1‑150 (MQ=255)
                          tGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTa                                                                                                                           <  2:1811566/150‑1 (MQ=255)
                          tGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTa                                                                                                                           >  1:3196042/1‑150 (MQ=255)
                            tCAGCCGTGAAACCATCGCCAAACAGAGTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTAtc                                                                                                                         <  2:1999947/150‑1 (MQ=255)
                                   tGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTg                                                                                                                  >  2:202040/1‑150 (MQ=255)
                                   tGAAACCATCGCCAAACAGATTGCGCTGGCGCTGAGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTATACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTg                                                                                                                  >  1:2246289/1‑150 (MQ=255)
                                    gAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTgg                                                                                                                 >  1:2330825/1‑150 (MQ=255)
                                    gAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTgg                                                                                                                 <  2:3202594/150‑1 (MQ=255)
                                     aaaCCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTggg                                                                                                                >  1:3590701/1‑150 (MQ=255)
                                      aaCCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTgggg                                                                                                               >  2:1069433/1‑150 (MQ=255)
                                      aaCCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTgggg                                                                                                               >  1:3876973/1‑150 (MQ=255)
                                      aaCCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTCACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTgggg                                                                                                               >  1:1615338/1‑150 (MQ=255)
                                        ccATCGCCAAACAGATTGAGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCg                                                                                                             <  2:241296/150‑1 (MQ=255)
                                          aTCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTa                                                                                                           >  1:4525/1‑150 (MQ=255)
                                             gCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTagag                                                                                                        <  2:121816/150‑1 (MQ=255)
                                                aaaCAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGcc                                                                                                     >  1:1365783/1‑150 (MQ=255)
                                                aaaCAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGcc                                                                                                     >  2:2083794/1‑150 (MQ=255)
                                                   cAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTc                                                                                                  >  1:756165/1‑150 (MQ=255)
                                                   cAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTc                                                                                                  >  1:839951/1‑150 (MQ=255)
                                                    aGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTcc                                                                                                 >  1:1881037/1‑150 (MQ=255)
                                                       ttGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTa                                                                                              <  2:2840650/150‑1 (MQ=255)
                                                        tGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTAt                                                                                             >  1:4129295/1‑150 (MQ=255)
                                                            ctggcgctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAAc                                                                                         <  2:3710084/150‑1 (MQ=255)
                                                                 gctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAAcgccgc                                                                                    <  1:1951975/150‑1 (MQ=255)
                                                                  ctgCGTGATGATGTGGCCGATCTGGAAGCAGCTGGATTTGGCATCATCCAGATTTACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCATCgccgcc                                                                                   <  2:1475908/150‑1 (MQ=255)
                                                                  ctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACgccgcc                                                                                   <  2:3685006/150‑1 (MQ=255)
                                                                  ctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACgccgcc                                                                                   <  1:3202697/150‑1 (MQ=255)
                                                                     cGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTg                                                                                <  2:720371/150‑1 (MQ=255)
                                                                      gTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTgg                                                                               >  1:1670586/1‑150 (MQ=255)
                                                                      gTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTgg                                                                               <  1:2759520/150‑1 (MQ=255)
                                                                       tgatgaAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGc                                                                              <  2:2091859/150‑1 (MQ=255)
                                                                             aGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGa                                                                        <  1:3768298/150‑1 (MQ=255)
                                                                             aGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGa                                                                        <  2:505899/150‑1 (MQ=255)
                                                                                ggCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGa                                                                     >  2:1717172/1‑150 (MQ=255)
                                                                                 gCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGac                                                                    >  1:4008284/1‑150 (MQ=255)
                                                                                 gCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGac                                                                    >  2:1185015/1‑150 (MQ=255)
                                                                                      tCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGAcacaca                                                               <  2:2186253/150‑1 (MQ=255)
                                                                                           aaGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATcc                                                          <  2:677180/150‑1 (MQ=255)
                                                                                           aaGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATcc                                                          >  2:3374411/1‑150 (MQ=255)
                                                                                               cGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCacac                                                      <  2:394464/150‑1 (MQ=255)
                                                                                                 cTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTc                                                    >  1:2532663/1‑150 (MQ=255)
                                                                                                  tGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCTTGGCGAAAGATGACACACAAATCCACACTca                                                   >  1:3674056/1‑150 (MQ=255)
                                                                                                  tGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTca                                                   >  1:1429235/1‑150 (MQ=255)
                                                                                                   ggAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcac                                                  >  2:4127603/1‑150 (MQ=255)
                                                                                                   ggAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcac                                                  >  2:2259388/1‑150 (MQ=255)
                                                                                                   ggAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcac                                                  <  2:2942319/150‑1 (MQ=255)
                                                                                                    gAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcaca                                                 <  1:573123/150‑1 (MQ=255)
                                                                                                    gAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcaca                                                 >  1:3730201/1‑150 (MQ=255)
                                                                                                     aaTTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAt                                                <  1:3004156/150‑1 (MQ=255)
                                                                                                     aaTTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAt                                                <  1:4197851/150‑1 (MQ=255)
                                                                                                      aTTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAtg                                               <  1:2604959/150‑1 (MQ=255)
                                                                                                       ttGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAtgt                                              >  2:4537237/1‑150 (MQ=255)
                                                                                                       ttGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAtgt                                              >  1:1706342/1‑150 (MQ=255)
                                                                                                          gCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGtt                                           >  1:3000499/1‑150 (MQ=255)
                                                                                                            atcatcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTAt                                         <  2:399059/150‑1 (MQ=255)
                                                                                                             tcatcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTAtt                                        <  1:3608529/150‑1 (MQ=255)
                                                                                                             tcatcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTAtt                                        >  2:1809699/1‑150 (MQ=255)
                                                                                                               atcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGc                                      <  1:3068574/150‑1 (MQ=255)
                                                                                                               atcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGc                                      >  1:2234899/1‑150 (MQ=255)
                                                                                                                    gATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGtt                                 <  1:2110858/150‑1 (MQ=255)
                                                                                                                    gATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGtt                                 >  1:1735191/1‑150 (MQ=255)
                                                                                                                     aTTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTc                                >  1:4085873/1‑150 (MQ=255)
                                                                                                                      ttGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCa                               >  2:4291480/1‑150 (MQ=255)
                                                                                                                       tGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCaa                              <  1:1406274/150‑1 (MQ=255)
                                                                                                                       tGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCaa                              <  1:1405981/150‑1 (MQ=255)
                                                                                                                        gACGAACCGGCGCTGCGCGGAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAAc                             <  2:2560137/150‑1 (MQ=255)
                                                                                                                        gACGAACCGGCGCTGCGCGAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAAc                             <  2:835792/150‑1 (MQ=255)
                                                                                                                        gACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAAc                             >  1:1992757/1‑150 (MQ=255)
                                                                                                                        gACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAAc                             >  1:926277/1‑150 (MQ=255)
                                                                                                                           gAACCGGCGCTTCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGAc                          <  2:250766/150‑1 (MQ=255)
                                                                                                                           gAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCACCGAc                          <  2:2306834/150‑1 (MQ=255)
                                                                                                                           gAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGAc                          <  2:636220/150‑1 (MQ=255)
                                                                                                                           gAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGAc                          >  2:1726431/1‑150 (MQ=255)
                                                                                                                            aaCCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGAca                         <  2:1656451/150‑1 (MQ=255)
                                                                                                                                   gcTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGa                  <  1:3886292/150‑1 (MQ=255)
                                                                                                                                     tGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACTTGTGTTATTGCGAGTTCAACGACATCATGGAtt                <  2:2991950/150‑1 (MQ=255)
                                                                                                                                       cgcgAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCg              <  2:1845677/150‑1 (MQ=255)
                                                                                                                                        gcgAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGa             >  1:3072092/1‑150 (MQ=255)
                                                                                                                                         cgAAGGTTTACCTCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGAt            >  1:3664993/1‑150 (MQ=255)
                                                                                                                                         cgAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGAt            <  2:3875194/150‑1 (MQ=255)
                                                                                                                                         cgAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGAt            >  1:1793274/1‑150 (MQ=255)
                                                                                                                                            aGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgc         >  1:680191/1‑150 (MQ=255)
                                                                                                                                             ggTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgcg        >  2:2748736/1‑150 (MQ=255)
                                                                                                                                              gTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgcgg       >  2:186837/1‑150 (MQ=255)
                                                                                                                                                ttACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgcggcg     >  2:2711067/1‑150 (MQ=255)
                                                                                                                                                ttACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgcggcg     >  1:1054164/1‑150 (MQ=255)
                                                                                                                                                  aCCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCt   >  2:3926399/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCTGCTGGTCGTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTG  >  CP009273/4008087‑4008383

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: