Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,318,202 T→C T72A (ACG→GCG)  trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,318,2020TC98.3% 197.0 / NA 58T72A (ACG→GCG) trpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/30);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGCCACGCCATGCTGACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGTGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGC  >  NC_000913/1318058‑1318346
                                                                                                                                                |                                                                                                                                                
ccGACGCCACGCCATGCTGACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGtt                                                                                                                                            >  2:305945/1‑151 (MQ=255)
   actccacgccaTGCTGACCGGCCCCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTTGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCg                                                                                                                                          <  2:213798/147‑1 (MQ=255)
        acgccaTGCTGACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCgg                                                                                                                                    <  1:95200/151‑1 (MQ=255)
               gCTGACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAAt                                                                                                                             <  2:443281/151‑1 (MQ=255)
                  gACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGt                                                                                                                           <  1:698255/150‑1 (MQ=255)
                     cGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTg                                                                                                                       <  2:391754/151‑1 (MQ=255)
                        cACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGc                                                                                                                     >  1:613604/1‑150 (MQ=255)
                            ggTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTg                                                                                                                >  1:803536/1‑151 (MQ=255)
                                  ggCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAgc                                                                                                          >  2:490744/1‑151 (MQ=255)
                                     gatgatTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGcgcg                                                                                                       >  2:515855/1‑151 (MQ=255)
                                        gatTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGtt                                                                                                    <  2:222624/151‑1 (MQ=255)
                                         atTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTg                                                                                                   >  2:739641/1‑151 (MQ=255)
                                           ttCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTgg                                                                                                  >  2:451228/1‑150 (MQ=255)
                                             cGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGAc                                                                                                <  2:332091/150‑1 (MQ=255)
                                             cGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACg                                                                                               <  1:413275/151‑1 (MQ=255)
                                             cGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACg                                                                                               >  1:297780/1‑151 (MQ=255)
                                             cGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTATTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACg                                                                                               >  2:827268/1‑151 (MQ=255)
                                              gggtttACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGc                                                                                              >  1:544667/4‑151 (MQ=255)
                                                 tttACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGccc                                                                                            <  2:389907/150‑1 (MQ=255)
                                                  ttACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCgg                                                                                          >  2:718745/1‑151 (MQ=255)
                                                   tACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCgg                                                                                          <  1:222615/150‑1 (MQ=255)
                                                    aCCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGc                                                                                         >  1:117538/1‑150 (MQ=255)
                                                               tCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTc                                                                               >  2:314110/1‑149 (MQ=255)
                                                                agacAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGTGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTcag                                                                             >  1:225368/4‑150 (MQ=255)
                                                                 gCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTcagca                                                                           <  2:228428/151‑1 (MQ=255)
                                                                    aGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGt                                                                        >  2:495864/1‑151 (MQ=255)
                                                                       aaCGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGt                                                                     <  2:240959/151‑1 (MQ=255)
                                                                        aCGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGtt                                                                    <  2:694961/151‑1 (MQ=255)
                                                                            cTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAAc                                                                >  2:234816/1‑151 (MQ=255)
                                                                             tGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAAc                                                                >  1:841393/1‑150 (MQ=255)
                                                                              gCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACt                                                               <  1:228052/150‑1 (MQ=255)
                                                                              gCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTg                                                              <  1:110955/151‑1 (MQ=255)
                                                                                     cACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTg                                                        >  2:435879/1‑150 (MQ=255)
                                                                                     cACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGa                                                       <  1:673931/151‑1 (MQ=255)
                                                                                          ggTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAAtt                                                   <  1:525299/150‑1 (MQ=255)
                                                                                          ggTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTc                                                  >  1:800564/1‑151 (MQ=255)
                                                                                          ggTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTc                                                  >  1:799323/1‑151 (MQ=255)
                                                                                             tAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAg                                                <  1:809214/150‑1 (MQ=255)
                                                                                              aGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAgg                                               <  2:385277/150‑1 (MQ=255)
                                                                                              aGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGa                                              <  1:728348/151‑1 (MQ=255)
                                                                                              aGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGa                                              <  2:762554/151‑1 (MQ=255)
                                                                                                ttttATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATc                                            >  2:152078/1‑151 (MQ=255)
                                                                                                 tttATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATc                                            <  1:218380/150‑1 (MQ=255)
                                                                                                    aTGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCtttt                                        <  2:123628/151‑1 (MQ=255)
                                                                                                    aTGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCtttt                                        <  1:41898/151‑1 (MQ=255)
                                                                                                     tGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCttttt                                       >  2:486606/1‑151 (MQ=255)
                                                                                                      gcgcgcCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTg                                      >  2:76759/1‑151 (MQ=255)
                                                                                                       cgcgcCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTgc                                     <  2:486596/151‑1 (MQ=255)
                                                                                                          gccgccGTGCAGCAAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTgcgc                                   >  1:66668/1‑150 (MQ=255)
                                                                                                                   cagcaAATCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTCGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACaaaa                         <  2:94037/151‑1 (MQ=255)
                                                                                                                       aaaTCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGctt                     <  2:667248/151‑1 (MQ=255)
                                                                                                                       aaaTCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGctt                     >  2:38627/1‑151 (MQ=255)
                                                                                                                       aaaTCTTCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGctt                     >  2:249375/1‑151 (MQ=255)
                                                                                                                            ttCACGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTcc                >  2:238682/1‑151 (MQ=255)
                                                                                                                                cGCTTGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCt            <  1:63272/151‑1 (MQ=255)
                                                                                                                                    tGAGATACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTgg          <  2:99570/149‑1 (MQ=255)
                                                                                                                                       gaTACAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCa      <  2:636345/150‑1 (MQ=255)
                                                                                                                                           cAGCGCGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGc  >  2:844073/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCGACGCCACGCCATGCTGACCGGCACCGGTTTCGGCGATGATTTCGGTTTTACCCATCCGCTTCGCCAGCAACGCCTGCCCCAGCACCTGGTTAGTTTTATGCGCGCCGCCGTGCAGCAAATCTTCACGCTTGAGATACAGCGTGGTGTTCGTCCCGGCTGTAATGTTCTGGCATTTGGTCAGCGCGGTTGGACGCCCGGCATAGTTTTTCAGCAGGTCGTTGAACTGAGCCTGAAATTCAGGATCTTTTTGCGCACTGACAAAAGCTTCTTCCAGCTGGCGCAGAGC  >  NC_000913/1318058‑1318346

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: