Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,940 G→T 10.3% intergenic (+72/‑95) 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,9400GT10.3% 180.8 / 10.2 66intergenic (+72/‑95)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base T (4/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG  >  NC_002947/660794‑661087
                                                                                                                                                    |                                                                                                                                                   
aCTGCGCGTGAGCGGACCCAGTATCTGGATTCCTTGCGCCTGTGGTTGAGCCGGTCTTCGTACCTGTCGTCGTAATTGAGGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGaa                                                                                                                                                   <  6:131142/149‑1 (MQ=255)
aCTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGa                                                                                                                                                    <  1:115402/148‑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                                                                                                            <  5:341750/149‑1 (MQ=255)
                                       cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                            >  7:334431/1‑149 (MQ=255)
                                       cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                            >  7:54662/1‑149 (MQ=255)
                                          ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGt                                                                                                          <  1:358971/148‑1 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                  <  6:314035/149‑1 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                  >  4:141905/1‑149 (MQ=255)
                                                 cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                  >  6:117637/1‑149 (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:345180/1‑149 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                  >  7:104727/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                                                 >  3:229191/1‑108 (MQ=255)
                                                                                                                                             tGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCaa                                                 <  4:229191/108‑1 (MQ=255)
                                                                                                                                             tGAGGGCTTCGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGa        <  1:298670/149‑1 (MQ=255)
                                                                                                                                             tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGa        <  3:141905/149‑1 (MQ=255)
                                                                                                                                              gAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGaa       >  8:55980/1‑149 (MQ=255)
                                                                                                                                               aGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAAt      >  2:266235/1‑149 (MQ=255)
                                                                                                                                               aGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAAt      <  1:386049/149‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAA‑‑GCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG  >  NC_002947/660794‑661087

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: