Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,821,897 T→A 9.6% L173Q (CTG→CAG)  metW → methionine biosynthesis protein MetW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,821,8970TA 9.6% 196.4 / 9.3 73L173Q (CTG→CAG) metWmethionine biosynthesis protein MetW
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base A (4/3);  total (40/33)
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.93e-01

GCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGC  >  NC_002947/5821753‑5822044
                                                                                                                                                 |                                                                                                                                                    
gCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCg                                                                                                                                                  >  3:273941/1‑149 (MQ=255)
    gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGt                                                                                                                                                <  4:261824/147‑1 (MQ=255)
     tACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                             <  6:221809/149‑1 (MQ=255)
       ccTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGAc                                                                                                                                            <  3:75806/148‑1 (MQ=255)
         tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCAc                                                                                                                                        >  8:218438/1‑149 (MQ=255)
           gCGACCAACGGCCTCATGCCGGTATCGGACTTCATGCCGTATACCTTGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCTACTTCGAACAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                                      <  8:194783/149‑1 (MQ=255)
           gCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCCAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCAc                                                                                                                                        <  2:175712/147‑1 (MQ=255)
             gACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTg                                                                                                                                     >  7:268950/1‑148 (MQ=255)
                cAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCAGTTACTGCACCTTCCCCGAATTCGAAGAACTTGTGCACCAACGCCGGGCCAAGGTACAAGAACCCCTGGCCGTCGACCACATGccc                                                                                                                                  >  1:302896/1‑147 (MQ=255)
                          aTGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtg                                                                                                                       <  7:218438/149‑1 (MQ=255)
                                     ggACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGc                                                                                                            >  6:231111/1‑149 (MQ=255)
                                      gACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCt                                                                                                           <  8:242419/149‑1 (MQ=255)
                                             tGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGcc                                                                                                     >  7:261940/1‑148 (MQ=255)
                                              gCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAAGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCt                                                                                                    >  4:238465/1‑148 (MQ=255)
                                                     aCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                             <  4:28112/101‑1 (MQ=255)
                                                     aCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                             >  3:28112/1‑101 (MQ=255)
                                                       cTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATAGCCTAATCTTCTAg                                                                                          <  8:218987/149‑1 (MQ=255)
                                                          gTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtcttc                                                                                             >  1:151728/1‑143 (MQ=255)
                                                              aacaCGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                        >  1:293005/1‑144 (MQ=255)
                                                               acccGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATc                                                                                  <  8:268950/146‑1 (MQ=255)
                                                               acacGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGAACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                        >  3:98072/1‑143 (MQ=255)
                                                                   gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCa                                                                              <  5:231111/149‑1 (MQ=255)
                                                                    ccGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCAt                                                                             <  3:366401/149‑1 (MQ=255)
                                                                      gAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                         >  8:37745/1‑135 (MQ=255)
                                                                      gAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                         <  7:37745/135‑1 (MQ=255)
                                                                          aTCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCt                                                                           >  6:389094/1‑145 (MQ=255)
                                                                          aTCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCt                                                                           <  5:389094/145‑1 (MQ=255)
                                                                          aTCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACcg                                                                       >  3:179345/1‑149 (MQ=255)
                                                                          aTCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAAGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACcg                                                                       <  3:238465/149‑1 (MQ=255)
                                                                              aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa                                                                                           <  8:104620/125‑1 (MQ=255)
                                                                              aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa                                                                                           >  7:104620/1‑125 (MQ=255)
                                                                              aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                        <  3:193285/128‑1 (MQ=255)
                                                                              aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                        >  4:193285/1‑128 (MQ=255)
                                                                                ttCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcag                                                                 >  6:100846/1‑149 (MQ=255)
                                                                                     gCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTga                                                                                      <  2:232657/124‑1 (MQ=255)
                                                                                     gCACCTTCGCCGACTTCGAAGAACTTGTGCACAAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTga                                                                                      >  1:232657/1‑124 (MQ=255)
                                                                                      cACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGTCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGGGTcagca                                                               <  2:302896/146‑1 (MQ=255)
                                                                                        ccTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCg                                                          <  8:96245/149‑1 (MQ=255)
                                                                                         cTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGc                                                                                                                                    <  2:215900/74‑1 (MQ=255)
                                                                                         cTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCgg                                                         <  4:273941/149‑1 (MQ=255)
                                                                                         cTTCGCCGACTTCGAAGAACTTGTGCACAAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGc                                                                                                                                    >  1:215900/1‑74 (MQ=255)
                                                                                             gCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCt                                                     <  2:293005/149‑1 (MQ=255)
                                                                                                gACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGc                                                   >  5:203092/1‑148 (MQ=255)
                                                                                                   ttCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGAGTGGGGTGGCCGGCTATGGCCTAATCTTCTa                                                                                           <  4:95543/105‑1 (MQ=255)
                                                                                                   ttCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGAGTGGGGTGGCCGGCTATGGCCTAATCTTCTa                                                                                           >  3:95543/1‑105 (MQ=255)
                                                                                                      gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                               >  2:266518/1‑114 (MQ=255)
                                                                                                      gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                               <  1:266518/114‑1 (MQ=255)
                                                                                                      gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCa                                            >  5:1365/1‑149 (MQ=255)
                                                                                                         gaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCa                                         <  8:261940/149‑1 (MQ=255)
                                                                                                           aCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGc                                       <  3:17257/149‑1 (MQ=255)
                                                                                                               gtgCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTcgc                                   >  1:14061/1‑149 (MQ=255)
                                                                                                                  cACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGt                                >  1:200168/1‑149 (MQ=255)
                                                                                                                     gAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCt                              >  5:200958/1‑148 (MQ=255)
                                                                                                                      aaCGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGa                            >  7:227125/1‑149 (MQ=255)
                                                                                                                         gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGc                         >  8:207500/1‑149 (MQ=255)
                                                                                                                         gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGCGGCCGGCTATGGCCTAATCTTCTAGGTGAGCTCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGc                         >  5:128962/1‑149 (MQ=255)
                                                                                                                             ggCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCAc                     >  2:309245/1‑149 (MQ=255)
                                                                                                                                  aGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggag                >  2:301161/1‑149 (MQ=255)
                                                                                                                                   ggTACTCGACCGCCTGGCGGTCGACCCCTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGGCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCAGGGCCTGCAAGAGGATCAGCTCGCGGTCTTACGCCCACCggagg               >  5:164416/1‑149 (MQ=255)
                                                                                                                                   ggTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggag                >  7:25677/1‑148 (MQ=255)
                                                                                                                                     tggtCGACCGCCAGGCGGTCGAGTACCTTGGCCCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggagg               >  4:120017/4‑148 (MQ=38)
                                                                                                                                      aCTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAAt            >  4:127189/1‑149 (MQ=255)
                                                                                                                                       gtCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTAACGCGTCAGCAGTCCGGGCCTGCAAGAGCAGGAGCTAGCGGTATGACGCCCACCGGAGGGAtg           >  8:241362/2‑148 (MQ=255)
                                                                                                                                       cTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATc           <  4:179345/149‑1 (MQ=255)
                                                                                                                                        tCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTc                                     >  2:326892/1‑122 (MQ=255)
                                                                                                                                        tCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTc                                     <  1:326892/122‑1 (MQ=255)
                                                                                                                                        tCGACCGCCAGGTGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCg          >  4:138783/1‑149 (MQ=255)
                                                                                                                                          gACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCa        >  8:258170/1‑149 (MQ=255)
                                                                                                                                          gACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCa        <  2:14061/149‑1 (MQ=255)
                                                                                                                                          gACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGGGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTATGACGCCCACCGGAGGAATCGc         <  7:241362/148‑1 (MQ=255)
                                                                                                                                            ccGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCAcc      <  6:1365/149‑1 (MQ=255)
                                                                                                                                               ccTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATg   >  4:163985/1‑149 (MQ=255)
                                                                                                                                                cTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGGGGAATCGCACCATGc  <  6:203092/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGC  >  NC_002947/5821753‑5822044

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