Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,900,119 T→G 5.8% L176R (CTG→CGG)  PP_3443 → glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,900,1190TG 5.8% 201.4 / 2.8 68L176R (CTG→CGG) PP_3443glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (2/2);  total (34/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.14e-01

CAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTCG  >  NC_002947/3899973‑3900259
                                                                                                                                                  |                                                                                                                                               
cAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTgg                                                                                                                                               >  4:161753/1‑149 (MQ=255)
  aGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGc                                                                                                                                             <  4:319636/149‑1 (MQ=255)
     aCCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTg                                                                                                                                          <  5:42052/149‑1 (MQ=255)
      ccTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTgg                                                                                                                                         <  1:122989/149‑1 (MQ=255)
      ccTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTgg                                                                                                                                         >  3:183869/1‑149 (MQ=255)
        tGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGCTCCGCCGCGCCCCGCTGGGCGTGGcc                                                                                                                                       >  7:386955/1‑149 (MQ=255)
         gAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGccc                                                                                                                                      >  7:145581/1‑149 (MQ=255)
              aTTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCtgtg                                                                                                                                  >  6:329237/1‑148 (MQ=255)
               ttCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCa                                                                                                                                >  4:92783/1‑149 (MQ=255)
                  ggAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATgg                                                                                                                             >  5:324862/1‑149 (MQ=255)
                  ggAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCACGCCCAGCGGGGCGCGGc                                                                                                                                        >  1:255962/1‑138 (MQ=255)
                  ggAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGATCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCACGCCCAGCGGGGCGCGGc                                                                                                                                        <  2:255962/138‑1 (MQ=255)
                   gAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATgg‑‑g                                                                                                                          >  5:196842/1‑149 (MQ=255)
                   gAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATgg‑‑g                                                                                                                          >  4:139882/1‑149 (MQ=255)
                   gAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATgg‑‑g                                                                                                                          >  4:145765/1‑149 (MQ=255)
                          aaTTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTa                                                                                                                   >  7:147113/1‑148 (MQ=255)
                          aaTTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTAc                                                                                                                  <  6:97139/149‑1 (MQ=255)
                          aaTTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTAc                                                                                                                  <  3:92783/149‑1 (MQ=255)
                               gACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTa                                                                                                             >  6:48405/1‑149 (MQ=255)
                                      ccTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACt                                                                                                              <  5:52270/141‑1 (MQ=255)
                                      ccTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACt                                                                                                              >  6:52270/1‑141 (MQ=255)
                                           gACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAAc                                                                                                   >  2:99373/1‑147 (MQ=255)
                                           gACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACga                                                                                                 <  5:48405/149‑1 (MQ=255)
                                            aCTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACgag                                                                                                <  7:223812/149‑1 (MQ=255)
                                              tACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACGGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGTGTGGCCCTGTGCATGG‑‑GTCC‑CTACAACTACCCGCTGAACGAGAc                                                                                              <  1:287058/149‑1 (MQ=255)
                                                cATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCt                                                                                            >  1:222387/1‑149 (MQ=255)
                                                       gACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTcaccac                                                                                     <  8:147113/149‑1 (MQ=255)
                                                          acCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCt                                                                                  >  7:222862/1‑149 (MQ=255)
                                                              tCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATc                                                                              <  3:161753/149‑1 (MQ=255)
                                                               cAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATcc                                                                             >  8:147886/1‑149 (MQ=255)
                                                                aaCGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATccc                                                                            <  4:183869/149‑1 (MQ=255)
                                                                     ccTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgcgc                                                                       >  6:258517/1‑149 (MQ=255)
                                                                     ccTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgcgc                                                                       >  3:30001/1‑149 (MQ=255)
                                                                      cTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCt                                                                      >  5:328995/1‑149 (MQ=255)
                                                                          aGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCACGCCCAGCGGGGCGCgg                                                                                                                                         >  4:258791/1‑81 (MQ=255)
                                                                          aGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCACGCCCAGCGGGGCGCgg                                                                                                                                         <  3:258791/81‑1 (MQ=255)
                                                                                      gcctcTCCAGCCGCTTCGAGCTGGAGCAAGGCACACTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACCAGACCTTCACAACGCTCATACCGGCGCTGATCATGGGCAACACg                                                      >  6:413611/1‑149 (MQ=255)
                                                                                      gccgcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACg                                                      >  2:337483/1‑149 (MQ=255)
                                                                                      gccgcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACg                                                      <  5:258517/149‑1 (MQ=255)
                                                                                       ccgcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACgg                                                     >  3:369545/1‑149 (MQ=255)
                                                                                         ggtgcAGGCGCTTCGAGCTGGAGCAAGGGACCCTCGGCAAGATGCGGCGGGCCCAGCTGGGCGTGGCTCTGTGCATGG‑‑GGGC‑CTACAACTACCCGGTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   <  6:328995/145‑1 (MQ=255)
                                                                                         gcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   <  8:145581/149‑1 (MQ=255)
                                                                                         gcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   <  7:147886/149‑1 (MQ=255)
                                                                                         gcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   <  1:337483/149‑1 (MQ=255)
                                                                                          cTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgcgc                                                                       >  8:148504/1‑128 (MQ=255)
                                                                                          cTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgcgc                                                                       <  7:148504/128‑1 (MQ=255)
                                                                                              aGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGgtgt                                               >  2:184942/1‑148 (MQ=255)
                                                                                                   cTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGc                                         <  8:71742/149‑1 (MQ=255)
                                                                                                       gAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCgg                                      <  2:344172/148‑1 (MQ=255)
                                                                                                        aGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGcc                                    <  6:324862/149‑1 (MQ=255)
                                                                                                          cTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGCGGCC‑CTACAAATACCCGCTGAACGCG‑CCATCCCCACGCTCATCCCGGCGATGATAATGGGCACCACGGGGGTGTTCAAGCCGGcaa                                   >  1:120755/1‑147 (MQ=255)
                                                                                                          cTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggt                                                    <  3:221806/131‑1 (MQ=255)
                                                                                                          cTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggt                                                    >  4:221806/1‑131 (MQ=255)
                                                                                                             gAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCACCACGGTGGTGTTCAAGCCGGCCAAGtt                               >  8:163512/1‑149 (MQ=255)
                                                                                                               gCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCg                             <  8:386955/149‑1 (MQ=255)
                                                                                                                 aaGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGc                           <  6:196842/149‑1 (MQ=255)
                                                                                                                    gCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTa                        <  3:145765/149‑1 (MQ=255)
                                                                                                                    gCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTa                        <  3:139882/149‑1 (MQ=255)
                                                                                                                     cACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   >  4:23903/1‑121 (MQ=255)
                                                                                                                     cACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACggtg                                                   <  3:23903/121‑1 (MQ=255)
                                                                                                                         cTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCACGCCGGCCAAGTTCGGCGTActgct                   >  8:128913/1‑149 (MQ=255)
                                                                                                                            ggCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGAt                <  5:329237/149‑1 (MQ=255)
                                                                                                                              ccAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCCGCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATc               >  1:288805/1‑149 (MQ=255)
                                                                                                                              ccAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATcc              <  4:9268/149‑1 (MQ=255)
                                                                                                                                 gATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGAc           <  4:30001/149‑1 (MQ=255)
                                                                                                                                   tCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGAccc         >  3:167336/1‑149 (MQ=255)
                                                                                                                                         gcgcCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTc   <  1:184942/149‑1 (MQ=255)
                                                                                                                                          cgcCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTCg  <  8:222862/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGG‑‑GGCC‑CTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTCG  >  NC_002947/3899973‑3900259

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