Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,122,363 C→G 5.9% G20A (GGC→GCC)  PP_2739 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,122,3630CG 5.9% 215.0 / 2.6 67G20A (GGC→GCC) PP_2739sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/2);  total (33/34)
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.70e-01

GATGCCCGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAAGGAG  >  NC_002947/3122215‑3122505
                                                                                                                                                     |                                                                                                                                              
gATGCCCGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGc                                                                                                                                                >  2:14636/1‑149 (MQ=255)
  tGCCCGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCa                                                                                                                                              <  6:126096/149‑1 (MQ=255)
     ccGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCg                                                                                                                                           >  1:288604/1‑149 (MQ=255)
        gCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTg                                                                                                                                        <  1:14636/149‑1 (MQ=255)
                     cTCGCCCTGAAGAACACGGCAGTCCTGATAGCGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCa                                                                                                                           >  8:80066/1‑149 (MQ=255)
                           cTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGtt                                                                                                                     >  5:155716/1‑149 (MQ=255)
                            tGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTg                                                                                                                    >  2:194667/1‑149 (MQ=255)
                            tGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTg                                                                                                                    >  2:12497/1‑149 (MQ=255)
                              aagaaAACGGCAGTCCTGATAGAGAATATCGTCGGCGGAGTAGCCGGTCAGGCGCTCGAAGGCCGGGGTGACGGAGATAAGGATGCTCTCATTGGCTTCCTGATCGGCGACAACGATGCCATCGTTGGAATGTTGGACCATCAGTTGct                                                                                                                  >  6:214219/1‑148 (MQ=255)
                              aagaaAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCa                                                                                                                  <  6:155716/149‑1 (MQ=255)
                                 aaaaCGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGa                                                                                                               <  1:260683/149‑1 (MQ=255)
                                   aaCGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGt                                                                                                             <  6:204207/149‑1 (MQ=255)
                                       gCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGGCATCGTTgtcg                                                                                                                                     >  7:105297/1‑118 (MQ=255)
                                       gCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGGCATCGTTgtcg                                                                                                                                     <  8:105297/121‑4 (MQ=255)
                                         aGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCg                                                                                                       <  2:288604/149‑1 (MQ=255)
                                          gTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGt                                                                                                      <  7:233717/149‑1 (MQ=255)
                                                aTAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCa                                                                                                >  3:210389/1‑149 (MQ=255)
                                                 tAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGGGTGGTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTAGACCATCAGTTGCAGGAGTTTTGCGTTTAtaa                                                                                                >  1:278549/1‑147 (MQ=255)
                                                 tAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGt                                                                                                      <  4:69387/142‑1 (MQ=255)
                                                 tAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGt                                                                                                      >  3:69387/1‑142 (MQ=255)
                                                   gagaATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATgg                                                                                             <  1:12497/149‑1 (MQ=255)
                                                    agaATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATggg                                                                                            >  5:303069/1‑149 (MQ=255)
                                                        tatCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATg                                                                                              <  8:179481/143‑1 (MQ=255)
                                                        tatCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATg                                                                                              >  7:179481/1‑143 (MQ=255)
                                                            gtcgGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCAcc                                                                                     >  2:90817/1‑148 (MQ=255)
                                                                 cgcAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAgg                                                                               <  4:210389/149‑1 (MQ=255)
                                                                    aGTAGCCGGTCAGGCGCTAGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCAt                                                                            >  2:3852/1‑149 (MQ=255)
                                                                               aGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGa                                                                                                                              >  2:34762/1‑88 (MQ=255)
                                                                               aGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGa                                                                                                                              <  1:34762/88‑1 (MQ=255)
                                                                                 gcgcTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGAtg                                                               <  1:90817/149‑1 (MQ=255)
                                                                                     tCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCg                                                                  <  2:180003/142‑1 (MQ=255)
                                                                                     tCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCg                                                                  >  1:180003/1‑142 (MQ=255)
                                                                                       gAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTgcgc                                                          <  6:303069/148‑1 (MQ=255)
                                                                                                    aCGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTg                                                                                                                                        >  7:129335/1‑57 (MQ=255)
                                                                                                    aCGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTg                                                                                                                                        <  8:129335/57‑1 (MQ=255)
                                                                                                      gTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCa                                                                                                                        >  6:108979/1‑71 (MQ=255)
                                                                                                      gTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCa                                                                                                                        <  5:108979/71‑1 (MQ=255)
                                                                                                      gTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCa                                          <  2:278549/149‑1 (MQ=255)
                                                                                                        aGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAg                                         >  6:287718/1‑148 (MQ=255)
                                                                                                          aTAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTg                                            >  1:306779/1‑143 (MQ=255)
                                                                                                          aTAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTg                                            <  2:306779/143‑1 (MQ=255)
                                                                                                           tAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGc                                     >  1:44700/1‑149 (MQ=255)
                                                                                                               ggATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCttt                                  >  2:192239/1‑149 (MQ=255)
                                                                                                                  tcctctCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTg                                <  1:122960/146‑1 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGAt                                                                >  5:188835/1‑113 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGAt                                                                <  6:188835/113‑1 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAg                                                     >  4:76609/1‑124 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTg                                >  2:122960/1‑145 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcgc                            >  3:83801/1‑149 (MQ=255)
                                                                                                                     tctcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcgc                            >  7:147878/1‑149 (MQ=255)
                                                                                                                       tcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCa                                          >  4:129969/1‑133 (MQ=255)
                                                                                                                       tcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCa                                          <  3:129969/133‑1 (MQ=255)
                                                                                                                       tcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcgc                            <  8:147878/147‑1 (MQ=255)
                                                                                                                       tcATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcgc                            <  2:44700/147‑1 (MQ=255)
                                                                                                                           tGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttg                      <  7:80066/149‑1 (MQ=255)
                                                                                                                               atCCTGCTCGGCGACAACGATGCCATCGTTGGTATGTTCGACCATCAGTTGCAGGAGGTTTGCGTTGATCTTGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGc                  <  4:74407/148‑1 (MQ=255)
                                                                                                                                  cTGCTCGGAGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCat                <  5:214219/148‑1 (MQ=255)
                                                                                                                                     cTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCACAAGCTTTTGGCGCCGTTTGTTGCatgat             <  1:194667/148‑1 (MQ=255)
                                                                                                                                      tCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCatgatg            >  8:149306/1‑148 (MQ=255)
                                                                                                                                        ggCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTcc                                               >  4:75857/1‑111 (MQ=255)
                                                                                                                                        ggCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTcc                                               <  3:75857/111‑1 (MQ=255)
                                                                                                                                        ggCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGaa          <  1:192239/148‑1 (MQ=255)
                                                                                                                                         gCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAAc         >  8:324940/1‑148 (MQ=255)
                                                                                                                                          ttccAACGATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCAtgatga           <  5:205490/142‑1 (MQ=255)
                                                                                                                                          ttccAACGATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCAtgatga           >  6:205490/4‑145 (MQ=255)
                                                                                                                                              aaCGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAag     >  4:149976/1‑147 (MQ=255)
                                                                                                                                               aCGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAaggag  <  4:83801/149‑1 (MQ=255)
                                                                                                                                               aCGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAaggag  <  5:321344/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GATGCCCGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCGTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAAGGAG  >  NC_002947/3122215‑3122505

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