Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,733,268 C→A 6.7% R56S (CGT→AGT)  PP_2394 → leucine‑rich repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,733,2680CA 6.7% 192.0 / 3.3 60R56S (CGT→AGT) PP_2394leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base A (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AACCATCGGCTGCGGTCGATACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTCACACACGA  >  NC_002947/2733122‑2733415
                                                                                                                                                  |                                                                                                                                                   
aaCCATCGGCTGCGGTCGATACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGt                                                                                                                                                   >  5:134197/1‑149 (MQ=255)
            cGGTCGATACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCACGGCAGGccagca                                                                                                                                        <  6:135068/148‑6 (MQ=255)
               tCGATACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGt                                                                                                                                      >  6:134367/1‑147 (MQ=255)
                  aTACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTg                                                                                                                                          >  7:146129/1‑140 (MQ=255)
                  aTACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTg                                                                                                                                          <  8:146129/140‑1 (MQ=255)
                     cccTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCt                                                                                                                               <  6:259511/148‑1 (MQ=255)
                      ccTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCCGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTgg                                                                                                                             >  1:273542/1‑149 (MQ=255)
                       cTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCAGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGgc                                                                                                                            >  2:271875/1‑149 (MQ=255)
                         tGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAgcg                                                                                                                                 >  6:333399/1‑142 (MQ=255)
                         tGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAgcg                                                                                                                                 <  5:333399/142‑1 (MQ=255)
                          gCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGACGCGCTGCCTGAAATGCTGGCCTGCAGTGCAGGCTTGGCGTCAGCGCTGgcgag                                                                                                                         >  6:273539/1‑147 (MQ=255)
                           ccATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCgg                                                                                                                        <  3:165800/149‑1 (MQ=255)
                             aTGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCt                                                                                                                      <  1:102602/149‑1 (MQ=255)
                             aTGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCGGCCGTGCAGGCTTGGCGGCAGCGCTGGCGCGGCt                                                                                                                      <  8:228758/149‑1 (MQ=255)
                                    aaGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAAc                                                                                                               <  4:52010/149‑1 (MQ=255)
                                     aGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAAcc                                                                                                              <  4:80793/149‑1 (MQ=255)
                                        aTGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTg                                                                                                            >  1:74064/1‑148 (MQ=255)
                                        aTGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCGGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTg                                                                                                            >  1:262354/1‑148 (MQ=255)
                                           ccAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAAt                                                                                                         <  6:134197/148‑1 (MQ=255)
                                           ccAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAAt                                                                                                         <  4:182545/148‑1 (MQ=255)
                                           ccAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCAGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAAt                                                                                                         <  1:271875/148‑1 (MQ=255)
                                                  cACTACATCGCACAGCGGTTGCCGGCCTTGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCTCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGc                                                                                                 <  4:166718/149‑1 (MQ=255)
                                                     tACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCg                                                                                              >  4:241346/1‑149 (MQ=255)
                                                       cATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGCACGGCAGGCCAGCAGTTCAgg                                                                                                                                              <  4:273290/99‑1 (MQ=38)
                                                       cATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGCACGGCAGGCCAGCAGTTCAgg                                                                                                                                              >  3:273290/1‑99 (MQ=38)
                                                       cATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCg                                                                                              <  6:6111/147‑1 (MQ=255)
                                                           gCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCt                                                                                         <  2:74064/148‑1 (MQ=255)
                                                            cacaGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAATCTGAATGCATTTGCCAGGCCCTtg                                                                                       >  4:81500/1‑149 (MQ=255)
                                                            cacaGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTtg                                                                                       >  7:69663/1‑149 (MQ=255)
                                                                 cGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTgca                                                                                  <  7:27499/149‑1 (MQ=255)
                                                                 cGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTgca                                                                                  <  5:160649/149‑1 (MQ=255)
                                                                   gTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTgcag                                                                                 >  3:44922/1‑148 (MQ=255)
                                                                        cGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAgcg                                                                           >  3:117542/1‑149 (MQ=255)
                                                                         ggCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAg                                                                             <  3:314367/146‑1 (MQ=255)
                                                                         ggCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAg                                                                             >  4:314367/1‑146 (MQ=255)
                                                                               gATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTg                                                                    >  6:75993/1‑149 (MQ=255)
                                                                                  gAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGcc                                                                  >  1:225924/1‑148 (MQ=255)
                                                                                            aGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCt                                                       >  8:233466/1‑149 (MQ=255)
                                                                                               aTGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGa                                                    <  4:117542/149‑1 (MQ=255)
                                                                                                tGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGAt                                                   >  5:88879/1‑149 (MQ=255)
                                                                                                  gCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGt                                                 <  3:241346/149‑1 (MQ=255)
                                                                                                     gAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGa                                              >  3:116582/1‑149 (MQ=255)
                                                                                                      aaTTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGAt                                             >  5:211552/1‑149 (MQ=255)
                                                                                                       aTTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAAt                                                                                                         <  8:167136/88‑1 (MQ=255)
                                                                                                       aTTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAAt                                                                                                         >  7:167136/1‑88 (MQ=255)
                                                                                                           aCCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCt                                        >  5:189825/1‑149 (MQ=255)
                                                                                                           aCCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCt                                        <  5:134367/149‑1 (MQ=255)
                                                                                                             ccTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGt                                      <  2:225924/149‑1 (MQ=255)
                                                                                                                         ggCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGt                          >  6:284407/1‑149 (MQ=255)
                                                                                                                                  tGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGc                              <  2:217517/136‑1 (MQ=255)
                                                                                                                                  tGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGc                              >  1:217517/1‑136 (MQ=255)
                                                                                                                                  tGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCa                            >  2:6919/1‑138 (MQ=255)
                                                                                                                                  tGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCa                            <  1:6919/138‑1 (MQ=255)
                                                                                                                                  tGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCt                 >  3:98097/1‑149 (MQ=255)
                                                                                                                                    aaCTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCt                 <  6:189825/147‑1 (MQ=255)
                                                                                                                                     aCTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTc              <  8:69663/149‑1 (MQ=255)
                                                                                                                                       tgctgGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTccc            >  7:232303/1‑149 (MQ=255)
                                                                                                                                         ctgGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCa                            >  2:136769/1‑131 (MQ=255)
                                                                                                                                         ctgGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCa                            <  1:136769/131‑1 (MQ=255)
                                                                                                                                         ctgGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTc          <  8:232303/149‑1 (MQ=255)
                                                                                                                                                 ccGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTCACACACg   <  6:211552/148‑1 (MQ=255)
                                                                                                                                                 ccGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTCACACACGa  <  5:75993/149‑1 (MQ=255)
                                                                                                                                                 ccGTGCAGGCTTGGCGGCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTCACACACg   <  2:273542/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACCATCGGCTGCGGTCGATACCCTTGCCATGGCACAAGCATGCCAGGACCACTACATCGCACAGCGGTTGCCGGCCTGGATGAAGCGGCTGAGCATGGCCGAATTCACCCTGCTGAGCGAGGCGCTGCCTGAACTGCTGGCCTGCCGTGCAGGCTTGGCGTCAGCGCTGGCGCGGCTAAGCAACCTGAATGCATTTGCCCGGCCCTTGCTGCAGCAAGCGCTCGGTGCCCACGGCGACCTTGATGTCGATCGCCTGTACTTTCGCCAGTGGTACACCTTTCCCTCACACACGA  >  NC_002947/2733122‑2733415

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