Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 815,853 C→T 5.7% A265A (GCG→GCA ybhK ← putative CofD superfamily transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,8530CT 5.7% 220.6 / 2.4 71A265A (GCG→GCAybhKputative CofD superfamily transferase
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (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 = 7.12e-01

TAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATC  >  NC_000913/815708‑816000
                                                                                                                                                 |                                                                                                                                                   
tAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCAccgcc                                                                                                                                                  >  6:213282/1‑149 (MQ=255)
tAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCAccgcc                                                                                                                                                  <  6:313847/149‑1 (MQ=255)
  ggAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCga                                                                                                                                                <  1:127141/149‑1 (MQ=255)
  ggAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCga                                                                                                                                                <  5:116826/149‑1 (MQ=255)
  ggAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCga                                                                                                                                                <  7:134666/149‑1 (MQ=255)
     aGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGAc                                                                                                                                             <  5:308237/149‑1 (MQ=255)
     aGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGAc                                                                                                                                             <  1:204354/149‑1 (MQ=255)
       ccGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGCCAt                                                                                                                                           >  3:300177/1‑149 (MQ=255)
       ccGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACAt                                                                                                                                           <  7:69288/149‑1 (MQ=255)
        cGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATc                                                                                                                                          >  4:125010/1‑149 (MQ=255)
        cGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATc                                                                                                                                          >  7:331551/1‑149 (MQ=255)
         gtaaTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCATGTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATcc                                                                                                                                         <  1:257160/147‑1 (MQ=255)
             ttttCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACtt                                                                                                                                     <  3:125010/149‑1 (MQ=255)
                tcgttCGGATTTCCGTTAACCTAATGCCTGTAACGCCTTTTCCAGCGCGTTATGTTACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGTATCACAATCCGCTCTTTCACCGCCGAGACATCGACtttt                                                                                                                                   <  8:252958/145‑1 (MQ=255)
                tCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCTGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACtttt                                                                                                                                   <  5:255756/148‑1 (MQ=255)
                  cTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCACGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGc                                                                                                                                >  5:170843/1‑149 (MQ=255)
                      ggATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTAGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcga                                                                                                                            <  4:21752/149‑1 (MQ=255)
                         tttCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgacga                                                                                                                         >  3:122697/1‑149 (MQ=255)
                           tCCGTTAACCAAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTAGCGGTCATGACGATACGGAATATCGCTGGCCTCCACTACCTCCAGGATCACACTCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgacgag                                                                                                                        >  1:287073/1‑147 (MQ=255)
                                   ccTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATc                                                                                                               <  4:275004/149‑1 (MQ=255)
                                   ccTAAAGCCTGTAACGCCTTATCCAGCGCGTAATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATc                                                                                                               <  8:269766/149‑1 (MQ=255)
                                    cTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCa                                                                                                              <  5:213282/149‑1 (MQ=255)
                                       aaGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCGGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATg                                                                                                           <  8:85186/149‑1 (MQ=255)
                                           ctctAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCATCCAGTACCGCCTGGCTCACAATCCGCTCTTTCACTGCCGAGACATCCACtttt                                                                                                                                   >  2:194475/4‑121 (MQ=255)
                                           cTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACTGCCGAGACATCCACtttt                                                                                                                                   >  2:199785/1‑121 (MQ=255)
                                           cTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACTGCCGAGACATCCACtttt                                                                                                                                   <  1:199785/121‑1 (MQ=255)
                                           cTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCAAGTAACTCCTGGATCACAATCCGCTCTTTCACTGCCGAGACATCCACtttt                                                                                                                                   <  1:194475/121‑1 (MQ=255)
                                             gTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTATTTCACCGCCGAGACATCCACTTTTGGCCCGACTATGACCGCATCAATGACCttt                                                                                                     >  2:353811/1‑149 (MQ=255)
                                                aCGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTa                                                                                                  >  2:140705/1‑149 (MQ=255)
                                                       ttCCAGCTCGTTATGTAACAACTGGCGGTCATGACGATACTGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACa                                                                                            <  6:170843/148‑1 (MQ=255)
                                                         ccAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATa                                                                                          <  8:331551/148‑1 (MQ=255)
                                                         ccAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATa                                                                                          <  8:122913/148‑1 (MQ=255)
                                                         ccAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATAc                                                                                         >  3:135895/1‑149 (MQ=255)
                                                           aGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGCCCGCATCAATGACCTTTTTACCAACATACTg                                                                                       >  1:22432/1‑149 (MQ=255)
                                                               cgTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTc                                                                                    >  1:33942/1‑148 (MQ=255)
                                                                 ttATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCAt                                                                                 <  3:55448/149‑1 (MQ=255)
                                                                     gTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAAtt                                                                             >  7:347925/1‑149 (MQ=255)
                                                                           acTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgc                                                                       <  6:131435/149‑1 (MQ=255)
                                                                             tGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgctt                                                                     >  2:34983/1‑149 (MQ=255)
                                                                                   tCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTc                                                               >  5:234882/1‑149 (MQ=255)
                                                                                     aTGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCCCAATCCGCTCTTTCACCGCCGCGCCATCCACTTTTGGCCCGCCGATGACCGCATCAATGACCTTTTTACCAACATACGGCTCCATAATTGCCAGCTTGCTTTc                                                               >  6:5961/1‑147 (MQ=255)
                                                                                     aTGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCCCCGCCGAGACATCCACTTTTGGCCCGACTATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCa                                                             <  2:294785/149‑1 (MQ=255)
                                                                                       gACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGCCCGCATCAATGCCCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGc                                                           >  3:246542/1‑149 (MQ=255)
                                                                                        aCGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGc                                                           <  4:135895/148‑1 (MQ=255)
                                                                                               ggAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAAt                                                    >  5:216808/1‑148 (MQ=255)
                                                                                                gAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTa                                                  <  8:96085/149‑1 (MQ=255)
                                                                                                 aaTATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAg                                                 >  4:4527/1‑149 (MQ=255)
                                                                                                  atatCGCTGTCCTCCAGTACATCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTACATAATTGCAAGCTTGCTTTCAAGATTCAAATTAGc                                                <  8:60535/149‑1 (MQ=255)
                                                                                                   tatCGCTTGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGcc                                               <  1:34983/149‑1 (MQ=255)
                                                                                                     tCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGc                                             <  8:347925/149‑1 (MQ=255)
                                                                                                       gCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCa                                            <  1:140705/148‑1 (MQ=255)
                                                                                                          ggCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTa                                        <  8:130796/149‑1 (MQ=255)
                                                                                                             cTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAaac                                     >  1:353792/1‑149 (MQ=255)
                                                                                                              tCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgacga                                                                                                                         >  4:162971/1‑64 (MQ=255)
                                                                                                              tCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgacga                                                                                                                         <  3:162971/64‑1 (MQ=255)
                                                                                                                   tACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGGAaactcacct                               >  8:343415/1‑149 (MQ=255)
                                                                                                                   tACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACAAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAaactcaact                               >  4:163946/1‑149 (MQ=255)
                                                                                                                     cctcctGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAactcaactca                             <  7:121631/149‑1 (MQ=255)
                                                                                                                     cctcctGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAactcaactca                             <  4:122697/149‑1 (MQ=255)
                                                                                                                      ctcctGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCTAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCAc                            >  4:124772/1‑149 (MQ=255)
                                                                                                                       tcctGGATCACAATCCGCTCTTTCACCGCCGAGCCATCCACTTTTGGCCCGCCGATGCCCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTGTCCAGCTGCAAATTAGCCGCAGGTAATCTCACCTCACg                           >  2:285484/1‑149 (MQ=255)
                                                                                                                          tGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGccc                        >  1:314602/1‑149 (MQ=255)
                                                                                                                           ggATCACACTCCGCTCTTTCACCGGCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGTCCTTTTTCCCAACATTCTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGGAGGTAAACTCAACTCACGcccc                       <  2:287073/149‑1 (MQ=255)
                                                                                                                           ggATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAAtt                                                   >  6:268848/1‑121 (MQ=255)
                                                                                                                           ggATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAAtt                                                   <  5:268848/121‑1 (MQ=255)
                                                                                                                            gATCACCATACGCTCTCTCAGCGCCGAGACATGCACTTTTGGCCCGTCGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGcccc                       <  8:178867/148‑1 (MQ=255)
                                                                                                                               cacaATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGAt                   >  4:149169/1‑149 (MQ=255)
                                                                                                                                  aaTCCGCTCTTTCACAGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATCACCTTTTTACCAACATACTGCTACATAATTGCCAGCTTGCTTTCCATCTTCAAATTAGCCGCAGCTAATCTCAACTCACGCCCCAGATTGc                >  5:240518/1‑149 (MQ=255)
                                                                                                                                        ctctTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCGTTCCAGCTTCAAATTAGCAGCAGGTAAACCCAACTCACGCCCCAGATTGCCGata           >  7:265338/1‑148 (MQ=255)
                                                                                                                                          ctTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGAtata         >  7:228122/1‑148 (MQ=255)
                                                                                                                                           tttCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATaaa       >  2:231963/1‑149 (MQ=255)
                                                                                                                                           tttCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATaaa       >  8:303687/1‑149 (MQ=255)
                                                                                                                                              cACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCa    >  7:42242/1‑149 (MQ=255)
                                                                                                                                                ccgccgAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCACCATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc  <  4:300177/149‑1 (MQ=255)
                                                                                                                                                ccgccgAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc  <  1:231963/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATC  >  NC_000913/815708‑816000

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