Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 373,757 G→T 5.3% S23* (TCA→TAA)  frmB ← S‑formylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273373,7570GT 5.3% 1356.3 / 30.4 341S23* (TCA→TAA) frmBS‑formylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base G (162/159);  new base T (12/8);  total (174/167)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

TGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAA  >  CP009273/373608‑373906
                                                                                                                                                     |                                                                                                                                                     
     tCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTgg                                                                                                                                                  >  1:1736092/1‑150 (MQ=255)
     tCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTgg                                                                                                                                                  >  1:100869/1‑150 (MQ=255)
       aGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGa                                                                                                                                                >  1:2739529/1‑150 (MQ=255)
       aGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGa                                                                                                                                                <  1:2684325/150‑1 (MQ=255)
       aGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGa                                                                                                                                                <  1:2305974/150‑1 (MQ=255)
           gcgTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGaataa                                                                                                                                            >  1:3288000/1‑150 (MQ=255)
           gcgTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGaataa                                                                                                                                            <  1:4573584/150‑1 (MQ=255)
               aaCGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGcc                                                                                                                                        <  2:3957483/150‑1 (MQ=255)
                aCGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCg                                                                                                                                       >  1:1310500/1‑150 (MQ=255)
                aCGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCg                                                                                                                                       <  2:623135/150‑1 (MQ=255)
                  gCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGAt                                                                                                                                     >  1:3666165/1‑150 (MQ=255)
                  gCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGAt                                                                                                                                     >  1:2517739/1‑150 (MQ=255)
                  gCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGAt                                                                                                                                     >  1:4899037/1‑150 (MQ=255)
                  gCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGAt                                                                                                                                     <  1:1967447/150‑1 (MQ=255)
                    tGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATAc                                                                                                                                   >  2:4873038/1‑150 (MQ=255)
                       aTCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACata                                                                                                                                <  1:2280867/150‑1 (MQ=255)
                       aTCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACata                                                                                                                                >  1:3932694/1‑150 (MQ=255)
                        tCCCCGATTTAGTAACGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACatat                                                                                                                               >  2:2347814/1‑150 (MQ=255)
                           ccGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATAtttt                                                                                                                            <  1:4502395/150‑1 (MQ=255)
                           ccGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATAtttt                                                                                                                            >  2:1993492/1‑150 (MQ=255)
                               tttAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCa                                                                                                                        <  1:2157118/150‑1 (MQ=255)
                                 tAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAgc                                                                                                                      >  2:2830446/1‑150 (MQ=255)
                                  aGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAgcc                                                                                                                     >  1:355889/1‑150 (MQ=255)
                                   gTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAgccg                                                                                                                    <  1:382282/150‑1 (MQ=255)
                                     aaTGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAgccgcc                                                                                                                  >  1:143144/1‑150 (MQ=255)
                                         aaaTTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCACAg                                                                                                              >  1:3232516/1‑150 (MQ=255)
                                         aaaTTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAg                                                                                                              <  2:4102201/150‑1 (MQ=255)
                                           aTTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGcg                                                                                                            >  2:3523989/1‑149 (MQ=255)
                                               tGCTCGGTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCCAAGCTGACa                                                                                                        >  2:4390289/1‑150 (MQ=255)
                                                 cTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATg                                                                                                      >  1:2850115/1‑150 (MQ=255)
                                                  tCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGt                                                                                                     <  2:1274972/150‑1 (MQ=255)
                                                    gTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGttt                                                                                                   >  2:343412/1‑150 (MQ=255)
                                                     ttGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGtttt                                                                                                  >  2:161632/1‑150 (MQ=255)
                                                      tGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGttttt                                                                                                 >  1:1060560/1‑150 (MQ=255)
                                                       gCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTc                                                                                                <  2:3743533/150‑1 (MQ=255)
                                                       gCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTc                                                                                                >  1:5063595/1‑150 (MQ=255)
                                                        cAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCa                                                                                               <  2:826787/150‑1 (MQ=255)
                                                          ggTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAAt                                                                                             <  2:1800251/150‑1 (MQ=255)
                                                          ggTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAAt                                                                                             >  1:4156335/1‑150 (MQ=255)
                                                          ggTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAAt                                                                                             <  1:1872825/150‑1 (MQ=255)
                                                           gTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATg                                                                                            <  2:3000146/150‑1 (MQ=255)
                                                            tCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGa                                                                                           >  1:991548/1‑150 (MQ=255)
                                                             cAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAg                                                                                          <  1:692627/150‑1 (MQ=255)
                                                              aGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGt                                                                                         >  2:4838196/1‑150 (MQ=255)
                                                               ggCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGtt                                                                                        <  1:834427/150‑1 (MQ=255)
                                                                gCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTc                                                                                       <  2:319036/150‑1 (MQ=255)
                                                                  cTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCa                                                                                     <  2:4674887/150‑1 (MQ=255)
                                                                      aaGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAg                                                                                 <  1:3113363/150‑1 (MQ=255)
                                                                        gCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTc                                                                               >  1:1483506/1‑150 (MQ=255)
                                                                          cAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCAt                                                                             <  1:660260/150‑1 (MQ=255)
                                                                          cAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCAt                                                                             <  2:3389434/150‑1 (MQ=255)
                                                                          cAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCAt                                                                             <  2:4142762/150‑1 (MQ=255)
                                                                              aCAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAGGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTAtc                                                                         >  2:335402/1‑150 (MQ=255)
                                                                              aCAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTAtc                                                                         <  2:2213397/150‑1 (MQ=255)
                                                                                aGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTAtctc                                                                       <  2:4315467/150‑1 (MQ=255)
                                                                                  cACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAg                                                                     <  1:139830/150‑1 (MQ=255)
                                                                                   aCCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAgg                                                                    >  2:2571179/1‑150 (MQ=255)
                                                                                    ccGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAggg                                                                   >  2:766270/1‑150 (MQ=255)
                                                                                    ccGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAggg                                                                   <  2:3845729/150‑1 (MQ=255)
                                                                                     cGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGt                                                                  >  1:3586805/1‑150 (MQ=255)
                                                                                        cAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTaa                                                               <  1:4296311/150‑1 (MQ=255)
                                                                                         aaTTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAAt                                                              <  1:712721/150‑1 (MQ=255)
                                                                                          aTTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATc                                                             >  1:52033/1‑150 (MQ=255)
                                                                                           tttttCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATcc                                                            <  1:4847698/150‑1 (MQ=255)
                                                                                               tCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCg                                                        >  2:868711/1‑150 (MQ=255)
                                                                                               tCAGTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCg                                                        <  2:2129971/150‑1 (MQ=255)
                                                                                                 aTTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGa                                                      >  2:4222503/1‑150 (MQ=255)
                                                                                                      cgGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCt                                                 >  1:2753121/1‑150 (MQ=255)
                                                                                                        gCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGc                                               <  1:1781423/150‑1 (MQ=255)
                                                                                                        gCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGc                                               >  1:4074093/1‑150 (MQ=255)
                                                                                                          ttttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGccc                                             <  2:143144/150‑1 (MQ=255)
                                                                                                          ttttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGccc                                             >  2:2451914/1‑150 (MQ=255)
                                                                                                           tttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                            <  2:922636/150‑1 (MQ=255)
                                                                                                           tttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                            <  1:1993492/150‑1 (MQ=255)
                                                                                                            ttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCgg                                           >  1:2393859/1‑150 (MQ=255)
                                                                                                              ggtggGAGAGCGACGCCGACATTCCTTTCACATTTCCGTGATTGGGACGACTGCCGATCCATATTTTGCCAGCCGCCACAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTACTCCCCCGGACGGCTGCCCgggg                                         <  2:3296997/150‑1 (MQ=255)
                                                                                                                tggGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGgag                                       <  2:2758150/150‑1 (MQ=255)
                                                                                                                  ggAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGCAATAATGCCGAGACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGagag                                     <  2:4606025/150‑1 (MQ=255)
                                                                                                                    agaTAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTa                                   >  2:1202610/1‑150 (MQ=255)
                                                                                                                    agaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTa                                   >  1:3488560/1‑150 (MQ=255)
                                                                                                                    agaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTa                                   >  1:3020880/1‑150 (MQ=255)
                                                                                                                      aTAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAtt                                 >  1:3177485/1‑150 (MQ=255)
                                                                                                                      aTAGACGCCGACATTCATTTCACATTTCAGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTAtt                                 >  1:1163613/1‑150 (MQ=255)
                                                                                                                       tAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTa                                >  1:3051004/1‑150 (MQ=255)
                                                                                                                           cGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGaa                            >  1:2033810/1‑150 (MQ=255)
                                                                                                                                aCATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAAccccc                       >  1:439596/1‑150 (MQ=255)
                                                                                                                                aCAGTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGGTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCGGCCCGGGGAGAGTATTAAGAAccccc                       <  2:1058930/150‑1 (MQ=255)
                                                                                                                                        ttcacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg               >  2:359353/1‑150 (MQ=255)
                                                                                                                                          cacatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTg             <  2:3640067/150‑1 (MQ=255)
                                                                                                                                           acatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTgg            <  1:2782044/150‑1 (MQ=255)
                                                                                                                                           acatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTgg            <  2:1362238/150‑1 (MQ=255)
                                                                                                                                           acatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTgg            <  2:2731920/150‑1 (MQ=255)
                                                                                                                                           acatttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTgg            <  2:5029896/150‑1 (MQ=255)
                                                                                                                                            catttcaGTTATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa           <  1:164089/150‑1 (MQ=255)
                                                                                                                                             atttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGAc          <  2:3380124/150‑1 (MQ=255)
                                                                                                                                              tttcaGTGATTGTGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACg         <  2:5128062/150‑1 (MQ=255)
                                                                                                                                              tttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACg         <  1:2882746/150‑1 (MQ=255)
                                                                                                                                              tttcaGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACg         <  1:4260819/150‑1 (MQ=255)
                                                                                                                                                  aGTGATTGGGAATAATGCCGAGACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGt     <  2:4998751/150‑1 (MQ=255)
                                                                                                                                                   gTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTa    >  2:1500566/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAA  >  CP009273/373608‑373906

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: ‑