Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,154 A→T 9.8% L139Q (CTG→CAG)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,1540AT 9.8% 154.4 / 5.8 59L139Q (CTG→CAG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base T (3/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGCAGCAG  >  NC_002947/274007‑274301
                                                                                                                                                     |                                                                                                                                                    
tCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcagcag                                                                                                                                                     <  8:283562/149‑1 (MQ=255)
  ggCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGc                                                                                                                                                   >  2:99089/1‑149 (MQ=255)
            tAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGCTCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCt                                                                                                                                         <  5:8410/149‑1 (MQ=255)
              ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTg                                                                                                                                        <  5:42903/148‑1 (MQ=255)
              ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGt                                                                                                                                       <  2:124592/149‑1 (MQ=255)
                  cgTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGAt                                                                                                                                   <  2:133607/149‑1 (MQ=255)
                       cAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGt                                                                                                                              >  2:42965/1‑149 (MQ=255)
                         gTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTggtcgggtcg                                                                                                                            >  2:255249/1‑149 (MQ=255)
                             cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGCTCCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGtt                                                                                                                          <  1:99089/149‑1 (MQ=255)
                                cGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGa                                                                                                                     >  5:201797/1‑149 (MQ=255)
                                  ccTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGc                                                                                                                                             <  1:19187/123‑1 (MQ=255)
                                  ccTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGc                                                                                                                                             >  2:19187/1‑123 (MQ=255)
                                          tGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCt                                                                                                                                         >  6:87963/1‑119 (MQ=255)
                                          tGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCt                                                                                                                                         <  5:87963/119‑1 (MQ=255)
                                          tGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCt                                                                                                           >  5:57532/1‑149 (MQ=255)
                                           gAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTg                                                                                                          >  5:48233/1‑149 (MQ=255)
                                           gAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTg                                                                                                          >  7:69303/1‑149 (MQ=255)
                                            aaGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTgg                                                                                                         >  1:137279/1‑149 (MQ=255)
                                              ggATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCATTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGgcg                                                                                                       >  1:73601/1‑149 (MQ=255)
                                                    tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCTGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                                                                                                   >  1:269279/1‑147 (MQ=255)
                                                            ttAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCt                                                                                          <  1:42965/148‑1 (MQ=255)
                                                                 gTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCt                                                                                          <  2:53378/143‑1 (MQ=255)
                                                                 gTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCt                                                                                          >  1:53378/1‑143 (MQ=255)
                                                                      cGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCa                                                                                     <  4:34056/143‑1 (MQ=255)
                                                                      cGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCa                                                                                     >  3:34056/1‑143 (MQ=255)
                                                                      cGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtt                                                                               <  2:137279/149‑1 (MQ=255)
                                                                       gAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCttt                                                                              <  8:69303/149‑1 (MQ=255)
                                                                        aGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCtttt                                                                             <  3:225020/149‑1 (MQ=255)
                                                                             gTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTAc                                                                         <  8:14436/148‑1 (MQ=255)
                                                                               cAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACttt                                                                      >  1:128478/1‑149 (MQ=255)
                                                                                aGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACtttt                                                                     >  6:64743/1‑149 (MQ=255)
                                                                                  ccTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTc                                                                   >  7:64015/1‑149 (MQ=255)
                                                                                      ggTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGt                                                               <  5:64743/149‑1 (MQ=255)
                                                                                            cTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGc                                                         >  8:218332/1‑149 (MQ=255)
                                                                                            cTTGTAGCCCACCCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTAGGCGATGGCAACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAgga                                                         >  6:178022/1‑148 (MQ=255)
                                                                                                   cccACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgtt                                                  >  7:135223/1‑149 (MQ=255)
                                                                                                    ccACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttg                                                 >  5:300642/1‑149 (MQ=255)
                                                                                                        gacGG‑‑GCTTGAGGGTGATGAACCCTGCCTTGGCCAGCAGCAGCTGGCGTCTGGCGCTGTTTTTCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTGAGTCAGCTGGGGTATGCTTTTGTTGTTTTTCGTGTTTTGGCCAACGTTGTtg‑a                                            <  6:352/147‑1 (MQ=255)
                                                                                                        gCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtg‑a                                            <  2:73601/149‑1 (MQ=255)
                                                                                                          cGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtg‑at                                           >  5:320899/1‑148 (MQ=255)
                                                                                                               gCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑atgatg                                       >  8:287634/1‑149 (MQ=255)
                                                                                                                 ttGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCAGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGcc                                     >  4:113192/1‑149 (MQ=255)
                                                                                                                   gAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtg‑a                                            <  5:18838/140‑1 (MQ=255)
                                                                                                                   gAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtg‑a                                            >  6:18838/1‑140 (MQ=255)
                                                                                                                      ggTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCggg                                                                                                                               >  6:125551/1‑55 (MQ=255)
                                                                                                                      ggTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCggg                                 <  6:48233/148‑1 (MQ=255)
                                                                                                                       gTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGc                               >  2:192255/1‑149 (MQ=255)
                                                                                                                       gTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGc                               <  6:201797/149‑1 (MQ=255)
                                                                                                                             acagcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCG‑GGTTGATCAGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAgg                         >  5:104129/5‑148 (MQ=37)
                                                                                                                             acagcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCG‑GGTTGATCAGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGc                        <  6:104129/145‑1 (MQ=37)
                                                                                                                              aCCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAgg                         >  2:40999/1‑148 (MQ=255)
                                                                                                                               cccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTc                                                                   <  6:143142/106‑1 (MQ=255)
                                                                                                                               cccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTc                                                                   >  5:143142/1‑106 (MQ=255)
                                                                                                                               agcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttg                                                 >  7:187448/3‑124 (MQ=38)
                                                                                                                                 cGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCATTTGTACTTTTTCGAGTTCAGGCCAACGTTGGTG‑ATGATGCCCGGGCCGTAAGGcac                      >  1:283533/1‑148 (MQ=255)
                                                                                                                                   gCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAg                   >  5:328262/1‑149 (MQ=255)
                                                                                                                                    ccTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAgg                  <  6:328262/149‑1 (MQ=255)
                                                                                                                                        ggCCTGCTGCTGCTGGCCAAGGCCG‑GGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGaa             <  6:122425/149‑1 (MQ=255)
                                                                                                                                        ggCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGCATGATGCCCGGGGCGTAAGGCACAAGGtt                <  1:192255/148‑1 (MQ=255)
                                                                                                                                         gCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGaa             <  1:40999/149‑1 (MQ=255)
                                                                                                                                           cagcagcagcagGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc           <  7:287634/149‑1 (MQ=255)
                                                                                                                                                    cagGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGcagcag  >  3:287430/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGG‑‑GCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCAGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTG‑ATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGCAGCAG  >  NC_002947/274007‑274301

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