Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 160,517 T→G 5.8% A63A (GCA→GCC) ‡ PP_0150 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947160,5170TG 5.8% 209.1 / 4.4 69A63A (GCA→GCC) ‡PP_0150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/39);  new base G (2/2);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGG  >  NC_002947/160371‑160659
                                                                                                                                                  |                                                                                                                                              
ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                              <  7:74418/149‑1 (MQ=255)
   acgcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGc                                                                                                                                           >  5:145249/1‑149 (MQ=255)
      cacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtc                                                                                                                                        >  6:328196/1‑149 (MQ=255)
      cacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtc                                                                                                                                        >  2:197219/1‑149 (MQ=255)
       acgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGTTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcg                                                                                                                                       <  1:167002/149‑1 (MQ=255)
          cgcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcg                                                                                                                                    >  5:295438/1‑149 (MQ=255)
           gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                   >  1:213552/1‑149 (MQ=255)
           gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                   >  7:319695/1‑149 (MQ=255)
           gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                   >  5:103096/1‑149 (MQ=255)
            cTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                  <  1:13556/149‑1 (MQ=255)
            cTTAACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgcc                                                                                                                                  <  4:293074/149‑3 (MQ=255)
                   gTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcg                                                                                                                           <  5:170031/149‑1 (MQ=255)
                             ggCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCAACACtt                                                                                                                 >  4:199719/1‑149 (MQ=255)
                              gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACttt                                                                                                                <  5:328196/149‑1 (MQ=255)
                              gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACttt                                                                                                                <  1:197219/149‑1 (MQ=255)
                               cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                               >  5:171728/1‑149 (MQ=255)
                                  cAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGa                                                                                                             <  2:135194/148‑1 (MQ=255)
                                   aGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGt                                                                                                           <  8:78980/149‑1 (MQ=255)
                                     ccAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcg                                                                                                         >  5:168640/1‑149 (MQ=255)
                                     ccAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcg                                                                                                         >  7:29604/1‑149 (MQ=255)
                                      cAATACCTGGTCGCTGGCGCCTCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACttt                                                                                                                <  6:207236/141‑1 (MQ=255)
                                      cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTgagctcg                                                                                                         >  5:207236/1‑142 (MQ=255)
                                      cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                        <  1:347724/149‑1 (MQ=255)
                                      cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                        <  6:336715/149‑1 (MQ=255)
                                      cAATACCTGGGCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCTGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                        <  8:191841/149‑1 (MQ=255)
                                            cTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATc                                                                                                   >  6:243480/1‑148 (MQ=255)
                                               gTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAgg                                                                                               >  8:142355/1‑149 (MQ=255)
                                                      gcgcCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTGTTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtggt                                                                                        <  7:306513/149‑1 (MQ=255)
                                                      gcgcCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtggt                                                                                        <  2:24222/149‑1 (MQ=255)
                                                        gccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGgtgt                                                                                      <  2:249534/149‑1 (MQ=255)
                                                        gccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGgtgt                                                                                      >  7:195390/1‑149 (MQ=255)
                                                         ccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTc                                                                                     <  4:26712/149‑1 (MQ=255)
                                                          cgcgcgCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCt                                                                                    <  5:180173/149‑1 (MQ=255)
                                                            ggcgcTTCCGCTACTCGCTCGCGCAGGGTTTTGAAGGCCGCCTCGATTTGTACGGCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTTGTGGTGTCTGa                                                                                  <  7:132766/148‑1 (MQ=255)
                                                             gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCCGCGTCGTCGTGGGCTTCGCGCACCACTTTGGGGTCGCGGTTCGAGGTGGTGGGGTCTGAg                                                                                 >  4:370973/1‑149 (MQ=255)
                                                             gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAg                                                                                 >  2:103715/1‑149 (MQ=255)
                                                                cTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTg                                                                              <  5:243480/149‑1 (MQ=255)
                                                                   ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGaa                                                                            <  5:13744/148‑1 (MQ=255)
                                                                   ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGaa                                                                            <  8:319695/148‑1 (MQ=255)
                                                                        aCTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGtcgat                                                                      >  8:187654/1‑149 (MQ=255)
                                                                         ctcgctcgCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGtcgatc                                                                     <  1:103715/149‑1 (MQ=255)
                                                                                gcgcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGc                                                              <  6:171728/149‑1 (MQ=255)
                                                                                gcgcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGc                                                              <  1:129315/149‑1 (MQ=255)
                                                                                  gcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGcc                                                             >  3:339759/1‑148 (MQ=255)
                                                                                    aGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCa                                                            >  2:640/1‑147 (MQ=255)
                                                                                    aGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCa                                                            >  4:87967/1‑147 (MQ=255)
                                                                                     gggTATTGAAGGCCGCCTCGAGTTGTCCGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                         <  5:294784/149‑1 (MQ=255)
                                                                                     gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                         <  8:110061/149‑1 (MQ=255)
                                                                                     gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                         <  4:139663/149‑1 (MQ=255)
                                                                                        tATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCATGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTAGATCGGACGGCCAAAGgca                                                      >  5:370955/1‑149 (MQ=255)
                                                                                          ttGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacg                                                    >  1:181158/1‑149 (MQ=255)
                                                                                              aGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                 >  4:256813/1‑148 (MQ=255)
                                                                                              aGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                 >  1:235780/1‑148 (MQ=255)
                                                                                                gccgccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgat                                              <  7:142355/149‑1 (MQ=255)
                                                                                                  cgccTGGAGTTGTACGCCGCGGATGGCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTTGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatga                                            <  3:238274/149‑1 (MQ=255)
                                                                                                   gccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatgat                                           <  3:256813/149‑1 (MQ=255)
                                                                                                    ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGatgatg                                          <  2:235780/149‑1 (MQ=255)
                                                                                                     cTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGt                                         <  8:77007/149‑1 (MQ=255)
                                                                                                           ttGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtgg                                   <  8:29604/149‑1 (MQ=255)
                                                                                                                  ccgcGGAGGTCACCGTTGCTGGCGACGAAGCGTGCAGCGTCGGCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGAGCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccac                            <  6:295438/149‑1 (MQ=255)
                                                                                                                      ggATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAgg                        <  8:195390/149‑1 (MQ=255)
                                                                                                                           tCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCGTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGGCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCt                    <  6:103096/148‑1 (MQ=255)
                                                                                                                                    cTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc          <  2:181158/149‑1 (MQ=255)
                                                                                                                                     ccacgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc          >  5:193182/4‑148 (MQ=255)
                                                                                                                                      cacgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                        <  1:316890/51‑1 (MQ=17)
                                                                                                                                      cacgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                        >  2:316890/3‑53 (MQ=17)
                                                                                                                                        cgacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCGGGTGGCCACCAGGGTCTGTGCCTGGGCCAgg      <  1:640/149‑1 (MQ=255)
                                                                                                                                           cgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCg   <  7:187654/149‑1 (MQ=255)
                                                                                                                                            gaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTAGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCgg  >  2:178169/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGG  >  NC_002947/160371‑160659

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