Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,448,094 Δ1 bp 5.4% coding (576/1500 nt) feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,448,0940T. 5.4% 237.6 / 6.3 74coding (576/1500 nt)feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (39/30);  new base . (2/2);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GACTTGTCGATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGTCACATCCTC  >  NC_000913/1447948‑1448242
                                                                                                                                                  |                                                                                                                                                    
gACTTGTCGATTCCCTTACCCCCGGGGGCGCGTTATCAGGCCGGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGtt                                                                                                                                                     <  2:236870/148‑1 (MQ=255)
  cTTGTCGATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCaa                                                                                                                                                  >  8:357235/1‑149 (MQ=255)
     gTCGATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCg                                                                                                                                               >  2:241514/1‑149 (MQ=255)
     gTCGATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGGTCAATCg                                                                                                                                               >  1:306586/1‑149 (MQ=255)
        gATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGa                                                                                                                                            >  4:220813/1‑149 (MQ=255)
        gATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGa                                                                                                                                            >  1:28174/1‑149 (MQ=255)
        gATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGa                                                                                                                                            >  1:88407/1‑149 (MQ=255)
         aTTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGAATGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGAt                                                                                                                                           >  2:227810/1‑149 (MQ=255)
             ccTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTaa                                                                                                                                        <  3:257313/148‑1 (MQ=255)
             ccTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGATGGCGTAGTGGCGGGAATTGAGCCATGGAACTATCCGTTGATGATTGGTATGTTGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAg                                                                                                                                       >  5:172435/1‑149 (MQ=255)
                  ccccAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTc                                                                                                                                  <  8:146076/149‑1 (MQ=255)
                          ggCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGGTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCa                                                                                                                          >  7:168682/1‑149 (MQ=255)
                             gcgTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTc                                                                                                                                  >  2:221248/1‑138 (MQ=255)
                             gcgTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTc                                                                                                                                  <  1:221248/138‑1 (MQ=255)
                               gTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCa                                                                                                                     <  3:7177/149‑1 (MQ=255)
                                  aTCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACt                                                                                                                   >  1:192559/1‑148 (MQ=255)
                                     aGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAACCCACGCCACTGAc                                                                                                                >  5:7151/1‑148 (MQ=255)
                                      ggCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGa                                                                                                              <  3:336436/149‑1 (MQ=255)
                                          tGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTg                                                                                                         >  7:329249/1‑149 (MQ=255)
                                          tGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTg                                                                                                         >  7:323711/1‑149 (MQ=255)
                                               gcgTAAAGAGCCGGTTGGCGTATTGGAGGGAATTTTGCCATGGAACTTTAATTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGAAGATGTTgc                                                                                                        <  2:196926/146‑1 (MQ=255)
                                               gcgTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGcgcg                                                                                                     <  6:85157/149‑1 (MQ=255)
                                               gcgTAAAGAGACTGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTCAAGATTGGTATGTGGAAGGTGATGCCAGCAATGGCAGCAGGCTGTTCAATCGTGATTAAGACTTCGGAAACCACGCCACAGACGATGTTga                                                                                                        >  1:196926/1‑145 (MQ=255)
                                                cgTAAAGAGCCGGTTTGCGTAGTGGCGGTAATTGTGCCATGGAACTTTCCGTTGATGATTGATATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGt                                                                                                    <  3:137819/149‑1 (MQ=255)
                                                 gTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGa                                                                                                              >  8:92997/1‑138 (MQ=255)
                                                 gTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGa                                                                                                              <  7:92997/138‑1 (MQ=255)
                                                 gTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTg                                                                                                   >  5:181545/1‑149 (MQ=255)
                                                  tAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTgg                                                                                                  >  5:112455/1‑149 (MQ=255)
                                                  tAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTgg                                                                                                  >  7:340748/1‑149 (MQ=255)
                                                    aaGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCg                                                                                                >  7:331097/1‑149 (MQ=255)
                                                    aaGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGAGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCg                                                                                                >  1:166420/1‑149 (MQ=255)
                                                            gTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGTGGAACTGGc                                                                                        <  2:28174/149‑1 (MQ=255)
                                                                  gTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGc                                                                                    >  2:310172/1‑147 (MQ=255)
                                                                     gTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGc                                                                               <  1:310172/149‑1 (MQ=255)
                                                                        gCGGGAATTGTTCCTTGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGAGGCGAGCACTGGCTGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAACCCACGCCACTGACGATGTTGCGCGTGGCGGAACTTGCCAGCGAGGCTg                                                                             <  1:233122/148‑1 (MQ=255)
                                                                        gCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCACTCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTg                                                                             <  7:235699/148‑1 (MQ=255)
                                                                        gCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTg                                                                             <  2:76176/148‑1 (MQ=255)
                                                                        gCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTg                                                                             <  2:88407/148‑1 (MQ=255)
                                                                          gggAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTa                                                                          <  1:298983/149‑1 (MQ=255)
                                                                                 gtgCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGa                                                                   <  8:323711/148‑1 (MQ=255)
                                                                                 gtgCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGa                                                                   <  8:329249/148‑1 (MQ=255)
                                                                                  tgCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGTTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCTCTGTCGATGTTGCGCGTGGCGGAACTGTCCAGCGTGGCTGGTATCCCTGa                                                                   <  6:172435/148‑1 (MQ=255)
                                                                                  tgCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGa                                                                   <  8:340748/148‑1 (MQ=255)
                                                                                     cATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGc                                                               >  1:42236/1‑149 (MQ=255)
                                                                                       tGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGt                                                             >  8:397377/1‑149 (MQ=255)
                                                                                         gAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGttt                                                           >  8:256802/1‑149 (MQ=255)
                                                                                          caCTTTCCGTTGATGATGGGTTTGTTGAAGGGGATGCCAGCACAGGCAGGAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATTTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGtttt                                                          <  6:242708/148‑1 (MQ=255)
                                                                                          aaCTTTCCGTTGATGATTGGTATGTGGAATGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGtttt                                                          <  3:220813/149‑1 (MQ=255)
                                                                                          aaCTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGa                                                                                  >  1:144936/1‑125 (MQ=255)
                                                                                          aaCTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGa                                                                                  <  2:144936/125‑1 (MQ=255)
                                                                                          aaCTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGtttt                                                          >  1:175235/1‑149 (MQ=255)
                                                                                             tttCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTaa                                                       >  4:66634/1‑149 (MQ=255)
                                                                                               tCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTaa                                                       <  6:112455/147‑1 (MQ=255)
                                                                                               tCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAAt                                                      <  8:331097/148‑1 (MQ=255)
                                                                                               tCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAAt                                                      <  5:189341/148‑1 (MQ=255)
                                                                                                 cGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATgtc                                                   <  1:241514/149‑1 (MQ=255)
                                                                                                 cGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATgt                                                    <  5:397368/148‑1 (MQ=255)
                                                                                                       tgatTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCAcc                                             <  5:20959/149‑1 (MQ=255)
                                                                                                       tgatTGGTATGTGGAAGGGGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCAcc                                             <  2:166420/149‑1 (MQ=255)
                                                                                                        gatTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCg                                            >  4:75486/1‑149 (MQ=255)
                                                                                                                      aGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTa                              >  5:323010/1‑149 (MQ=255)
                                                                                                                       ggTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTAt                             >  3:359682/1‑149 (MQ=255)
                                                                                                                           aTGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACGGACGATGTTGCGCGTGGCGGACCTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATgcgg                         >  3:21223/1‑149 (MQ=255)
                                                                                                                            tGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTa                              >  8:379075/1‑143 (MQ=255)
                                                                                                                            tGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTa                              <  7:379075/143‑1 (MQ=255)
                                                                                                                            tGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATgcggc                        >  3:290102/1‑149 (MQ=255)
                                                                                                                            tGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATgcggc                        >  7:360269/1‑149 (MQ=255)
                                                                                                                              ccAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCTGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATgcggcgc                      >  3:71285/1‑149 (MQ=255)
                                                                                                                               cAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATgcggcgcg                     >  5:375112/1‑149 (MQ=255)
                                                                                                                                     tGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATgtcgt                                                 >  4:55340/1‑115 (MQ=255)
                                                                                                                                     tGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCt                >  2:174136/1‑148 (MQ=255)
                                                                                                                                         agcagGCTGTTCAATCGTGATTAAGCCTTCTGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGt           <  4:71285/149‑1 (MQ=255)
                                                                                                                                          gcagGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGt           >  7:161158/1‑148 (MQ=255)
                                                                                                                                            agGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGTca         >  8:345081/1‑148 (MQ=255)
                                                                                                                                             ggCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGTcaca       >  5:122345/1‑149 (MQ=255)
                                                                                                                                                  ttCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGTCACATCCTc  <  2:42236/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GACTTGTCGATTCCCTTACCCCAGGGGGCGCGTTATCAGGCCTGGACGCGTAAAGAGCCGGTTGGCGTAGTGGCGGGAATTGTGCCATGGAACTTTCCGTTGATGATTGGTATGTGGAAGGTGATGCCAGCACTGGCAGCAGGCTGTTCAATCGTGATTAAGCCTTCGGAAACCACGCCACTGACGATGTTGCGCGTGGCGGAACTGGCCAGCGAGGCTGGTATCCCTGATGGCGTTTTTAATGTCGTCACCGGGTCAGGTGCTGTATGCGGCGCGGCCCTGACGTCACATCCTC  >  NC_000913/1447948‑1448242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: