Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,452,259 C→T 5.9% G462R (GGG→AGG)  paaZ ← oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,452,2590CT 5.9% 216.4 / 3.3 68G462R (GGG→AGG) paaZoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (37/27);  new base T (2/2);  total (39/29)
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.55e-01

CTGATAGCGGCAAGCATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTT  >  NC_000913/1452114‑1452402
                                                                                                                                                 |                                                                                                                                               
cTGATAGCGGCAAGCATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCgg                                                                                                                                               <  4:61955/148‑1 (MQ=255)
   aTAGCGGCAAGCATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGt                                                                                                                                              <  6:127048/146‑1 (MQ=255)
         gCAAGCATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCtt                                                                                                                                      <  3:60894/148‑1 (MQ=255)
               aTCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCg                                                                                                                                <  1:359998/148‑1 (MQ=255)
                tcgtcgGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCTGGTGGATTCTTTTGCCGa                                                                                                                               >  7:128305/1‑148 (MQ=255)
                          aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCAtt                                                                                                                      >  1:125827/1‑147 (MQ=255)
                          aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCAtt                                                                                                                      >  5:343743/1‑147 (MQ=255)
                          aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCAtt                                                                                                                      >  8:26985/1‑147 (MQ=255)
                          aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCAtt                                                                                                                      >  8:20558/1‑147 (MQ=255)
                           cccTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGa                                                                                                                    <  5:32510/148‑1 (MQ=255)
                            ccTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAg                                                                                                                   >  2:332196/1‑148 (MQ=255)
                             cTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAgg                                                                                                                  >  8:315140/1‑148 (MQ=255)
                               gAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                >  4:354560/1‑148 (MQ=255)
                                    gCGGTTCGCTGCATGTAATGTTTCACCGCGCGTAAACCGCCTAATTCTTCACCGCCGCCTGCGCCCCCGGGCCCCCCATGTACCAGTTGTGGCAGTGGGCAGCCATGCCTGGGGGATAGTTTTGCGGACTCTTCATGGCGGAGCGGaa                                                                                                           >  6:374076/1‑148 (MQ=255)
                                     cGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAAt                                                                                                          <  1:332196/148‑1 (MQ=255)
                                        ttCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                       <  8:58575/148‑1 (MQ=255)
                                         tCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATc                                                                                                               <  4:81928/139‑1 (MQ=255)
                                         tCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATc                                                                                                               >  3:81928/1‑139 (MQ=255)
                                          cGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGc                                                                                                      <  6:74848/147‑1 (MQ=255)
                                          cGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGc                                                                                                      <  2:125827/147‑1 (MQ=255)
                                              gCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTGATTGAGGATCTGAATTCGCCCATg                                                                                                 >  2:341683/1‑148 (MQ=255)
                                                 tGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCg                                                                                               >  5:406096/1‑147 (MQ=255)
                                                 tGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCg                                                                                               >  2:45640/1‑147 (MQ=255)
                                                   tAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTa                                                                                             <  7:26985/147‑1 (MQ=255)
                                                   tAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTa                                                                                             <  7:20558/147‑1 (MQ=255)
                                                    aaTGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                >  1:102337/1‑127 (MQ=255)
                                                    aaTGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                <  2:102337/127‑1 (MQ=255)
                                                    aaTGTATCACCGCTCGTAAACCGCCTAATTCTCCACCGCCTCCTGCGCGACAACGCCCACCATGTACCAGTTGTGCCAGTCGGGAGCCAGGCCCGGTGGATTCTTTAGCCGACTCTTCATTGAGGATCGGAATTCGCGCAAGCGTACg                                                                                           >  7:189141/1‑148 (MQ=255)
                                                          tCACCGCTCGTAAACCGCCTATTTCTTCACCGCCTCATGCGCGACAAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcg                                                                                      <  6:343743/147‑1 (MQ=255)
                                                           cACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgcg                                                                                    >  2:153769/1‑148 (MQ=255)
                                                            aCCGCTCGTAATCCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGGCTCGTCATTGAGGAGCTGAATTCGCCCATGCGGACGTGACGCGt                                                                                   >  2:127424/1‑148 (MQ=255)
                                                             ccGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTc                                                                                  >  1:271335/1‑148 (MQ=255)
                                                               gCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCgg                                                                                >  5:83805/1‑148 (MQ=255)
                                                               gCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCgg                                                                                >  5:293437/1‑148 (MQ=255)
                                                                cTCGTAAACCGCCTTATTCTTCACGGCCTCATGCGCGACCTGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgc                                                                                     <  6:28860/142‑1 (MQ=255)
                                                                 acGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgcg                                                                                    >  5:28860/2‑142 (MQ=255)
                                                                         cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                       >  3:188736/1‑115 (MQ=255)
                                                                         cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                       <  4:188736/115‑1 (MQ=255)
                                                                         cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTg                                                                      <  1:153769/148‑1 (MQ=255)
                                                                          gCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGa                                                                     >  8:363491/1‑148 (MQ=255)
                                                                            cTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCAAGTACCACTTGTGGCAGAGGGGACCCATGCCCGGTGGAGTCTTTTGCCGACTATTCATTGAGCATCTGAATTCGCCAATGCGTACGTGCCGCGTCGGCAATAAACTGCcg                                                                   >  2:406293/1‑148 (MQ=255)
                                                                             tAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATTCCCGGTGGATTCTTTTGc                                                                                                                                  <  3:61601/84‑1 (MQ=255)
                                                                              aaTTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGTGAGCCATGCCCCGTGGATTCTTTTgac                                                                                                                                 >  4:61601/1‑82 (MQ=255)
                                                                              aaTTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATATACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGAcgcg                                                                 >  3:132264/1‑148 (MQ=255)
                                                                                ttcttcACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGTAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCa                                                               <  8:104808/148‑1 (MQ=255)
                                                                                      aCCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCTGTGGGGAGCCATGCCTGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTg                                                          <  5:374076/147‑1 (MQ=255)
                                                                                           ctcctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGa                                                      >  6:200241/1‑146 (MQ=255)
                                                                                           ctcctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGa                                                      <  5:200241/146‑1 (MQ=255)
                                                                                             cctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAg                                                  >  6:248758/1‑148 (MQ=255)
                                                                                              ctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGc                                                 <  6:293437/148‑1 (MQ=255)
                                                                                                    gACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCAc                                           >  2:4456/1‑148 (MQ=255)
                                                                                                           cccACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCCTGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGtt                                    >  7:104372/1‑148 (MQ=255)
                                                                                                                caTGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCg                                >  2:228526/1‑147 (MQ=255)
                                                                                                                  tGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCaa                             <  3:170016/148‑1 (MQ=255)
                                                                                                                    tACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAg                            <  5:248758/147‑1 (MQ=255)
                                                                                                                       cAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTa                        >  5:53821/1‑148 (MQ=255)
                                                                                                                           tgtgGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGt                                     >  7:193377/1‑131 (MQ=255)
                                                                                                                           tgtgGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGt                                     <  8:193377/131‑1 (MQ=255)
                                                                                                                                   tGGGGAGCCATGCCCGGTGGATTCTTTTGCCTACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGc                                                         <  6:136668/103‑1 (MQ=255)
                                                                                                                                   tGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGc                                                         >  5:136668/1‑103 (MQ=255)
                                                                                                                                    ggggAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGca            >  8:260081/1‑147 (MQ=255)
                                                                                                                                     gggAGCCATGCCTGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGcaca          <  8:128305/148‑1 (MQ=255)
                                                                                                                                     gggAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGcac           >  5:196436/1‑147 (MQ=255)
                                                                                                                                           cATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAg    >  7:366646/1‑148 (MQ=255)
                                                                                                                                            aTGCGCGGTGGCTTCTTTTGCCGACTGTTCATTCTGGCTCTCATGTCGCCAAGGCGAACGTGCCGCGGCGGCCATAAACTGCCGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCGGCAAGGCTACCGCCGCCTGCACAAGCCAGt   <  3:269578/148‑1 (MQ=255)
                                                                                                                                             tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                  <  4:134642/116‑1 (MQ=255)
                                                                                                                                             tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                  >  3:134642/1‑116 (MQ=255)
                                                                                                                                             tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGtt  >  5:185429/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTGATAGCGGCAAGCATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTT  >  NC_000913/1452114‑1452402

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: