Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,942 A→C 9.4% intergenic (+74/‑93) speA → / → PP_0568 biosynthetic arginine decarboxylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,9420AC 9.4% 185.5 / 10.3 64intergenic (+74/‑93)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base C (4/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGGC  >  NC_002947/660794‑661089
                                                                                                                                                      |                                                                                                                                                   
aCTGCGCGTGAGCGGACCCAGTATCTGGATTCCTTGCGCCTGTGGTTGAGCCGGTCTTCGTACCTGTCGTCGTAATTGAGGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGaa                                                                                                                                                     <  6:131142/149‑1 (MQ=255)
 cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAg                                                                                                                                                    <  7:30445/149‑1 (MQ=255)
     gcgTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCt                                                                                                                                                >  6:97007/1‑149 (MQ=255)
          aGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc                                                                                                                                             >  6:150847/1‑147 (MQ=255)
              gACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGt                                                                                                                                       >  2:102329/1‑149 (MQ=255)
                cccAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTca                                                                                                                                     <  7:170042/149‑1 (MQ=255)
                 ccAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTcc                                                                                                                                     >  7:56936/1‑147 (MQ=255)
                    gTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAAc                                                                                                                                 <  5:169068/149‑1 (MQ=255)
                    gTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGAGTCAAAAc                                                                                                                                 >  8:364197/1‑149 (MQ=255)
                            aTGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAAATGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaaa                                                                                                                          >  1:84604/1‑148 (MQ=255)
                                cTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCTGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATcc                                                                                                                     >  1:341364/1‑149 (MQ=255)
                                cTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATcc                                                                                                                     <  7:156681/149‑1 (MQ=255)
                                    cgcCTGGGCTTTACCCGGTCTTTTTACCTGTCTTCTTTATTTACGGGGCCTCTCCTGCAATACCCCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCg                                                                                                                    <  1:125642/148‑1 (MQ=255)
                                    cgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGcc                                                                                                                 <  1:102329/149‑1 (MQ=255)
                                     gcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGccc                                                                                                                <  4:358777/149‑1 (MQ=255)
                                       cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                              >  7:334431/1‑149 (MQ=255)
                                       cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                              <  5:341750/149‑1 (MQ=255)
                                       cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                              >  7:54662/1‑149 (MQ=255)
                                          ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGt                                                                                                            <  1:358971/148‑1 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                    >  4:141905/1‑149 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                    >  6:117637/1‑149 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                    <  6:314035/149‑1 (MQ=255)
                                                  ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAc                                                                                                   <  6:18095/149‑1 (MQ=255)
                                                      tCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAAt                                                                                               <  8:334431/149‑1 (MQ=255)
                                                          cGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGcc                                                                                            <  2:341364/148‑1 (MQ=255)
                                                          cGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGcc                                                                                            <  5:97007/148‑1 (MQ=255)
                                                           gTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGc                                                                                          >  7:87587/1‑149 (MQ=255)
                                                              ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGTGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAc                                                                                       <  4:58538/149‑1 (MQ=255)
                                                              ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAc                                                                                       <  1:149953/149‑1 (MQ=255)
                                                                tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                                                                     >  6:169135/1‑149 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                    >  7:104727/1‑149 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                    >  7:345180/1‑149 (MQ=255)
                                                                        tAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGc                                                                             <  7:167579/149‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAA‑‑GCCTGTGCTTCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCa                                                                     >  1:48713/1‑148 (MQ=255)
                                                                                  ggCCGCTCCTGCAAAA‑‑GCCTGTGCTTCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                   <  2:48713/149‑1 (MQ=255)
                                                                                  ggCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCCCCAATGCCGCTACCTTTGAATGCAATGTTCa                                                                     >  1:309694/1‑147 (MQ=255)
                                                                                  ggCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                   <  8:54662/149‑1 (MQ=255)
                                                                                            gCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGa                                                          >  4:267664/1‑148 (MQ=255)
                                                                                             cAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGata                                                        >  3:150106/1‑149 (MQ=255)
                                                                                              aaaa‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGatat                                                       >  2:386049/1‑149 (MQ=255)
                                                                                                        gtgCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTc                                               >  6:138123/1‑149 (MQ=255)
                                                                                                         tgctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTcc                                              <  3:267664/149‑1 (MQ=255)
                                                                                                              ccTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg                                           <  8:345180/147‑1 (MQ=255)
                                                                                                                   tgtgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCa                                     <  5:138123/148‑1 (MQ=255)
                                                                                                                     tgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCg                                  >  8:148423/1‑149 (MQ=255)
                                                                                                                       ggAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGtt                                <  1:377414/149‑1 (MQ=255)
                                                                                                                        tcgCGGCATTGCGTCGCGATATGATGGCTAAGCCATCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGttt                               <  3:194601/147‑1 (MQ=255)
                                                                                                                          gCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTTGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTAc                             <  7:364197/149‑1 (MQ=255)
                                                                                                                               cTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGc                        >  2:47673/1‑149 (MQ=255)
                                                                                                                                ttGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCTGCCCTGGTCCATCACAACAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTACCAACGATTACGTcgcg                       >  5:270928/1‑149 (MQ=255)
                                                                                                                                    gTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTc                    >  2:248915/1‑148 (MQ=255)
                                                                                                                                        gctggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAgtg               >  2:298670/6‑149 (MQ=255)
                                                                                                                                          tggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCTCCGTCCAGGCGCCCATCGTTTACGTGGCGGTCGAGTGTc             >  6:306136/4‑149 (MQ=255)
                                                                                                                                          tggtGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTAAGTGGCGGTCTAgtt               >  6:311360/4‑146 (MQ=255)
                                                                                                                                             tGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCaa                                                   <  4:229191/108‑1 (MQ=255)
                                                                                                                                             tGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCaa                                                   >  3:229191/1‑108 (MQ=255)
                                                                                                                                             tGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGa          <  1:298670/149‑1 (MQ=255)
                                                                                                                                             tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGa          <  3:141905/149‑1 (MQ=255)
                                                                                                                                              gAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGaa         >  8:55980/1‑149 (MQ=255)
                                                                                                                                               aGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAAt        <  1:386049/149‑1 (MQ=255)
                                                                                                                                               aGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAAt        >  2:266235/1‑149 (MQ=255)
                                                                                                                                                 ggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATcc      >  4:326784/1‑149 (MQ=255)
                                                                                                                                                  gCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCg     <  7:55980/149‑1 (MQ=255)
                                                                                                                                                  gCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCg     <  8:87587/149‑1 (MQ=255)
                                                                                                                                                   cGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCACACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCgg    <  4:254471/149‑1 (MQ=255)
                                                                                                                                                     aaGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGgc  <  2:84604/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGGC  >  NC_002947/660794‑661089

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