Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,655,518 T→C 7.8% L138P (CTC→CCC) ‡ PP_4118 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,655,5180TC 7.8% 184.5 / 5.0 65L138P (CTC→CCC) ‡PP_4118membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/3);  total (33/32)
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.17e-01

TTCGTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655370‑4655665
                                                                                                                                                     |                                                                                                                                                   
ttcgtgttcgGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCt                                                                                                                                                     >  4:58042/1‑149 (MQ=255)
   gtgttcgGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCgg                                                                                                                                                  <  4:81423/149‑1 (MQ=255)
           ccTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGc                                                                                                                                          <  2:345226/149‑1 (MQ=255)
            cTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCt                                                                                                                                         >  1:167486/1‑149 (MQ=255)
               gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        <  6:131826/147‑1 (MQ=255)
               gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGa                                                                                                                                       <  5:134962/148‑1 (MQ=255)
               gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGa                                                                                                                                       <  4:162913/148‑1 (MQ=255)
                ccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATc                                                                                                                                     >  3:20977/1‑149 (MQ=255)
                ccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGGACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATc                                                                                                                                     <  3:125226/149‑1 (MQ=255)
                     ttccttTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCgg                                                                                                                                 >  2:226865/1‑148 (MQ=255)
                        cttTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCCTTGCTGATCCCGGTAg                                                                                                                              >  7:143031/1‑148 (MQ=255)
                              gcCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTa                                                                                                                       >  8:160268/1‑149 (MQ=255)
                                ccACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACg                                                                                                                     >  3:269256/1‑149 (MQ=255)
                                  aCCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGcc                                                                                                                   >  5:142634/1‑149 (MQ=255)
                                       gCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTg                                                                                                                <  4:12490/147‑1 (MQ=255)
                                         ttACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCCGGCACTGATCATCAGCGCAATCACTGGCATCGCGCTGCTCGGCTTACTGGCCAGCAACGGCCTCGGCCCGTTGGTGATCCCGGTAGCACTCAACGCCCAGGc                                                                                                              >  1:249324/1‑148 (MQ=255)
                                            ccTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATc                                                                                                          <  3:58042/148‑1 (MQ=255)
                                                     cTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGc                                                                                                >  8:167437/1‑149 (MQ=255)
                                                          gCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCAt                                                                                                           >  1:66797/1‑133 (MQ=255)
                                                          gCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCAt                                                                                                           <  2:66797/133‑1 (MQ=255)
                                                              cATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGgcgcg                                                                                       >  5:244671/1‑149 (MQ=255)
                                                               aTTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        >  2:59605/1‑99 (MQ=38)
                                                               aTTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        <  1:59605/99‑1 (MQ=38)
                                                                    ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAATGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTgg                                                                                 >  6:65090/1‑149 (MQ=255)
                                                                    ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTAc                                                                                                                      <  2:116016/112‑1 (MQ=255)
                                                                    ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTgg                                                                                 >  2:269441/1‑149 (MQ=255)
                                                                     cTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTTGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGc                                                                                <  4:20977/149‑1 (MQ=255)
                                                                     cTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGc                                                                                <  4:269256/149‑1 (MQ=255)
                                                                      tGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGccc                                                                                                                  <  1:141965/114‑1 (MQ=255)
                                                                      tGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGccc                                                                                                                  >  2:141965/1‑114 (MQ=255)
                                                                           cccgcccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCgg                                                                          >  2:221056/1‑149 (MQ=255)
                                                                                ccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCgggg                                                                        >  1:216855/1‑146 (MQ=255)
                                                                                ccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGa                                                                      <  2:216855/148‑1 (MQ=255)
                                                                                     cc‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATGCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGAGCACTGGCCTGCGGGGGAGTGgca                                                                <  8:143031/149‑1 (MQ=255)
                                                                                            tGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTg                                                          >  3:273901/1‑149 (MQ=255)
                                                                                              aTCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTggg                                                        >  3:38549/1‑149 (MQ=255)
                                                                                                 ttCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCACGGCCCGTTGAGGATCCCGGTAACACACAACGCCCAGGCCAACAGCGCAATGAAAGGGCGCGCACGGGCAAGCGGGGGAGTGGCCGCAAAGGGTGc                                                     >  2:309104/1‑149 (MQ=255)
                                                                                                  tCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGcc                                                    >  3:65921/1‑149 (MQ=255)
                                                                                                       gcAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCt                                               <  7:167437/149‑1 (MQ=255)
                                                                                                        cAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCg                                                   >  8:304006/1‑144 (MQ=255)
                                                                                                        cAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCg                                                   <  7:304006/144‑1 (MQ=255)
                                                                                                           tCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           <  6:142634/149‑1 (MQ=255)
                                                                                                              cTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgtg                                        <  4:65921/149‑1 (MQ=255)
                                                                                                               tGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGgcag                                                               <  1:339310/125‑1 (MQ=255)
                                                                                                               tGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGgcag                                                               >  2:339310/1‑125 (MQ=255)
                                                                                                                ggCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGtt                                      <  4:273901/149‑1 (MQ=255)
                                                                                                                   aTCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTtctc                                   >  7:72269/1‑149 (MQ=255)
                                                                                                                       cgcTGCTCGGCGTACTGGCCAGCAATGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGAc                               <  5:65090/149‑1 (MQ=255)
                                                                                                                            ctCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTTGCCATCAGCGCCATGCTCTGGCGCGCACTTGCCTGCGGGTGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGc                            <  2:167486/147‑1 (MQ=255)
                                                                                                                               ggCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGAt                       <  4:38549/149‑1 (MQ=255)
                                                                                                                               ggCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGAt                       >  6:150424/1‑149 (MQ=255)
                                                                                                                                gCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATc                      <  8:72269/149‑1 (MQ=255)
                                                                                                                                 cGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTTGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     <  1:269441/149‑1 (MQ=255)
                                                                                                                                       gatcAGCAACGGGCCGAGGCCGTTGCTGGTCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTcc                                  >  7:39433/4‑130 (MQ=255)
                                                                                                                                       gatcAGCAACGGGCCGAGGCCGTTGCTGGTCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTcc                                  <  8:39433/127‑1 (MQ=255)
                                                                                                                                         ccAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           <  7:253394/119‑1 (MQ=255)
                                                                                                                                         ccAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           >  8:253394/1‑119 (MQ=255)
                                                                                                                                          cAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgt                                         >  1:335731/1‑120 (MQ=255)
                                                                                                                                          cAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgt                                         <  2:335731/120‑1 (MQ=255)
                                                                                                                                           aGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCg           >  4:304316/1‑149 (MQ=255)
                                                                                                                                           aGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCAGCGGACGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCg           >  1:110691/1‑149 (MQ=255)
                                                                                                                                             cAACGGCCTCGGCCCGTTGGTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     >  6:251445/1‑137 (MQ=255)
                                                                                                                                             cAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     <  5:251445/137‑1 (MQ=255)
                                                                                                                                             cAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGACCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCt         <  1:309104/149‑1 (MQ=255)
                                                                                                                                                   ccTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTtcat     >  6:90405/1‑147 (MQ=255)
                                                                                                                                                    cTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAg  >  4:333118/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCGTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655370‑4655665

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: ‑