Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,458,747 A→T 6.9% intergenic (+39/‑229) pcaJ → / → trkH‑I 3‑oxoadipate CoA‑transferase subunit A/potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,458,7470AT 6.9% 165.2 / 3.0 58intergenic (+39/‑229)pcaJ/trkH‑I3‑oxoadipate CoA‑transferase subunit A/potassium transporter
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base T (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCC  >  NC_002947/4458603‑4458894
                                                                                                                                                |                                                                                                                                                     
ggCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatc                                                                                                                                                    >  6:357059/1‑148 (MQ=255)
 gCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatcaa                                                                                                                                                  >  4:108009/1‑149 (MQ=255)
 gCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatca                                                                                                                                                   <  3:315916/148‑1 (MQ=255)
   cGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatcaa                                                                                                                                                  >  4:77231/1‑147 (MQ=255)
   cGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATa                                                                                                                                                <  8:57772/149‑1 (MQ=255)
          gAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCg                                                                                                                                          >  5:73618/1‑148 (MQ=255)
            aGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                       <  1:93018/149‑1 (MQ=255)
             gTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                       >  2:132477/1‑148 (MQ=255)
             gTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGcac                                                                                                                                       >  1:58006/1‑148 (MQ=255)
                    cTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTATGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCAGCACAAGcg‑ccg                                                                                                                              >  3:139920/1‑149 (MQ=255)
                     tGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcaga                                                                                                                                      <  5:222696/141‑3 (MQ=255)
                     tGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcaga                                                                                                                                      >  6:222696/1‑139 (MQ=255)
                        aGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGa                                                                                                                           <  4:9670/148‑1 (MQ=255)
                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAgg                                                                                                                      >  8:225672/1‑148 (MQ=255)
                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTAGACGAGCTGCAGAAAATCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAgg                                                                                                                      >  1:50344/1‑148 (MQ=255)
                                   tGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCAcg                                                                                                              <  3:77231/149‑1 (MQ=255)
                                     gTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCAcgcg                                                                                                            <  6:73618/149‑1 (MQ=255)
                                        gAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTg                                                                                                         <  3:262379/149‑1 (MQ=255)
                                         aaaTCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCATAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGc                                                                                                        <  1:274289/149‑1 (MQ=255)
                                           aTCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCtt                                                                                                      <  1:329548/149‑1 (MQ=255)
                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTgaga                                                                                                  <  3:117307/148‑1 (MQ=255)
                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                    <  6:185961/146‑1 (MQ=255)
                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                    <  2:58006/146‑1 (MQ=255)
                                                cgcgGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTga                                                                                                    <  1:132477/146‑1 (MQ=255)
                                                         cGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTa                                                                                         >  2:187669/1‑148 (MQ=255)
                                                          gACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACa                                                                                       <  1:187669/149‑1 (MQ=255)
                                                           aCTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAg                                                                                      <  3:78904/149‑1 (MQ=255)
                                                            cTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGc                                                                                     >  2:129677/1‑149 (MQ=255)
                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                               <  4:139920/149‑1 (MQ=255)
                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                               <  5:357059/149‑1 (MQ=255)
                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                               <  4:369212/149‑1 (MQ=255)
                                                                  cGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                               <  7:114376/149‑1 (MQ=255)
                                                                          acttgatCAGCGGCACGCCACTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGa                                                                       >  2:375621/7‑149 (MQ=255)
                                                                            aaaCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATg                                                                     <  7:225672/149‑1 (MQ=255)
                                                                              aCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATgg                                                                    >  3:268447/1‑148 (MQ=255)
                                                                               cTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGaa                                                                  <  4:99369/149‑1 (MQ=255)
                                                                                tCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAAc                                                                 >  7:127263/1‑149 (MQ=255)
                                                                                         tGCCGCTGATCAAGTGATCCAGGAAGACGTAGGGCGTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGgaga                                                        <  5:236512/149‑1 (MQ=255)
                                                                                                gATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAAc                                                 >  3:35253/1‑149 (MQ=255)
                                                                                                 aTCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAAcc                                                >  6:90859/1‑149 (MQ=255)
                                                                                                   cAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACtctc                                                             >  7:287056/1‑134 (MQ=255)
                                                                                                   cAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACtctc                                                             <  8:287056/134‑1 (MQ=255)
                                                                                                     aGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAAt                                            >  4:339633/1‑149 (MQ=255)
                                                                                                     aGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAAt                                            <  2:50344/149‑1 (MQ=255)
                                                                                                        gATCCAGGAAGAGTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGAGTACGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTAAACCCAAt                                            >  1:109222/1‑147 (MQ=255)
                                                                                                        gATCCAGGAAGAATTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGTACGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGATTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGa                                          >  1:114040/1‑149 (MQ=255)
                                                                                                             aGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATc                                    >  6:112180/1‑149 (MQ=255)
                                                                                                                 aGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGc                                <  1:129677/149‑1 (MQ=255)
                                                                                                                   aCTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCgg                              >  3:228088/1‑149 (MQ=255)
                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCaa                                                  <  6:252115/125‑1 (MQ=255)
                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCaa                                                  >  5:252115/1‑125 (MQ=255)
                                                                                                                       aGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAg                           <  3:108009/148‑1 (MQ=255)
                                                                                                                          gCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAg                       <  3:339633/149‑1 (MQ=255)
                                                                                                                           cTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTAt                                     <  6:8165/134‑1 (MQ=255)
                                                                                                                           cTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTAt                                     >  5:8165/1‑134 (MQ=255)
                                                                                                                            tttCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGa                      >  4:26327/1‑148 (MQ=255)
                                                                                                                            tttCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGa                      >  2:179375/1‑148 (MQ=255)
                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGGAATCACGCGGTGCTTGAGCCTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcg                 >  1:283685/1‑149 (MQ=255)
                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                >  7:351475/1‑149 (MQ=255)
                                                                                                                                  tGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcggc               <  4:116386/149‑1 (MQ=255)
                                                                                                                                    cttgtGCGGCACTATTGATTGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                >  7:104472/5‑146 (MQ=255)
                                                                                                                                    cttgtGCGGCACTATTGATTGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcgg                <  8:104472/142‑1 (MQ=255)
                                                                                                                                               caatcaatAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGcc  <  1:375621/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCG‑CCGCGAGCAGGG‑AGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCC  >  NC_002947/4458603‑4458894

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