Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,443,196 A→T 47.0% Q474L (CAG→CTG)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,1960AT47.0% 58.4 / inf 234Q474L (CAG→CTG) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (66/58);  new base T (57/53);  total (123/111)
Fisher's exact test for biased strand distribution p-value = 8.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345
                                                                                                                                                     |                                                                                                                                                     
aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                                                                                                                       >  2:2415787/1‑150 (MQ=255)
    gACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCa                                                                                                                                                   <  1:727030/150‑1 (MQ=255)
       acGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGa                                                                                                                                                >  2:2334672/1‑150 (MQ=255)
        cGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGaa                                                                                                                                               >  2:3066008/1‑150 (MQ=255)
        cGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGaa                                                                                                                                               >  2:3552303/1‑150 (MQ=255)
         tgAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAAc                                                                                                                                              <  1:3668621/149‑1 (MQ=255)
         ggAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAAc                                                                                                                                              >  2:209069/1‑150 (MQ=255)
          gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACg                                                                                                                                             >  2:2205995/1‑150 (MQ=255)
              atGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcc                                                                                                                                         <  2:577968/150‑1 (MQ=255)
                   ttAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAt                                                                                                                                    >  2:3082377/1‑150 (MQ=255)
                     aaCTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAttt                                                                                                                                  >  1:1216400/1‑150 (MQ=255)
                        tGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGCAATATAGCGATCGCTTACAGGCAGGGAAGTTGAGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTGTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTgg                                                                                                                               >  1:1376998/1‑150 (MQ=255)
                          ccAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcc                                                                                                                             >  1:2103814/1‑150 (MQ=255)
                             gtgtCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCg                                                                                                                          >  1:1317250/1‑150 (MQ=255)
                               gtCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                                                                        <  2:491066/150‑1 (MQ=255)
                               gtCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                                                                        >  1:3828798/1‑150 (MQ=255)
                               gtCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                                                                        >  1:2553439/1‑150 (MQ=255)
                               gtCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                                                                        >  1:2553494/1‑150 (MQ=255)
                                   ggACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggc                                                                                                                    >  1:3554236/1‑150 (MQ=255)
                                    gACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCCGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggct                                                                                                                   <  2:1887538/150‑1 (MQ=255)
                                    gACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActggct                                                                                                                   <  1:3250154/150‑1 (MQ=255)
                                     aCGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActggctg                                                                                                                  <  1:1723380/150‑1 (MQ=255)
                                     aCGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCATACAGGCAGGGACGGTGTGGGTAAACAGCCATACATTAATTGACGCTAACTTACCGTGTGGTGGGATGAAGCTGTCAGGACCGGGCCGTGATTTTGGCCCCGActggctg                                                                                                                  >  1:245374/1‑150 (MQ=255)
                                        cAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                                                               >  2:3393515/1‑150 (MQ=255)
                                        cAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                                                               <  1:1987254/150‑1 (MQ=255)
                                         aaaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                                                              >  2:2019700/1‑150 (MQ=255)
                                         aaaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                                                              <  1:685246/150‑1 (MQ=255)
                                              atctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         >  2:4081131/2‑150 (MQ=255)
                                              atctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         <  2:4142564/149‑1 (MQ=255)
                                               tctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgt                                                                                                        >  1:341824/1‑150 (MQ=255)
                                                    cAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                                                                   >  2:2155176/1‑150 (MQ=255)
                                                    cAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                                                                   <  1:1156031/150‑1 (MQ=255)
                                                     aGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaaa                                                                                                  >  1:1762607/1‑150 (MQ=255)
                                                       gCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                                                                <  1:4078460/150‑1 (MQ=255)
                                                       gCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                                                                <  2:3359365/150‑1 (MQ=255)
                                                          ctGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                                                                             <  2:887451/150‑1 (MQ=255)
                                                          ctGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                                                                             >  1:2262855/1‑150 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           <  1:3298664/150‑1 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           >  2:2897166/1‑150 (MQ=255)
                                                               atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgt                                                                                        >  2:738947/1‑150 (MQ=255)
                                                                tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                                                                       >  2:2559697/1‑150 (MQ=255)
                                                                tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                                                                       <  1:3973403/150‑1 (MQ=255)
                                                                 ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                                                                      <  1:2140435/150‑1 (MQ=255)
                                                                 ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                                                                      >  1:1330824/1‑150 (MQ=255)
                                                                  taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                                                                     >  2:791421/1‑150 (MQ=255)
                                                                      gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                                                                 >  1:538127/1‑150 (MQ=255)
                                                                         cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACgg                                                                              >  2:3012469/1‑150 (MQ=255)
                                                                           cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTa                                                                            >  1:830258/1‑150 (MQ=255)
                                                                            ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                                                                           >  2:3274994/1‑150 (MQ=255)
                                                                            ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCTACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                                                                           <  2:11133/150‑1 (MQ=255)
                                                                              aCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                                                                         <  2:950442/150‑1 (MQ=255)
                                                                                aggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAAt                                                                       >  2:2632049/1‑150 (MQ=255)
                                                                                aggcaggGACGGTGTGGGAAAACAGCCATACCTGAATTGACGCTAACTAACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATATTGGCCCCGACTGGCTGGACGGTTGGTGTGAACCTAAGTCGGTGTGTGTACGGTATTAAt                                                                       >  2:1734526/1‑150 (MQ=255)
                                                                                 ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATc                                                                      >  2:2691646/1‑150 (MQ=255)
                                                                                    aggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCTTAATCTgg                                                                   >  2:199740/1‑150 (MQ=255)
                                                                                     gggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGt                                                                  <  2:1051156/150‑1 (MQ=255)
                                                                                     gggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGt                                                                  <  1:2943591/150‑1 (MQ=255)
                                                                                     gggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGt                                                                  <  1:1570974/150‑1 (MQ=255)
                                                                                     gggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGt                                                                  <  1:3042316/150‑1 (MQ=255)
                                                                                       gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                                                                <  1:1934385/150‑1 (MQ=255)
                                                                                          ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt                                                             <  2:2056730/150‑1 (MQ=255)
                                                                                          ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt                                                             >  2:3575980/1‑150 (MQ=255)
                                                                                          ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt                                                             <  2:4048438/150‑1 (MQ=255)
                                                                                            tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                                                           <  1:1680561/150‑1 (MQ=255)
                                                                                              tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa                                                         >  2:634420/1‑150 (MQ=255)
                                                                                              tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa                                                         <  1:945553/150‑1 (MQ=255)
                                                                                                 gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt                                                      <  2:1616615/150‑1 (MQ=255)
                                                                                                  tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTa                                                     >  2:2282231/1‑150 (MQ=255)
                                                                                                  tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTa                                                     <  1:2495565/150‑1 (MQ=255)
                                                                                                   aaaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaa                                                    >  1:4102531/1‑150 (MQ=255)
                                                                                                    aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa                                                   >  2:1468090/1‑150 (MQ=255)
                                                                                                    aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa                                                   >  2:243218/1‑150 (MQ=255)
                                                                                                       aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaaa                                                <  1:1776951/150‑1 (MQ=255)
                                                                                                       aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTATGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaaa                                                >  2:3302008/1‑150 (MQ=255)
                                                                                                        gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc                                               <  1:3954762/150‑1 (MQ=255)
                                                                                                         ccATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACg                                              <  1:2418475/150‑1 (MQ=255)
                                                                                                            tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCa                                           <  2:1257204/150‑1 (MQ=255)
                                                                                                            tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCa                                           <  2:2829248/150‑1 (MQ=255)
                                                                                                            tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCa                                           >  1:123508/1‑150 (MQ=255)
                                                                                                              ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCAcc                                         <  2:3650625/150‑1 (MQ=255)
                                                                                                               cTTAATTGACGCTAACTTACCGTTTGGTGGGATGCAGCAGTCTGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                        <  1:2098639/150‑1 (MQ=255)
                                                                                                                  aaTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGt                                     >  1:1910099/1‑150 (MQ=255)
                                                                                                                        cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                               >  2:3244321/1‑150 (MQ=255)
                                                                                                                        cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                               >  1:2678954/1‑150 (MQ=255)
                                                                                                                          cTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACTGCACCTGGTGCCGtttt                             >  2:3824650/1‑150 (MQ=255)
                                                                                                                          cTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttt                             <  1:748402/150‑1 (MQ=255)
                                                                                                                          cTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttt                             >  2:3824611/1‑150 (MQ=255)
                                                                                                                               ttACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGt                        <  1:3620354/150‑1 (MQ=255)
                                                                                                                               ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGt                        <  1:4145931/150‑1 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTc                       >  1:876190/1‑150 (MQ=255)
                                                                                                                                 aCCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                      <  1:1731776/150‑1 (MQ=255)
                                                                                                                                 aCCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGGTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                      <  2:1556952/150‑1 (MQ=255)
                                                                                                                                   cGTTTGGTGGGATGAAGCTGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGGGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGa                    <  1:3664250/150‑1 (MQ=255)
                                                                                                                                   cGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGa                    <  1:561321/150‑1 (MQ=255)
                                                                                                                                     tttGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaaa                  <  1:1502673/150‑1 (MQ=255)
                                                                                                                                       tggtggGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACa                >  1:2995511/1‑150 (MQ=255)
                                                                                                                                                 aaGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc      <  1:2205995/150‑1 (MQ=255)
                                                                                                                                                 aaGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc      >  1:1101919/1‑150 (MQ=255)
                                                                                                                                                     aGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttc  >  2:3129303/1‑150 (MQ=255)
                                                                                                                                                     aGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttc  >  2:694379/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345

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: