Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 496,096 T→A 7.8% intergenic (‑63/‑33) PP_0408 ← / → PP_0409 hypothetical protein/two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947496,0960TA 7.8% 175.4 / 3.4 61intergenic (‑63/‑33)PP_0408/PP_0409hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (3/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAG  >  NC_002947/495954‑496239
                                                                                                                                              |                                                                                                                                               
gCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATgg                                                                                                                                           >  7:38947/1‑149 (MQ=255)
  aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                         >  7:206047/1‑149 (MQ=255)
  aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                         >  2:236028/1‑149 (MQ=255)
     ggcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTc                                                                                                                                      <  7:48534/149‑1 (MQ=255)
       cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                    <  2:334934/149‑3 (MQ=255)
                   cGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAg                                                                                                                        >  5:239209/1‑149 (MQ=255)
                    gAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                       <  4:100913/149‑1 (MQ=255)
                         gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGta                                                                                                                    >  2:229888/1‑147 (MQ=255)
                          aaGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGta                                                                                                                    >  6:140106/1‑146 (MQ=255)
                          aaGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAg                                                                                                                 >  7:169309/1‑149 (MQ=255)
                                  aGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct                                                                                                                                                 <  6:319314/109‑1 (MQ=255)
                                  aGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct                                                                                                                                                 >  5:319314/1‑109 (MQ=255)
                                   ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                        <  5:122985/149‑1 (MQ=255)
                                     ccAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGt                                                                                                      <  1:283005/149‑1 (MQ=255)
                                       aGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                    >  3:107901/1‑149 (MQ=255)
                                                aTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAAGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGt                                                                                           >  2:116173/1‑149 (MQ=255)
                                                 ttGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGtt                                                                                          <  2:289537/149‑1 (MQ=255)
                                                    ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTc                                                                                                 >  7:251167/1‑139 (MQ=255)
                                                    ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTc                                                                                                 <  8:251167/139‑1 (MQ=255)
                                                    ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTgg                                                                                        <  8:206047/148‑1 (MQ=255)
                                                     cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGtt                                                                                      <  2:251669/149‑1 (MQ=255)
                                                     cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGtt                                                                                      >  8:238682/1‑149 (MQ=255)
                                                     cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGtt                                                                                      >  1:183950/1‑149 (MQ=255)
                                                       gCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTAt                                                                                    <  2:183950/149‑1 (MQ=255)
                                                         gtcgtGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCg                                                                                  >  6:349556/1‑149 (MQ=255)
                                                             tGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGGATGGCCGACTTAGCCAAGTCGAGCGTATAGAAGAAGCGAGTGCGTCGCCAGTTGGTTATCGgctg                                                                              >  7:48189/1‑149 (MQ=255)
                                                               cGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgctgct                                                                            >  2:59266/1‑149 (MQ=255)
                                                                    aaGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAATCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTgg                                                                       <  1:59266/149‑1 (MQ=255)
                                                                    aaGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTgg                                                                       <  6:239209/149‑1 (MQ=255)
                                                                         tgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGc                                                                   <  7:238682/148‑1 (MQ=255)
                                                                                 aGTCAACCCAAGGCCAGAATGAAACCGGCCGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGAtt                                                           <  5:274440/148‑1 (MQ=255)
                                                                                 aGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGAtt                                                           <  1:229888/148‑1 (MQ=255)
                                                                                  gTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGc                                                         >  5:135644/1‑149 (MQ=255)
                                                                                     aaCCCAAGTCCAGAATGAAACCGGCAGCCTAATAGTATTTGCGCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCt                                                       <  8:169309/148‑1 (MQ=255)
                                                                                       cccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTggg                                                    <  6:135644/149‑1 (MQ=255)
                                                                                        ccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                   >  5:217212/1‑149 (MQ=255)
                                                                                        ccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                   <  2:185883/149‑1 (MQ=255)
                                                                                               cAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                    >  3:198866/1‑93 (MQ=255)
                                                                                               cAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                    <  4:198866/93‑1 (MQ=255)
                                                                                                    tGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                       >  1:240432/1‑149 (MQ=255)
                                                                                                    tGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                       >  3:128056/1‑149 (MQ=255)
                                                                                                    tGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                       >  6:323368/1‑149 (MQ=255)
                                                                                                      aaaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgcgc                                     >  7:123701/1‑149 (MQ=255)
                                                                                                      aaaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgcgc                                     >  3:99991/1‑149 (MQ=255)
                                                                                                               gcgaaaTAGTATTTGCCCGTCACGGCGTTGTTGCATGGCCGTCTTTGCGACTTCGAGCGTATACATGAAGCGAGTGCGTCGCCTGTTGGTTATGGGCTGCTTCTGGCTGCCCTTGATTGCCTGGGCCAAGGCCGAGGCGCTGCCctctc                            <  8:48189/145‑2 (MQ=255)
                                                                                                                   aaTAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggt                        >  3:185941/1‑149 (MQ=255)
                                                                                                                   aaTAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTAGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggt                        >  1:82900/1‑149 (MQ=255)
                                                                                                                   aaTAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTAGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggt                        >  1:87842/1‑149 (MQ=255)
                                                                                                                     tAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGc                      >  6:271260/1‑149 (MQ=255)
                                                                                                                       gTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                    <  4:107901/149‑1 (MQ=255)
                                                                                                                        tATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                    >  5:321054/1‑148 (MQ=255)
                                                                                                                         aTTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTgg                   <  5:68215/148‑1 (MQ=255)
                                                                                                                         aTTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTgg                   <  5:140106/148‑1 (MQ=255)
                                                                                                                         aTTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                    >  4:202003/1‑147 (MQ=255)
                                                                                                                          tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGa                  <  8:38947/148‑1 (MQ=255)
                                                                                                                            tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                <  4:99991/148‑1 (MQ=255)
                                                                                                                             gCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGcc               <  6:321054/148‑1 (MQ=255)
                                                                                                                               ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGa            <  3:202003/149‑1 (MQ=255)
                                                                                                                                aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc           >  1:344432/4‑149 (MQ=255)
                                                                                                                                 gTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAcc          <  2:240432/149‑1 (MQ=255)
                                                                                                                                 aagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc           >  7:24644/3‑148 (MQ=255)
                                                                                                                                    aCGGCCATGCTGCATGGCCGTCTTTGCCACGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagc       <  5:323368/149‑1 (MQ=255)
                                                                                                                                         cATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagcag  <  1:236028/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAG  >  NC_002947/495954‑496239

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