Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 362,933 T→C 6.4% K184E (AAG→GAG)  lcdH ← L‑carnitine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947362,9330TC 6.4% 182.5 / 2.2 62K184E (AAG→GAG) lcdHL‑carnitine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base C (2/2);  total (28/34)
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.56e-01

CCATGAACGACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGCCGACGACGCAG  >  NC_002947/362789‑363075
                                                                                                                                                |                                                                                                                                              
ccATGAACGACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTgcg                                                                                                                                             >  7:277052/1‑148 (MQ=255)
  aTGAACGACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGca                                                                                                                                           >  1:14626/1‑148 (MQ=255)
     aaCGACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACAt                                                                                                                                        <  2:14626/148‑1 (MQ=255)
        gACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTAGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAg                                                                                                                                    >  8:172236/1‑149 (MQ=255)
           cAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCgg                                                                                                                                 <  6:321514/149‑1 (MQ=255)
           cAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCACTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCgg                                                                                                                                 <  6:307980/149‑1 (MQ=255)
             gcgcAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTCGCGCACATGCAGCGGgc                                                                                                                               >  2:51394/1‑149 (MQ=255)
               gcAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGgcgc                                                                                                                             <  3:155008/149‑1 (MQ=255)
               gcAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGgcg                                                                                                                              >  7:130137/1‑148 (MQ=255)
                 aaGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCAt                                                                                                                           >  7:274069/1‑149 (MQ=255)
                     cTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGcc                                                                                                                        >  1:165553/1‑148 (MQ=255)
                      tGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGcc                                                                                                                        >  7:102539/1‑147 (MQ=255)
                        cgcCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGa                                                                                                                      <  1:141048/147‑1 (MQ=255)
                         gcCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTg                                                                                                                   <  2:165553/149‑1 (MQ=255)
                         gcCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTg                                                                                                                   >  1:129888/1‑149 (MQ=255)
                             aaGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGgt                                                                                                               >  3:292868/1‑149 (MQ=255)
                             aaGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTCGCGCACATGCAGCGGGCGCATGCCGAGTGCGgt                                                                                                               <  1:51394/149‑1 (MQ=255)
                               gCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGgtgt                                                                                                             <  8:277052/149‑1 (MQ=255)
                               gCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTCGCGCACATGCAGCGGGCGCATGCCGAGTGCGgtg                                                                                                              >  6:316874/1‑148 (MQ=255)
                                 gAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAg                                                                                                           <  8:102539/149‑1 (MQ=255)
                                  aaTCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGa                                                                                                          <  8:274069/149‑1 (MQ=255)
                                         tcgtcgATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTTGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGtttt                                                                                                   <  7:172236/149‑1 (MQ=255)
                                           gtcgATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTg                                                                                                  <  4:292868/148‑1 (MQ=255)
                                                  tCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCtt                                                                                           <  8:130137/148‑1 (MQ=255)
                                                    gCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCaa                                                                                        <  3:157803/149‑1 (MQ=255)
                                                    gCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCaa                                                                                        >  4:157803/1‑149 (MQ=255)
                                                    gCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCaa                                                                                        >  6:212149/1‑149 (MQ=255)
                                                    gCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCaa                                                                                        <  5:125760/149‑1 (MQ=255)
                                                     cccGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCaa                                                                                        >  4:292002/1‑148 (MQ=255)
                                                      ccGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATa                                                                                      <  1:306133/149‑1 (MQ=255)
                                                      ccGTGGTGGCAACTCCGTCGTTGACCAGGGGCAAGGCCTCGCGCCATAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGGAATTCCTTGCGCACATGCATCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATa                                                                                      <  6:261280/149‑1 (MQ=255)
                                                           gtggCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGTGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTc                                                                                 >  8:48573/1‑149 (MQ=255)
                                                              gCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTCGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCggg                                                                              <  5:316874/149‑1 (MQ=255)
                                                                aaCTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGgtgt                                                                                                             <  8:128028/116‑1 (MQ=255)
                                                                aaCTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGgtgt                                                                                                             >  7:128028/1‑116 (MQ=255)
                                                                  cTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAg                                                                          <  2:129888/149‑1 (MQ=255)
                                                                       tcgtTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCg                                                                     <  3:292002/149‑1 (MQ=255)
                                                                         gtTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCgg                                                                    <  2:157778/148‑1 (MQ=255)
                                                                               cAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTAccgc                                                             >  4:305801/1‑149 (MQ=255)
                                                                               cAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTAccgc                                                             >  1:231421/1‑149 (MQ=255)
                                                                                      aagcacTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGAAGAAGCCTGGCACTTCCTTGCGCACATGCAGCGCGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGAt                                                      >  5:46403/6‑149 (MQ=255)
                                                                                      aaGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGAt                                                      <  3:305801/149‑1 (MQ=255)
                                                                                              gcgcCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGAGGTGTAGATGGTTTTTGCAGCTTCAATAGCATCGCGTCAGGTGCGACTACAGCCCACGATATcaac                                                >  7:294659/1‑147 (MQ=255)
                                                                                                  cAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTgg                                                                                                 <  4:2077/94‑1 (MQ=255)
                                                                                                  cAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTgg                                                                                                 >  3:2077/1‑94 (MQ=255)
                                                                                                          ttCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTcc                                                  >  3:52043/1‑133 (MQ=255)
                                                                                                          ttCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTcc                                                  <  4:52043/133‑1 (MQ=255)
                                                                                                                   gCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCCCATGCATCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATATTTTTTTTTTATGTTAGGCTTACGACGTAGATCTACACCAGCGGCAGCAGGTAAACTGGGt                         <  2:255623/149‑1 (MQ=255)
                                                                                                                   gCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATATTTTCTTTTTATGTTCGGCTACCGCCGAAGATCTCCACCAGCGGCAGCAGGTAAACTGGGt                         <  2:237297/149‑1 (MQ=255)
                                                                                                                       tCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTg                                                                            <  6:45198/94‑1 (MQ=255)
                                                                                                                       tCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTg                                                                            >  5:45198/1‑94 (MQ=255)
                                                                                                                            gATGAAGCGTGGAACTTCCTTGCGCTCATGCAGCGGGCGTATGCCGAGTGGGGTGTAGATGGTTTTGGCAGCTTCAATAGTTTTTGTTTATGTTCGACTTACGACGAAGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGt                <  8:294659/149‑1 (MQ=255)
                                                                                                                               gAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGc             >  2:316425/1‑149 (MQ=255)
                                                                                                                                 aGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGt                >  8:162569/1‑144 (MQ=255)
                                                                                                                                 aGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGt                <  7:162569/144‑1 (MQ=255)
                                                                                                                                  gCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGCcga          <  5:212149/149‑1 (MQ=255)
                                                                                                                                      ggCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTggg                          <  2:83137/129‑1 (MQ=255)
                                                                                                                                      ggCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGGTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTggg                          >  1:83137/1‑129 (MQ=255)
                                                                                                                                        cagttCCTTGCGCTCTTCGAGCGGGCGCTTGCGGATGGGTGTGCTGATGCTTTTGGCAGCTTCATTAGCCTGTGTTGAGGTGAGGCTACCGCCTACGATCTCCAGCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGCCGACGACGc    <  6:46403/146‑1 (MQ=255)
                                                                                                                                          cTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCATCAATAGCATCGGGTGAGGTGCGGCTACCGCCGCCGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGCCGACGACGCAg  >  6:178039/1‑149 (MQ=255)
                                                                                                                                           ttCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTccacca                                              >  2:181902/1‑104 (MQ=255)
                                                                                                                                           ttCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGATTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTccacca                                              <  1:181902/104‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CCATGAACGACCAGCGCAAGCCTGCGCCAAAGCGAATCGCGTCGTCGATTTCGCCCGTGGTGGCAACTCCGTCGTTGACCAGGTGCAAGGCCTCGCGCCAGAGCGCTTCGAGCAGGCGGTCGGCGATGAAGCCTGGCACTTCCTTGCGCACATGCAGCGGGCGCATGCCGAGTGCGGTGTAGATGGTTTTGGCAGCTTCAATAGCCTCGGGTGAGGTGCGGCTACCGCCGACGATCTCCACCAGCGGCAGCAGGTAAACTGGGTTGAACGGGTGGCCGACGACGCAG  >  NC_002947/362789‑363075

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