Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,731,502 C→A 7.7% intergenic (‑59/+35) lpdG ← / ← sucB dihydrolipoyl dehydrogenase/2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,731,5020CA 7.7% 173.9 / 6.4 66intergenic (‑59/+35)lpdG/sucBdihydrolipoyl dehydrogenase/2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (3/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-02

GTCTTCAGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAG  >  NC_002947/4731355‑4731647
                                                                                                                                                      |                                                                                                                                                  
gTCTTCAGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTCTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGc                                                                                                                                                    <  1:303311/149‑1 (MQ=255)
  cTTCAGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGaa                                                                                                                                                 <  6:17967/149‑1 (MQ=255)
      aGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCa                                                                                                                                             >  3:302180/1‑149 (MQ=255)
      aGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCGTCATGGGAAGCCa                                                                                                                                             >  4:48207/1‑149 (MQ=255)
      aGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGAACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTACTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGGATCTTTGGCAGCCa                                                                                                                                             >  4:401059/1‑149 (MQ=255)
        aCCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAgg                                                                                                                                           >  2:96097/1‑149 (MQ=255)
         ccAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGc                                                                                                                                          <  1:252621/149‑1 (MQ=255)
               tGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAg                                                                                                                                            >  1:303693/1‑141 (MQ=255)
               tGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAg                                                                                                                                            <  2:303693/141‑1 (MQ=255)
                 ggCAGCCTTGATTGCAGACACATAGCCGCGTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCAca                                                                                                                                  <  8:328539/149‑2 (MQ=255)
                 ggCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGAGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGTGCCGGGCCTGGCTTCGCAAGAAGCCAGGccgcgc                                                                                                                                     <  7:342136/148‑5 (MQ=37)
                        ttGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCa                                                                                                                          >  8:127143/1‑149 (MQ=255)
                            tGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGcc                                                                                                                       <  8:277175/148‑1 (MQ=255)
                                aGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAgg                                                                                                                  <  7:505/149‑1 (MQ=255)
                                aGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAgg                                                                                                                  <  4:302180/149‑1 (MQ=255)
                                             ccTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCcag                                                                                                     <  5:300446/149‑1 (MQ=255)
                                                  gCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCcagcagca                                                                                                <  7:127143/149‑1 (MQ=255)
                                                   ccTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCcagcagcag                                                                                               >  6:402374/1‑149 (MQ=255)
                                                       cACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGc                                                                                             <  1:103168/147‑1 (MQ=255)
                                                         ccGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGa                                                                                         >  6:335441/1‑149 (MQ=255)
                                                           gATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACg                                                                                       >  1:5454/1‑149 (MQ=255)
                                                                ccACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCTAACCGCGC‑‑GGCCGGGCGTGTTTCGCAGCCACGCCCCGGACCGGACGCCCCCACAGGGGTTTATTTGCAACAGCACGCCGGCCGGCTcg                                                                                  >  6:34461/1‑148 (MQ=255)
                                                                  aCTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTc                                                                                <  5:402374/149‑1 (MQ=255)
                                                                   cTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGCTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTc                                                                                <  1:96097/149‑1 (MQ=255)
                                                                   cTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTc                                                                                >  4:251898/1‑148 (MQ=255)
                                                                       gTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCAACGCCCGGCCCGGTACGCCCCCAGATGGTTAGATTTGCAGCAGCAAGCGGGACGGATCTTCcagca                                                                           >  7:204383/1‑149 (MQ=255)
                                                                          gAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAggg                                                                        >  2:324152/1‑148 (MQ=255)
                                                                             tttCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCt                                                                     <  8:248338/149‑1 (MQ=255)
                                                                               tCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCtt                                                                    <  5:335441/148‑1 (MQ=255)
                                                                                   ggTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCg                                                                                                                                       <  1:395098/78‑1 (MQ=255)
                                                                                   ggTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCcagca                                                                           <  3:335036/137‑1 (MQ=255)
                                                                                   ggTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCcagca                                                                           >  4:335036/1‑137 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGtt                                                                       >  7:149347/1‑136 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGtt                                                                       <  8:149347/136‑1 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGAc                                                           >  8:219649/1‑148 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGAc                                                           >  5:389774/1‑148 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTACAGCAGCACGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGAc                                                           >  6:341487/1‑148 (MQ=255)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCGCAAGAAGCCAg                                                                                                                                            >  8:93275/1‑68 (MQ=37)
                                                                                        tAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCGCAAGAAGCCAg                                                                                                                                            <  7:93275/68‑1 (MQ=37)
                                                                                         aaaaGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCtt                                                                    >  3:303140/1‑138 (MQ=255)
                                                                                         aaaaGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCtt                                                                    <  4:303140/138‑1 (MQ=255)
                                                                                         aaaaGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGAc                                                           >  7:53664/1‑147 (MQ=255)
                                                                                          aaaGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAg                                                        >  8:317913/1‑149 (MQ=255)
                                                                                              aTTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAg                                                    >  1:91964/1‑149 (MQ=255)
                                                                                                         agctacaagctacaAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCtt                                          >  3:257524/1‑149 (MQ=255)
                                                                                                         agctacaagctacaAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCt                                           <  1:206460/148‑1 (MQ=255)
                                                                                                              caagctacaagcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCtt                                      <  7:219649/148‑1 (MQ=255)
                                                                                                                 gctacaagcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGcc                                   <  3:259576/148‑1 (MQ=255)
                                                                                                                 gctacaagcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACCGGACGCCGCCAGAGGGGTAGATGTCCAGCACCAGGCCGGCCGGAGATTCCCGCACGGTCGTGATAGGTACCAGGAAGGTTACCGCTTCCGTGCCg                                  >  4:259576/1‑149 (MQ=255)
                                                                                                                     caagcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGa                              <  2:5454/149‑1 (MQ=255)
                                                                                                                     caagcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGa                              >  4:385620/1‑149 (MQ=255)
                                                                                                                       agcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATc                            <  7:317913/149‑1 (MQ=255)
                                                                                                                        gcCGCGAACCGCGC‑‑GGCCTGGCTTCTTTCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGCCGGATCTTCCAGCAGGTTCTTGATAGTGACAAGGAAGGTTACCGCTTCCTTGCCGTCGCTcc                           >  6:382020/1‑148 (MQ=255)
                                                                                                                        gcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCca                                                                                                      <  6:73219/74‑1 (MQ=255)
                                                                                                                        gcCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCca                                                                                                      >  5:73219/1‑74 (MQ=255)
                                                                                                                           gcgAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCCCGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGc                        >  3:82316/1‑149 (MQ=255)
                                                                                                                             gAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCgg                      <  8:53664/149‑1 (MQ=255)
                                                                                                                              aaCCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGt                     <  3:22359/149‑1 (MQ=255)
                                                                                                                                 cgcgc‑‑gGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACTAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGAt                  <  3:251898/149‑1 (MQ=255)
                                                                                                                                    gc‑‑gGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCCGCCAGCGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTGGATCAGGCGGTGATCg                >  6:91000/1‑148 (MQ=255)
                                                                                                                                    gc‑‑gGCCTGGATTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGt               <  7:328452/149‑1 (MQ=255)
                                                                                                                                     c‑‑gGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTa              <  6:389774/149‑1 (MQ=255)
                                                                                                                                         gCCTGGCTTCGCAAGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTAc             >  7:322475/1‑148 (MQ=255)
                                                                                                                                          ccTGGCTTCGCAAGAAGCCAGGCC‑‑GCGCGGTTCGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGa           >  5:393159/1‑147 (MQ=38)
                                                                                                                                                   ttGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAg  >  2:139129/1‑149 (MQ=255)
                                                                                                                                                   ttGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAg  >  5:371820/1‑149 (MQ=255)
                                                                                                                                                     gCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCa   <  5:382020/146‑1 (MQ=255)
                                                                                                                                                      cGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTg                                     <  3:384548/111‑1 (MQ=255)
                                                                                                                                                      cGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTg                                     >  4:384548/1‑111 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTCTTCAGACCAAGTTGGGCAGCCTTGATGGCAGCCACATAGCCGCCTGGGCCTGCACCGATTACCACTACGTCGAATTTCTGGGTCATAAAAGATTCCTTATCAGCTACAAGCTACAAGCCGCGAACCGCGC‑‑GGCCTGGCTTCTTGCGAAGCCAGGCCCGGCACGGTACGCAGCCAGAGGGTTAGATGTCCAGCAGCAGGCGGGACGGATCTTCCAGCAGGTTCTTGATAGTGACCAGGAAGGTTACCGCTTCCTTGCCGTCGATCAGGCGGTGATCGTACGACAGCGCCAG  >  NC_002947/4731355‑4731647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: