Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 405,012 C→A 7.8% E350* (GAG→TAG)  PP_0337 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947405,0120CA 7.8% 189.6 / 2.4 60E350* (GAG→TAG) PP_0337sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCTGAGGTTCTCCAGTGCCGGTAGCTCACCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTTGGCTGGCCGCTTTGCTCCAGGC  >  NC_002947/404866‑405158
                                                                                                                                                  |                                                                                                                                                  
gCTGAGGTTCTCCAGTGCCGGTAGCTCACCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGAc                                                                                                                                                 >  1:783453/1‑150 (MQ=255)
  tGAGGTTCTCCAGTGCCGGTAGCTCACCCAGCTGGCAGCCCTGCACTTCGTCGGTACTTTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGGCGAGGTTCACCAGCAGGAAGGCACTGGTGCTCTGCTCGACCg                                                                                                                                               >  2:567499/1‑150 (MQ=255)
         ctcCAGTGCCGGTAGCTCACCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTg                                                                                                                                        >  2:1024228/1‑150 (MQ=255)
                           aCCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCt                                                                                                                      >  2:1012418/1‑150 (MQ=255)
                           aCCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGAGTGATGGAGGTGAACATGTCTATGACGTATTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCt                                                                                                                      >  1:1121835/1‑150 (MQ=255)
                              cAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACa                                                                                                                                   >  1:730711/1‑134 (MQ=255)
                              cAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACa                                                                                                                                   <  2:730711/134‑1 (MQ=255)
                                   ggcggcCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCgg                                                                                                              <  2:149264/150‑1 (MQ=255)
                                          cTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCcgtcg                                                                                                       >  2:260431/1‑150 (MQ=255)
                                             cACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCgttt                                                                                                    <  1:1012418/150‑1 (MQ=255)
                                             cACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCgttt                                                                                                    <  1:367711/150‑1 (MQ=255)
                                               cTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCgtttgt                                                                                                  <  1:260431/150‑1 (MQ=255)
                                                ttCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCgtttgtg                                                                                                 <  2:783453/150‑1 (MQ=255)
                                                  cgtcgGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACCCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTAGACCGCACGCTGGAACAGGTGCAGGGCGCTTGCCGCGGCCCGTCgtttgtggt                                                                                               <  2:408861/150‑1 (MQ=255)
                                                    tcgGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCgtttgtggttt                                                                                             >  2:68235/1‑150 (MQ=255)
                                                           tgtactgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCgg                                                                                                              <  1:301878/126‑1 (MQ=255)
                                                           tgtactgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCgg                                                                                                              >  2:301878/1‑126 (MQ=255)
                                                            gtactgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCAc                                                                                     <  1:1095369/150‑1 (MQ=255)
                                                             tcctgGGTGGTGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACt                                                                                    >  2:638069/3‑150 (MQ=255)
                                                             tactgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACt                                                                                    >  2:440159/1‑150 (MQ=255)
                                                              actgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACtt                                                                                   <  1:440159/150‑1 (MQ=255)
                                                                tgGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTgg                                                                                 <  1:1024228/150‑1 (MQ=255)
                                                                           ggTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTc                                                                      >  1:1052243/1‑150 (MQ=255)
                                                                              gAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGc                                                                   <  1:638069/150‑1 (MQ=255)
                                                                                  cTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGc                                                               >  1:422867/1‑150 (MQ=255)
                                                                                           aCGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCa                                                      >  1:1084367/1‑150 (MQ=255)
                                                                                           aCGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTACTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCa                                                      >  1:1151553/1‑150 (MQ=255)
                                                                                             gTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCaaa                                                    >  1:620979/1‑150 (MQ=255)
                                                                                                 tCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTg                                                                         >  1:547431/1‑125 (MQ=255)
                                                                                                 tCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTg                                                                         <  2:547431/125‑1 (MQ=255)
                                                                                                 tCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCa                                                >  1:80387/1‑150 (MQ=255)
                                                                                                  ccccccGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCaa                                               <  2:1201589/147‑1 (MQ=255)
                                                                                                   cACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAAt                                              >  2:1102058/1‑150 (MQ=255)
                                                                                                    aCCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCaacc                                             <  2:1181218/150‑3 (MQ=255)
                                                                                                    aCCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAAtc                                             <  1:1102058/150‑1 (MQ=255)
                                                                                                      ccGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCa                                                            <  1:1118052/133‑1 (MQ=255)
                                                                                                      ccGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCa                                                            >  2:1118052/1‑133 (MQ=255)
                                                                                                      ccGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAAtctc                                           >  2:590517/1‑150 (MQ=255)
                                                                                                       cGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCg                                          <  2:422867/150‑1 (MQ=255)
                                                                                                         gCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTc                                        <  2:745258/150‑1 (MQ=255)
                                                                                                            gTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCACAGGCAATCTCGTCCGa                                     >  2:292154/1‑150 (MQ=255)
                                                                                                             tCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGc                                                 <  1:670952/137‑1 (MQ=255)
                                                                                                             tCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGc                                                 >  2:670952/1‑137 (MQ=255)
                                                                                                             tCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGAc                                    >  2:329858/1‑150 (MQ=255)
                                                                                                               gcgGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTg                                  >  1:992584/1‑150 (MQ=255)
                                                                                                                cgGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGc                                 <  1:1158578/150‑1 (MQ=255)
                                                                                                                 ggTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGcc                                >  1:69402/1‑150 (MQ=255)
                                                                                                                  gTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGAAAAAGTGGGGGAACTGGTCGGCCAGTTCACCGTCAAAGGCAATCTCGTCCGACTGCCa                               <  1:482220/150‑1 (MQ=255)
                                                                                                                       ccAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGACAATCTCGTCCGACTGCCAGGccc                          >  1:979302/1‑150 (MQ=255)
                                                                                                                          gcagGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtg                       <  2:1121835/150‑1 (MQ=255)
                                                                                                                           cagGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtgg                      <  2:1084367/150‑1 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGccc                          >  1:918927/1‑145 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGccc                          <  2:918927/145‑1 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTAGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtgg                      <  2:1251172/149‑1 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTAGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtgg                      <  2:1250963/149‑1 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTAGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtgg                      >  1:1251172/1‑149 (MQ=255)
                                                                                                                            agGAAGGCACTGGCGCTCTGCTAGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtgg                      >  1:1250963/1‑149 (MQ=255)
                                                                                                                              gAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTtggctg                   >  2:1021819/1‑150 (MQ=255)
                                                                                                                                     cTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTCCTTGCGTGTGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTTGGCTGGCCGCtt            <  1:1119620/150‑1 (MQ=255)
                                                                                                                                               gCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTTGGCTGGCCGCTTTGCTCCAGGc  >  1:812322/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGAGGTTCTCCAGTGCCGGTAGCTCACCCAGCTGGCGGCCCTGCACTTCGTCGGTACTGTACTGGGTGATGGCGGTGAAGCTGGGGTTGACGTACTCCACCCGGCCGTCGCGGTTCACCAGCAGGAAGGCACTGGCGCTCTGCTCGACCGCACGCTGGAACAGGTGCAGGGCGCTGGCCGCGGCCCGTCGTTTGTGGTTTGTGATCACTTGGGCGAACTGGTCGGCCAGCTCACCGGCAAAGGCAATCTCGTCCGACTGCCAGGCCCTTGGCTGGCCGCTTTGCTCCAGGC  >  NC_002947/404866‑405158

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