Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,447,499 C→T 6.0% intergenic (‑26/+9) aroC ← / ← prmB chorismate synthase/N5‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,4990CT 6.0% 204.9 / 2.2 66intergenic (‑26/+9)aroC/prmBchorismate synthase/N5‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATG  >  NC_000913/2447352‑2447646
                                                                                                                                                   |                                                                                                                                                   
aGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCg                                                                                                                                                    >  5:111228/1‑149 (MQ=255)
  gTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtg                                                                                                                                                  >  5:81031/1‑149 (MQ=255)
     cGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGtt                                                                                                                                                <  7:136276/148‑1 (MQ=255)
                  ggAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTtat                                                                                                                                  <  6:111228/149‑1 (MQ=255)
                         cTGGCGGAAGACCATCGACGATGCAGCCGAGCGGCAGCGCGTGCGTTTCGCCGAAGGTGGTTACGCGATAGAGTTGTCCATTTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATcgc                                                                                                                           <  3:323672/149‑1 (MQ=255)
                            gCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGaa                                                                                                                        <  7:217028/149‑1 (MQ=255)
                             cGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGaaa                                                                                                                       <  6:81031/149‑1 (MQ=255)
                               gAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAAt                                                                                                                      <  6:182609/148‑1 (MQ=255)
                               gAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAAt                                                                                                                      >  5:182609/1‑148 (MQ=255)
                                  cacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtc                                                                                                                  <  1:337073/149‑1 (MQ=255)
                                  cacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtc                                                                                                                  <  2:75043/149‑1 (MQ=255)
                                  cacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtc                                                                                                                  <  6:310154/149‑1 (MQ=255)
                                     cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                                                                                                                   >  6:219158/1‑145 (MQ=255)
                                     cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                                                                                                                   <  5:219158/145‑1 (MQ=255)
                                     cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCg                                                                                                               >  7:122435/1‑149 (MQ=255)
                                     cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCg                                                                                                               >  5:344741/1‑149 (MQ=255)
                                     cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCg                                                                                                               >  8:205638/1‑149 (MQ=255)
                                          acgaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCgg                                                                                                          >  6:223994/1‑149 (MQ=255)
                                               gCAGACGAGCGCCAGCCCGTGCGATACGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGATATTGTTGTGTTTGCGTGTTTACTAAATCTTTATAAATCACGAAATGTTCTCGTGAGGCAAt                                                                                                      >  2:272912/1‑148 (MQ=255)
                                                            aGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctgc                                                                                               <  8:252693/142‑1 (MQ=255)
                                                            aGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctgc                                                                                               >  7:252693/1‑142 (MQ=255)
                                                            aGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTgg                                                                                        <  6:344741/149‑1 (MQ=255)
                                                               ccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGa                                                                                     <  4:12453/149‑1 (MQ=255)
                                                               ccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGa                                                                                     <  8:48441/149‑1 (MQ=255)
                                                                      aTTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATaaa                                                                              <  5:223994/149‑1 (MQ=255)
                                                                         cGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAaca                                                                            <  7:234395/148‑1 (MQ=255)
                                                                                  tggtTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATcgcc                                                                  >  8:73771/1‑149 (MQ=255)
                                                                                         gcgAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTa                                                           >  7:301632/1‑149 (MQ=255)
                                                                                         gcgAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTa                                                           <  5:102591/149‑1 (MQ=255)
                                                                                          cgAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAACTGTTCTCGTGCGTCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTAt                                                          <  3:229556/149‑1 (MQ=255)
                                                                                           gAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGt                                                                                                              <  6:258394/96‑1 (MQ=255)
                                                                                           gAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGt                                                                                                              >  5:258394/1‑96 (MQ=255)
                                                                                           gAAAGAGTTGTCCAATTGTGTTTCCAGCAATCAATTCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATTTTCTCGTTCGGAAATAAGCTGATCTTTTGTGAGCATAAACACACAAAAGACGCCGTTAt                                                          >  8:96525/1‑148 (MQ=255)
                                                                                            aaaGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtctc                                                                                                                >  3:305400/1‑93 (MQ=255)
                                                                                            aaaGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtctc                                                                                                                <  4:305400/93‑1 (MQ=255)
                                                                                               gagTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAAc                                                     <  7:205638/149‑1 (MQ=255)
                                                                                                 gTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACAcc                                                                        <  3:256470/128‑1 (MQ=255)
                                                                                                 gTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACAcc                                                                        >  4:256470/1‑128 (MQ=255)
                                                                                                    gTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGt                                                                                       <  8:341753/110‑1 (MQ=255)
                                                                                                    gTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGt                                                                                       >  7:341753/1‑110 (MQ=255)
                                                                                                    gTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTcca                                                 <  8:122435/148‑1 (MQ=255)
                                                                                                       cAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGTGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTccagc                                               >  6:256280/1‑147 (MQ=255)
                                                                                                       cAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGTGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTccagc                                               >  6:256281/1‑147 (MQ=255)
                                                                                                         aTTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAATTGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAgg                                           >  8:72650/1‑149 (MQ=255)
                                                                                                             tgtTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGTGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTgta                                       <  5:256281/149‑3 (MQ=255)
                                                                                                             tgtTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGTGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGaa                                       <  5:256280/149‑1 (MQ=255)
                                                                                                              gtTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGaa                                       <  7:73771/148‑1 (MQ=255)
                                                                                                                  ccAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTcca                                                 >  2:117595/1‑134 (MQ=255)
                                                                                                                  ccAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTcca                                                 <  1:117595/134‑1 (MQ=255)
                                                                                                                   cAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTta                                                                                                                                   >  6:350496/1‑51 (MQ=255)
                                                                                                                   cAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTta                                                                                                                                   <  5:350496/51‑1 (MQ=255)
                                                                                                                   cAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGaa                                  <  5:83498/148‑1 (MQ=255)
                                                                                                                        aTCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCg                            >  4:219653/1‑149 (MQ=255)
                                                                                                                           aCGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAAc                                      <  3:279348/136‑1 (MQ=255)
                                                                                                                           aCGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAAc                                      >  4:279348/1‑136 (MQ=255)
                                                                                                                                  cGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATTAATGGCGAATTGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTTAGAATAAAGACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTc                  <  1:272912/149‑1 (MQ=255)
                                                                                                                                  cGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCTGATATTGTTc                  <  2:237321/149‑1 (MQ=255)
                                                                                                                                  cGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTc                  >  6:362057/1‑149 (MQ=255)
                                                                                                                                       ttgttgTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTccagccag                                            <  6:329501/118‑1 (MQ=255)
                                                                                                                                         gttgtGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCAt               >  2:384220/1‑145 (MQ=255)
                                                                                                                                          ttgtGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTAt                                                          >  5:22582/1‑101 (MQ=255)
                                                                                                                                          ttgtGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTAt                                                          <  6:22582/101‑1 (MQ=255)
                                                                                                                                           tgtgTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATg         <  5:362057/149‑1 (MQ=255)
                                                                                                                                             tgtTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTa       <  3:219653/149‑1 (MQ=255)
                                                                                                                                                 tGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCAt   >  4:138902/1‑149 (MQ=255)
                                                                                                                                                  gCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATg  >  5:55930/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATG  >  NC_000913/2447352‑2447646

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