Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,794,737 G→C 5.6% A587G (GCC→GGC)  fcs ← feruloyl‑CoA‑synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,794,7370GC 5.6% 232.0 / 4.4 72A587G (GCC→GGC) fcsferuloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/40);  new base C (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTGGCATCCTC  >  NC_002947/3794595‑3794885
                                                                                                                                               |                                                                                                                                                    
tGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCTAGCCCTTGTCGGTGATTTCACCGGCGtcgat                                                                                                                                               >  1:832088/1‑150 (MQ=255)
 gCACAGCCCGCCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatc                                                                                                                                              <  1:2297372/150‑1 (MQ=255)
   acaGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatcga                                                                                                                                            >  1:253953/1‑150 (MQ=255)
        ccGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggcg                                                                                                                                       >  2:1531186/1‑150 (MQ=255)
             cGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCg                                                                                                                                  <  1:1531186/150‑1 (MQ=255)
                 cTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCca                                                                                                                              >  1:1747351/1‑150 (MQ=255)
                  tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCcag                                                                                                                             >  2:533053/1‑150 (MQ=255)
                  tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCcag                                                                                                                             >  1:1522333/1‑150 (MQ=255)
                  tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCcag                                                                                                                             >  1:2113561/1‑150 (MQ=255)
                    aaGGCCGCAGCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCcagca                                                                                                                           >  2:1678807/1‑150 (MQ=255)
                    aaGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCgg                                                                                                                                                      <  2:223835/123‑1 (MQ=255)
                    aaGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCgg                                                                                                                                                      >  1:223835/1‑123 (MQ=255)
                    aaGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCcagca                                                                                                                           <  2:1522333/150‑1 (MQ=255)
                       gCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCg                                                                                                                        >  2:647605/1‑150 (MQ=255)
                        ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGa                                                                                                                       >  2:1088574/1‑150 (MQ=255)
                        ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGa                                                                                                                       >  1:2226552/1‑150 (MQ=255)
                        ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGa                                                                                                                       <  2:1622228/150‑1 (MQ=255)
                        ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGa                                                                                                                       <  1:619952/150‑1 (MQ=255)
                          gCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACa                                                                                                                     <  1:134167/150‑1 (MQ=255)
                           cACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGGGACAg                                                                                                                    <  1:2017218/150‑1 (MQ=255)
                              cTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCa                                                                                                                 >  2:931871/1‑150 (MQ=255)
                                   cGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGAGTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGa                                                                                                            <  2:2374462/150‑1 (MQ=255)
                                    gTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGAt                                                                                                           <  1:93693/150‑1 (MQ=255)
                                      gCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATAc                                                                                                         <  1:1678807/150‑1 (MQ=255)
                                      gCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCCCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGCCTCGGCCAGCAGCGACAGCCATTCGATAc                                                                                                         <  2:246757/150‑1 (MQ=255)
                                           gCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTg                                                                                                    <  1:931871/150‑1 (MQ=255)
                                              tCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCg                                                                                                 >  1:1827209/1‑150 (MQ=255)
                                              tCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCg                                                                                                 >  1:1826779/1‑150 (MQ=255)
                                                ggggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                                                                                <  2:1100011/117‑1 (MQ=255)
                                                ggggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                                                                                >  1:1100011/1‑117 (MQ=255)
                                                   gTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGCCGTCGATCGAcgccggtga                                                                                                                                   <  2:2339887/111‑4 (MQ=255)
                                                   gTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGCCGTCGATCGAcgccggtga                                                                                                                                   >  1:2339887/1‑108 (MQ=255)
                                                      ttCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggcg                                                                                                                                       >  2:1823579/1‑104 (MQ=255)
                                                      ttCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAcggcg                                                                                                                                       <  1:1823579/104‑1 (MQ=255)
                                                      ttCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTg                                                                                         <  2:2226552/150‑1 (MQ=255)
                                                       tCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTgg                                                                                        <  1:881203/150‑1 (MQ=255)
                                                                 aCAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGt                                                                              <  1:1586397/150‑1 (MQ=255)
                                                                  cAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGtt                                                                             <  1:387522/150‑1 (MQ=255)
                                                                     cgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                                                                          <  2:2336307/150‑1 (MQ=255)
                                                                      gcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAg                                                                         <  2:2195586/150‑1 (MQ=255)
                                                                         tCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCACACTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCCAgcgc                                                                      >  1:571931/1‑150 (MQ=255)
                                                                          cgcccTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCt                                                                     <  2:2180157/147‑1 (MQ=255)
                                                                            aCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTcc                                                                   <  2:641710/150‑1 (MQ=255)
                                                                            aCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTcc                                                                   <  1:1765025/150‑1 (MQ=255)
                                                                              cTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccag                                                                 >  2:87658/1‑150 (MQ=255)
                                                                                  tcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTg                                                                                         <  2:2174779/119‑1 (MQ=255)
                                                                                  tcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTg                                                                                         >  1:2174779/4‑122 (MQ=255)
                                                                                     gcgcgcCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCCGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCa                                                          <  2:1747351/150‑1 (MQ=255)
                                                                                       gcgcCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGc                                                        <  1:229930/150‑1 (MQ=255)
                                                                                         gccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGTTTCACGCGCTCCAGCCAGTCGGCGa                                                      >  1:1192612/1‑150 (MQ=255)
                                                                                          ccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGaa                                                     <  2:1529220/150‑1 (MQ=255)
                                                                                          ccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGaa                                                     <  2:253953/150‑1 (MQ=255)
                                                                                                  gca‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGCCGTCGATCGAcgccgg                                                                                                                                      <  2:586700/61‑1 (MQ=38)
                                                                                                  gca‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGCCGTCGATCGAcgccgg                                                                                                                                      >  1:586700/1‑61 (MQ=38)
                                                                                                           cgcTGATTGATCGGGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGgtg                                     <  2:1826779/150‑1 (MQ=255)
                                                                                                           cgcTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGgtg                                     <  2:1827209/150‑1 (MQ=255)
                                                                                                             cTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTc                                   <  1:3298/150‑1 (MQ=255)
                                                                                                              tgattgatCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTgg                                                                                                   >  2:833679/1‑85 (MQ=255)
                                                                                                              tgattgatCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTgg                                                                                                   <  1:833679/85‑1 (MQ=255)
                                                                                                              tgattgatCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCg                                  >  2:716761/1‑150 (MQ=255)
                                                                                                                attgatCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGtt                                <  2:1762580/150‑1 (MQ=255)
                                                                                                                attgatCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGtt                                <  1:647605/150‑1 (MQ=255)
                                                                                                                    atCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGcc                            >  1:2263490/1‑150 (MQ=255)
                                                                                                                     tCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCa                           <  1:87658/150‑1 (MQ=255)
                                                                                                                     tCGAGCACTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGccc                           >  1:141346/1‑149 (MQ=255)
                                                                                                                       gAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGc                         >  2:1962651/1‑150 (MQ=255)
                                                                                                                             ttGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCcgc                   >  2:1822107/1‑150 (MQ=255)
                                                                                                                             ttGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGCCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCcgc                   <  1:20528/150‑1 (MQ=255)
                                                                                                                                   gTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAAc                                                    <  2:596582/111‑1 (MQ=255)
                                                                                                                                   gTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAAc                                                    >  1:596582/1‑111 (MQ=255)
                                                                                                                                   gTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTAACTACAAGCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCg             >  1:149887/1‑150 (MQ=255)
                                                                                                                                      aTTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTg          <  1:1088574/150‑1 (MQ=255)
                                                                                                                                      aTTTCACCGGCGTCGATCGACGGCGGCCCGCACCGGAGCTACAGCCATTCGATAGGGCTGGCTTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTg          <  2:571931/150‑1 (MQ=255)
                                                                                                                                       tttCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTgg         <  1:647394/150‑1 (MQ=255)
                                                                                                                                              ggCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGACATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTGGCATCctc  >  2:384610/1‑150 (MQ=255)
                                                                                                                                               gCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc                                                                    >  1:1435103/1‑83 (MQ=255)
                                                                                                                                               gCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc                                                                    <  2:1435103/83‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCACAGCCCGGCCGCACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGCGAACCAACTGCGCACGGTGTCGTTGGCCAGCACCCGCGCATCGCTGGCATCCTC  >  NC_002947/3794595‑3794885

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