Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 160,518 G→C 5.7% A63G (GCA→GGA) ‡ PP_0150 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947160,5180GC 5.7% 234.0 / 3.4 70A63G (GCA→GGA) ‡PP_0150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (2/2);  total (35/35)
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.71e-01

AGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGC  >  NC_002947/160370‑160660
                                                                                                                                                      |                                                                                                                                              
aGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCTCCGTTGCTGGCGACGAAGCTTg                                                                                                                                                <  2:224406/149‑1 (MQ=255)
 ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  5:63361/149‑1 (MQ=255)
 ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  3:293033/149‑1 (MQ=255)
 ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  5:79177/149‑1 (MQ=255)
        acgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcg                                                                                                                                        >  6:333438/1‑149 (MQ=255)
            gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                    >  8:79338/1‑149 (MQ=255)
                cacaGTGCGAGTATGGCTTCAGCCAATACCTTGTCGCTGGCGGCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGACCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   <  1:253989/148‑1 (MQ=255)
                cacaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                  <  3:40061/147‑1 (MQ=255)
                    gTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcg                                                                                                                            >  5:181449/1‑149 (MQ=255)
                    gTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcg                                                                                                                            >  4:145804/1‑149 (MQ=255)
                     tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCa                                                                                                                                              >  5:57466/1‑130 (MQ=255)
                     tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCa                                                                                                                                              <  6:57466/130‑1 (MQ=255)
                     tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcgc                                                                                                                           <  5:64059/149‑1 (MQ=255)
                         aGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGcac                                                                                                                       <  6:135844/149‑1 (MQ=255)
                               gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGAGCTTCGCGCACCACttt                                                                                                                 >  6:201238/1‑149 (MQ=255)
                                cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGGTGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                >  7:146187/1‑149 (MQ=255)
                                 ttCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                >  3:289719/1‑148 (MQ=255)
                                       cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                         <  2:256316/149‑1 (MQ=255)
                                       cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                         <  6:181449/149‑1 (MQ=255)
                                        aaTACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGGTGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgcg                                                                                                        <  8:146187/149‑1 (MQ=255)
                                        aaTACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgcg                                                                                                        <  3:145804/149‑1 (MQ=255)
                                              tGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATc                                                                                                    >  6:77954/1‑147 (MQ=255)
                                                gTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAg                                                                                                 <  5:333438/148‑1 (MQ=255)
                                                     tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                           >  6:161854/1‑149 (MQ=255)
                                                          ccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTc                                                                                      <  5:161854/149‑1 (MQ=255)
                                                          ccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTc                                                                                      <  6:108197/149‑1 (MQ=255)
                                                           cgcgcgCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                    <  1:16269/102‑2 (MQ=255)
                                                           cgcgcgCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                    >  2:16269/1‑101 (MQ=255)
                                                            gcgcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCt                                                                                     >  5:165062/1‑148 (MQ=255)
                                                              gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAg                                                                                  <  7:25765/149‑1 (MQ=255)
                                                                gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtc                                                                                                                                      >  5:19278/1‑95 (MQ=255)
                                                                gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtc                                                                                                                                      <  6:19278/95‑1 (MQ=255)
                                                                gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATAGTGGTGGTGGTGTCTGAGGt                                                                                >  8:256140/1‑149 (MQ=255)
                                                                 cTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTg                                                                               <  4:289719/149‑1 (MQ=255)
                                                                     cGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAg                                                                            <  4:134901/148‑1 (MQ=255)
                                                                     cGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGt                                                                           <  3:190061/149‑1 (MQ=255)
                                                                           tcgctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                               >  6:87567/1‑123 (MQ=255)
                                                                           tcgctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                               <  5:87567/123‑1 (MQ=255)
                                                                             gctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGaa                                                                   >  7:78041/1‑149 (MQ=255)
                                                                                 gcgcAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   <  7:86616/81‑1 (MQ=255)
                                                                                 gcgcAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   >  8:86616/1‑81 (MQ=255)
                                                                                       ggTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGa                                                                                   <  7:99333/123‑1 (MQ=255)
                                                                                       ggTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGa                                                                                   >  8:99333/1‑123 (MQ=255)
                                                                                           ttGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacg                                                     >  2:294077/1‑149 (MQ=255)
                                                                                             gAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgca                                                   >  3:71954/1‑149 (MQ=255)
                                                                                              aaGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGGGAAGTCGATCGAACGGCCAAAGgcacgcac                                                  >  8:126955/1‑149 (MQ=255)
                                                                                                ggCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                >  6:81456/1‑149 (MQ=255)
                                                                                                  ccg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatg                                              <  6:165062/149‑1 (MQ=255)
                                                                                                   cg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                <  5:208034/146‑1 (MQ=255)
                                                                                                   cg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                >  6:208034/1‑146 (MQ=255)
                                                                                                    g‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatgat                                            <  8:78041/149‑1 (MQ=255)
                                                                                                       ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGAGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGatgatg                                           <  5:201238/149‑1 (MQ=255)
                                                                                                        cTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGt                                          >  4:303045/1‑149 (MQ=255)
                                                                                                                gTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                  >  7:202922/1‑149 (MQ=255)
                                                                                                                 tACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtgg                                 >  4:306187/1‑149 (MQ=255)
                                                                                                                    gccgcGGATGTCAACGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCCACAATTTGGAGTCGAGGATAGAGGTGGTTGTTTGTGAGGTGTATTAGAAAAAACAGCCAAAAGCACGAACCAAGATGTTGCAGGTGGcca                              >  8:108052/1‑149 (MQ=255)
                                                                                                                        cgGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAg                          <  3:306187/149‑1 (MQ=255)
                                                                                                                         ggATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAgg                         <  4:71954/149‑1 (MQ=255)
                                                                                                                               cACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGGTCGAGGTGGTGGGGTGGGAGGTGAAGTCGATCGCACGGCCAAAGGCACGCACGATGATGTTTCTGGTGGCCCCCAgagatcgg                   >  2:319951/1‑142 (MQ=255)
                                                                                                                               cACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccacca                           <  1:319951/141‑1 (MQ=255)
                                                                                                                                aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                  >  3:26333/1‑149 (MQ=255)
                                                                                                                                aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                  >  4:298975/1‑149 (MQ=255)
                                                                                                                                     tGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTgg             <  8:202922/149‑1 (MQ=255)
                                                                                                                                        ccacgacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc           >  5:333474/4‑148 (MQ=255)
                                                                                                                                            gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt      >  6:241502/1‑149 (MQ=255)
                                                                                                                                             acgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTAGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTc     >  8:268551/1‑149 (MQ=255)
                                                                                                                                              cgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCg    <  5:241502/149‑1 (MQ=255)
                                                                                                                                               gaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCgg   >  1:22148/1‑149 (MQ=255)
                                                                                                                                                aaGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGc  <  5:81456/149‑1 (MQ=255)
                                                                                                                                                    ttGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtg                    <  7:87071/127‑1 (MQ=255)
                                                                                                                                                    ttGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtg                    >  8:87071/1‑127 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGC  >  NC_002947/160370‑160660

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: