Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 784,333 G→A 11.2% R531R (CGC→CGT)  PP_0672 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947784,3330GA11.2% 165.7 / 7.8 62R531R (CGC→CGT) PP_0672sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (4/3);  total (31/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 = 1.00e+00

GTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAGCTCGGTGCGGT  >  NC_002947/784185‑784475
                                                                                                                                                    |                                                                                                                                              
gTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTgg                                                                                                                                                >  8:205430/1‑149 (MQ=255)
 ttGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGc                                                                                                                                               >  2:360666/1‑149 (MQ=255)
 ttGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGc                                                                                                                                               >  2:37799/1‑149 (MQ=255)
     tCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGTGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCgg                                                                                                                                             >  2:278952/1‑147 (MQ=255)
             aCGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCt                                                                                                                                   >  8:115357/1‑149 (MQ=255)
             aCGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCt                                                                                                                                   >  3:342037/1‑149 (MQ=255)
               gTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGTGGAAGCGCAGCTgg                                                                                                                                 <  7:254312/149‑1 (MQ=255)
               gTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTgg                                                                                                                                 <  6:387304/149‑1 (MQ=255)
                  gCGACTGGCTGGTCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCa                                                                                                                               <  2:2642/148‑1 (MQ=255)
                   cGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGt                                                                                                                             >  2:96761/1‑149 (MQ=255)
                         gctggcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCg                                                                                                                       <  1:140380/149‑1 (MQ=255)
                            ggcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGtt                                                                                                                    >  6:268139/1‑149 (MQ=255)
                             gcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTg                                                                                                                   >  2:34917/1‑149 (MQ=255)
                             gcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTg                                                                                                                   <  2:281958/149‑1 (MQ=255)
                              ccTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTgg                                                                                                                  >  2:24968/1‑149 (MQ=255)
                               cTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGt                                                                                                                 >  6:235445/1‑149 (MQ=255)
                               cTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGCCACCGCGACGGAAg                                                                                                                                          >  2:44335/1‑124 (MQ=255)
                               cTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGCCACCGCGACGGAAg                                                                                                                                          <  1:44335/124‑1 (MQ=255)
                                  gATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCa                                                                                                               <  8:333971/148‑1 (MQ=255)
                                  gATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCa                                                                                                               <  7:205430/148‑1 (MQ=255)
                                           ggTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGtt                                                                                                       >  5:74155/1‑147 (MQ=255)
                                           ggTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGtt                                                                                                       <  6:74155/147‑1 (MQ=255)
                                               acaTCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTcgcg                                                                                                 <  2:57245/149‑1 (MQ=255)
                                               acaTCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCATCGCGGTGGCGGAAGCGCa                                                                                                                                      <  7:52876/112‑1 (MQ=255)
                                               acaTCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCATCGCGGTGGCGGAAGCGCa                                                                                                                                      >  8:52876/1‑112 (MQ=255)
                                                   ccTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGt                                                                                             >  8:45778/1‑149 (MQ=255)
                                                         ccACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTggcgg                                                                                       <  6:17030/149‑1 (MQ=255)
                                                         ccACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGT‑GCGCTGTTggcggc                                                                                      >  6:312493/1‑149 (MQ=255)
                                                           acacCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCTGCGGTTGGTCATGACAGTTGTCGCGCTGTTGgcggcg                                                                                     >  2:181290/1‑149 (MQ=255)
                                                           acacCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGCCACCGCGACGGAAgctgg                                                                                                                                      <  1:271502/100‑4 (MQ=38)
                                                           acacCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGCCACCGCGACGGAAgctgg                                                                                                                                      >  2:271502/1‑97 (MQ=38)
                                                             acCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGgcggcg                                                                                     >  2:323192/1‑147 (MQ=255)
                                                                gatgatGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTg                                                                                <  4:342037/149‑1 (MQ=255)
                                                                  tgatgCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGCcgcg                                                                                                 >  5:252292/1‑130 (MQ=255)
                                                                  tgatgCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGCcgcg                                                                                                 <  6:252292/130‑1 (MQ=255)
                                                                    atgCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCt                                                                             <  7:45778/148‑1 (MQ=255)
                                                                     tgCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTgg                                                                           <  1:360666/149‑1 (MQ=255)
                                                                     tgCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTgg                                                                           <  1:37799/149‑1 (MQ=255)
                                                                             aCCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTgc                                                                   >  8:302279/1‑149 (MQ=255)
                                                                                   ccTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTc                                                             <  1:278952/149‑1 (MQ=255)
                                                                                        ggCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCgg                                                        <  8:25634/149‑1 (MQ=255)
                                                                                          cGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCgggg                                                      <  4:268657/149‑1 (MQ=255)
                                                                                              gCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCTGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTgc                                                   <  5:268139/148‑1 (MQ=255)
                                                                                                   tCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTgcagcag                                              >  1:380833/1‑148 (MQ=255)
                                                                                                     aGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGcagcagcag                                           <  1:323192/149‑1 (MQ=255)
                                                                                                        aCCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTc                                        >  7:346672/1‑149 (MQ=255)
                                                                                                          cTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTcc                                       >  5:147695/1‑148 (MQ=255)
                                                                                                              tCGCGGATGTCGTAGCAGCGCCCGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCGGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAg                                    >  3:201679/1‑147 (MQ=255)
                                                                                                                gcgGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGc                                >  6:166765/1‑149 (MQ=255)
                                                                                                                gcgGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGAGTCCTCGGGGTGCAGCAGCAGCTCCCAGaa                                  <  1:34917/147‑1 (MQ=255)
                                                                                                                gcgGATGTCCGAGCAGCGCCAGCGTCCGTCGCGGGGGCGGAAGCGCAGCTGGCAGTGCAGCGGCTGGGCATGACCCTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGaa                                  <  4:201679/147‑1 (MQ=255)
                                                                                                                       tCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGc                          <  1:24968/148‑1 (MQ=255)
                                                                                                                              agcGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCtgt                  <  7:302279/149‑1 (MQ=255)
                                                                                                                                 ctgcgCTTCCGCCACCGCGACGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTcc                                       <  2:343295/121‑1 (MQ=255)
                                                                                                                                 ctgcgCTTCCGCCACCGCGACGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTcc                                       >  1:343295/5‑125 (MQ=255)
                                                                                                                                  tgcgCTTCCGCCACCGCGACGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGc                >  6:277241/4‑147 (MQ=255)
                                                                                                                                   cAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAg             >  6:335081/1‑149 (MQ=255)
                                                                                                                                     gCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAGCt           <  2:380833/149‑1 (MQ=255)
                                                                                                                                      cTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGt                             >  4:286695/1‑130 (MQ=255)
                                                                                                                                      cTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGt                             <  3:286695/130‑1 (MQ=255)
                                                                                                                                      cTTCAGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGCCGGCGCAGTGCCAGGTAGAGTGCGCAGTCGTCGGGGTGCAGCAGCAGCTCCCAGACGCGGTCGCACATCTCTGCCAGCTc          >  8:336307/1‑149 (MQ=255)
                                                                                                                                           gTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAGCTCGGTGc     >  6:18441/1‑149 (MQ=255)
                                                                                                                                              gcgGTGGCGGAAGCGCAGCTGGCAGTGCAGCTGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAGCTCGGTGCGGt  <  8:346672/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTGGCGGCGCAGTGCCTGGTAGTGTGCGGCGTCCTCGGGGTGCAGCAGCAGCTCCCAGAAGCGGTCGCCCATCTGTGCCAGCTCGGTGCGGT  >  NC_002947/784185‑784475

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