Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,265,851 T→C 5.6% A280A (GCT→GCC prfA → peptide chain release factor RF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,265,8510TC 5.6% 220.8 / 2.8 71A280A (GCT→GCCprfApeptide chain release factor RF‑1
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/2);  total (35/36)
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.85e-01

GGGCGGGTGGTCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTG  >  NC_000913/1265703‑1265998
                                                                                                                                                    |                                                                                                                                                   
gggCGGGTGGTCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCTGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACgct                                                                                                                                                     <  7:8207/149‑1 (MQ=255)
gggCGGGTGGTCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACgct                                                                                                                                                     >  3:335973/1‑149 (MQ=255)
       tggtCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGaaa                                                                                                                                              >  8:29714/1‑149 (MQ=255)
        ggtCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAAt                                                                                                                                             >  4:185875/1‑149 (MQ=255)
          tCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAAt                                                                                                                                             >  1:165956/1‑147 (MQ=255)
          tCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAAt                                                                                                                                             <  2:165956/147‑1 (MQ=255)
          tCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATgg                                                                                                                                           >  3:9771/1‑149 (MQ=255)
           cAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACgctgct                                                                                                                                                  <  8:7432/141‑1 (MQ=255)
           cAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACgctgct                                                                                                                                                  >  7:7432/1‑141 (MQ=255)
                 gTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACg                                                                                                                                    >  1:392713/1‑149 (MQ=255)
                  ttAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGc                                                                                                                                   <  4:335973/149‑1 (MQ=255)
                     acacCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCaa                                                                                                                                >  2:377330/1‑149 (MQ=255)
                       accaccGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACa                                                                                                                              <  4:9771/149‑1 (MQ=255)
                          accGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGc                                                                                                                           <  2:400929/149‑1 (MQ=255)
                            cGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCg                                                                                                                         <  3:373689/149‑1 (MQ=255)
                            cGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCg                                                                                                                         <  7:93705/149‑1 (MQ=255)
                            cGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCg                                                                                                                         <  6:185409/149‑1 (MQ=255)
                              aTTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCg                                                                                                                         >  7:417626/1‑147 (MQ=255)
                                tCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATgg                                                                                                                                           <  4:275537/127‑1 (MQ=255)
                                tCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATgg                                                                                                                                           >  3:275537/1‑127 (MQ=255)
                                  ggCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGt                                                                                                                   <  2:22906/149‑1 (MQ=255)
                                  ggCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGt                                                                                                                   <  2:16221/149‑1 (MQ=255)
                                   gCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGc                                                                                                                                   >  7:73303/1‑132 (MQ=255)
                                   gCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGc                                                                                                                                   <  8:73303/132‑1 (MQ=255)
                                    cAATTCGTATTACTCACTTGCCGACCTGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCt                                                                                                                 <  2:373915/149‑1 (MQ=255)
                                       ttCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTAcc                                                                                                              >  5:200109/1‑149 (MQ=255)
                                       ttCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTAcc                                                                                                              >  2:313663/1‑149 (MQ=255)
                                         cGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTAccc                                                                                                             <  5:1626/148‑1 (MQ=255)
                                                     ttGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACctgct                                                                                                >  7:1628/1‑149 (MQ=255)
                                                        ccGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTggg                                                                                             <  4:77894/149‑1 (MQ=255)
                                                         cGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTggg                                                                                             >  8:98322/1‑148 (MQ=255)
                                                                           aaTGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGac                                                                          <  1:377330/149‑1 (MQ=255)
                                                                                 aGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaac                                                                    >  4:38462/1‑149 (MQ=255)
                                                                                 aGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaac                                                                    <  8:1628/149‑1 (MQ=255)
                                                                                   gACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTggg                                                                                             >  6:283371/1‑122 (MQ=255)
                                                                                   gACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTggg                                                                                             <  5:283371/122‑1 (MQ=255)
                                                                                   gACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaaccg                                                                  <  1:313663/149‑1 (MQ=255)
                                                                                     cGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaaccgta                                                                >  6:66164/1‑149 (MQ=255)
                                                                                      gAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTAc                                                               >  2:266195/1‑149 (MQ=255)
                                                                                      gAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTAc                                                               >  2:207862/1‑149 (MQ=255)
                                                                                        aCGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACt                                                              >  6:228790/1‑148 (MQ=255)
                                                                                        aCGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACt                                                              >  5:404734/1‑148 (MQ=255)
                                                                                         cGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGTGGAGTGGCGATCGCAGCGACCGTAACCGTACTTa                                                            <  1:207862/149‑1 (MQ=255)
                                                                                         cGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTa                                                            <  2:392713/149‑1 (MQ=255)
                                                                                          gTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTa                                                            >  7:50810/1‑148 (MQ=255)
                                                                                             cacaACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACaa                                                         <  3:26244/148‑1 (MQ=255)
                                                                                              acaacaTAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAAc                                                        <  6:200109/148‑1 (MQ=255)
                                                                                                aacaTAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTc                                                     >  5:194370/1‑149 (MQ=255)
                                                                                                aacaTAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTc                                                     <  3:185875/149‑1 (MQ=255)
                                                                                                aacaTAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTc                                                     >  4:344320/1‑149 (MQ=255)
                                                                                                  caTAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTccc                                                   <  8:417626/149‑1 (MQ=255)
                                                                                                     aaaaaCAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCg                                                  >  3:28192/1‑147 (MQ=255)
                                                                                                      aaaaCAAAGCTAACGCACTTTCTGTTCTCGGTGCTCGCATCGACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAg                                               <  2:86559/149‑1 (MQ=255)
                                                                                                             aGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGcgcg                                        <  7:29714/149‑1 (MQ=255)
                                                                                                                 aaaGCACTTTCTGTTCTCGGTGCTCGCATCCTCGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTAc                                    >  2:49103/1‑149 (MQ=255)
                                                                                                                  aaGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTAcc                                   >  2:39793/1‑149 (MQ=255)
                                                                                                                     cACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGAt                                <  1:266195/149‑1 (MQ=255)
                                                                                                                          tCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCAcc                            <  7:98322/148‑1 (MQ=255)
                                                                                                                          tCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCAcc                            <  4:28192/148‑1 (MQ=255)
                                                                                                                             gTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCAt                        >  3:101486/1‑149 (MQ=255)
                                                                                                                              ttCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAgggg                                            <  4:31154/128‑1 (MQ=255)
                                                                                                                              ttCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAgggg                                            >  3:31154/1‑128 (MQ=255)
                                                                                                                                 tcGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTa                                     <  5:42384/132‑1 (MQ=255)
                                                                                                                                 tcGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTa                                     >  6:42384/1‑132 (MQ=255)
                                                                                                                                  cGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACt                                                              <  6:319919/106‑1 (MQ=255)
                                                                                                                                  cGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACt                                                              >  5:319919/1‑106 (MQ=255)
                                                                                                                                   ggTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCt                  >  7:178340/1‑149 (MQ=255)
                                                                                                                                   ggTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCt                  <  6:194370/149‑1 (MQ=255)
                                                                                                                                        tCGCATCCACGTTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGc             >  5:312091/1‑149 (MQ=255)
                                                                                                                                          gCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGctc           >  1:48295/1‑149 (MQ=255)
                                                                                                                                               cACGCCGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTATCCTGCTGGGGAGTGGCGATCGCAGCGaccgtaa                                                                     <  6:251401/86‑1 (MQ=255)
                                                                                                                                               cACGCCGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaa                                                                     >  5:251401/1‑86 (MQ=255)
                                                                                                                                               cACGCCGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaa                                                                     >  5:251395/1‑86 (MQ=255)
                                                                                                                                               cACGCCGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGaccgtaa                                                                     <  6:251395/86‑1 (MQ=255)
                                                                                                                                                aCGTTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGc     <  6:312091/149‑1 (MQ=255)
                                                                                                                                                  gctgctGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCt   >  3:381613/1‑149 (MQ=255)
                                                                                                                                                  gctgctGAAATGGCAAAACGCCAACAGGCAGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCt   >  8:44174/1‑149 (MQ=255)
                                                                                                                                                  gctgctGAAANGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGCCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAAGCTGACGCTCAACCCCCt   >  8:50705/1‑149 (MQ=255)
                                                                                                                                                   ctgctgAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTg  <  5:66164/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGCGGGTGGTCAGCACGTTAACACCACCGATTCGGCAATTCGTATTACTCACTTGCCGACCGGGATTGTTGTTGAATGTCAGGACGAACGTTCACAACATAAAAACAAAGCTAAAGCACTTTCTGTTCTCGGTGCTCGCATCCACGCTGCTGAAATGGCAAAACGCCAACAGGCCGAAGCGTCTACCCGTCGTAACCTGCTGGGGAGTGGCGATCGCAGCGACCGTAACCGTACTTACAACTTCCCGCAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTG  >  NC_000913/1265703‑1265998

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