Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 496,096 T→A 12.5% 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,0960TA12.5% 170.7 / 11.6 69intergenic (‑63/‑33)PP_0408/PP_0409hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base A (3/6);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGG  >  NC_002947/495955‑496242
                                                                                                                                             |                                                                                                                                                  
actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGc                                                                                                                                             >  1:14162/3‑149 (MQ=255)
actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGc                                                                                                                                             >  5:222843/3‑149 (MQ=255)
 aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                            <  1:16203/149‑1 (MQ=255)
  tCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCg                                                                                                                                           <  5:161449/149‑1 (MQ=255)
   cggcggCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                          >  5:223097/1‑149 (MQ=255)
     gcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTca                                                                                                                                        <  3:166423/149‑2 (MQ=255)
      cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACCTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCtt                                                                                                                                       <  2:105661/149‑1 (MQ=255)
      cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCtt                                                                                                                                       <  4:152360/149‑1 (MQ=255)
      cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                       <  6:223097/149‑3 (MQ=255)
      cggcCAGGGCGCCGAGTGGAACGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                       <  2:181793/149‑3 (MQ=255)
       ggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCttt                                                                                                                                      <  3:200177/149‑1 (MQ=255)
           aGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGcc                                                                                                                                   >  2:266694/1‑148 (MQ=255)
              gcgcCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGt                                                                                                                               <  3:57125/149‑1 (MQ=255)
                 ccGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGa                                                                                                                            <  6:247740/149‑1 (MQ=255)
                   gAGTGTTAGGACCGAGGCAAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGTACATGAGTCAACCAAAGGCCAGAATGAAAGCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                          <  6:173520/149‑1 (MQ=255)
                        gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGtata                                                                                                                     <  4:246280/149‑1 (MQ=255)
                          aGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCGAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGa                                                                                                                   <  3:260737/149‑1 (MQ=255)
                             aCCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTCCATGGCCGTCTATGCCAAGTCGAGCGTATAGATg                                                                                                                 >  5:258083/1‑148 (MQ=255)
                              ccGAGGCCAGGCAGCACATTGCTAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGACATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                 >  7:201565/1‑147 (MQ=255)
                              ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                          <  5:90600/122‑1 (MQ=255)
                              ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                          <  5:95153/122‑1 (MQ=255)
                              ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                          >  6:95153/1‑122 (MQ=255)
                              ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                          >  6:90600/1‑122 (MQ=255)
                                gAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGc                                                                                                             >  6:184541/1‑149 (MQ=255)
                                    ccAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCGAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGt                                                                                                         >  7:62782/1‑149 (MQ=255)
                                            cacaTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGcc                                                                                                 >  6:50211/1‑149 (MQ=255)
                                             acaTTGCCAACGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCt                                                                                                >  1:238007/1‑149 (MQ=255)
                                               aTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGt                                                                                              >  5:172984/1‑149 (MQ=255)
                                                ttGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCATAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGtt                                                                                             <  5:235761/149‑1 (MQ=255)
                                                ttGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGtt                                                                                             >  5:49363/1‑149 (MQ=255)
                                                  gCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTg                                                                                            <  6:49363/148‑1 (MQ=255)
                                                   ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGt                                                                                          >  5:75497/1‑149 (MQ=255)
                                                    cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGtt                                                                                         <  1:287222/149‑1 (MQ=255)
                                                     aGCTTCGTGCGATAAAGTGTGGACATCAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTa                                                                                        <  8:62782/149‑1 (MQ=255)
                                                      gCGTCTTGCGATAAAGTGTTGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTAt                                                                                       <  4:294920/149‑1 (MQ=255)
                                                       cgtcgtGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATc                                                                                      >  7:12953/1‑149 (MQ=255)
                                                         tcgtGCGATAAAGTGTGGACATGAGTCACCCCAAGGCCAGAATGAACCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGGTATCgg                                                                                    <  6:228011/149‑1 (MQ=255)
                                                            tGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgctg                                                                                 >  7:38516/1‑149 (MQ=255)
                                                               gATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCtt                                                                              >  2:20356/1‑149 (MQ=255)
                                                                aTAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTTTAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtg                                                                             >  4:105856/1‑149 (MQ=255)
                                                                    aGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGc                                                                         >  1:220513/1‑149 (MQ=255)
                                                                    aGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCAAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGc                                                                         >  8:249898/1‑149 (MQ=255)
                                                                       gtgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGc                                                                      >  3:92071/1‑149 (MQ=255)
                                                                             aTGAGTCAACCCAAGTCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGa                                                                <  6:222843/149‑1 (MQ=255)
                                                                             aTGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGa                                                                <  1:5502/149‑1 (MQ=255)
                                                                                  tCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTc                                                                                                                              >  6:264358/1‑82 (MQ=255)
                                                                                  tCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTc                                                                                                                              <  5:264358/82‑1 (MQ=255)
                                                                                  tCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGcc                                                           >  1:125407/1‑149 (MQ=255)
                                                                                    aaCCCAAGTCCAGAATGAAACCGGCAGCCCAATATTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTg                                                         <  5:45874/149‑1 (MQ=255)
                                                                                     aCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAAGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                        <  1:266694/149‑1 (MQ=255)
                                                                                      cccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTggg                                                       <  2:238007/149‑1 (MQ=255)
                                                                                       ccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                      <  8:12953/149‑1 (MQ=255)
                                                                                          aGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCaa                                                   <  6:172984/149‑1 (MQ=255)
                                                                                            gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                        <  1:88995/142‑1 (MQ=255)
                                                                                            gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                        >  2:88995/1‑142 (MQ=255)
                                                                                              cAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCTAGTGCGTCGCCTGTTTGTTATCTGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGccc                                               >  6:116606/1‑149 (MQ=255)
                                                                                                  attaaaCCGGCAGCCTAATAGTATTTGCCCGTGACGTCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAgg                                           <  8:201565/146‑1 (MQ=255)
                                                                                                      aaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCt                                       <  2:125407/149‑1 (MQ=255)
                                                                                                      aaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGGTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCt                                       <  4:231869/149‑1 (MQ=255)
                                                                                                                cTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgt                              <  7:249898/148‑1 (MQ=255)
                                                                                                                   aTAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggtg                          >  4:224213/1‑149 (MQ=255)
                                                                                                                    tAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGc                         >  8:184754/1‑149 (MQ=255)
                                                                                                                                gTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAcc             <  6:75497/149‑1 (MQ=255)
                                                                                                                                 tGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc              >  4:171571/1‑147 (MQ=255)
                                                                                                                                 tGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACca            <  5:50211/149‑1 (MQ=255)
                                                                                                                                    cGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagca         <  5:184541/149‑1 (MQ=255)
                                                                                                                                     ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcag        <  8:38516/149‑1 (MQ=255)
                                                                                                                                        cATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAg                    <  7:234081/134‑1 (MQ=255)
                                                                                                                                        cATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAg                    >  8:234081/1‑134 (MQ=255)
                                                                                                                                           gctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCATCTGTGGTGCTGGAGCCTGACCAGCAGCAGTgg  >  8:69909/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGG  >  NC_002947/495955‑496242

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: