Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,670,473 T→A 7.0% intergenic (+62/‑23) metF → / → PP_4978 5,10‑methylenetetrahydrofolate reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,670,4730TA 7.0% 243.3 / 7.3 86intergenic (+62/‑23)metF/PP_49785,10‑methylenetetrahydrofolate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/38);  new base A (4/2);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

GGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGCTAGGCAGC  >  NC_002947/5670325‑5670618
                                                                                                                                                    |                                                                                                                                                 
ggcggcGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAt                                                                                                                                                   <  8:327838/149‑1 (MQ=255)
ggcggaGCGCCAGGCCTTCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAt                                                                                                                                                   <  6:113383/149‑1 (MQ=255)
   ggcgCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCa                                                                                                                                                <  8:282248/149‑1 (MQ=255)
         ccAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCg                                                                                                                                          <  3:340356/149‑1 (MQ=255)
          cAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCgg                                                                                                                                         <  7:333719/149‑1 (MQ=255)
          cAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCgg                                                                                                                                         >  5:1452/1‑149 (MQ=255)
             gCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGtt                                                                                                                                      <  7:100211/149‑1 (MQ=255)
                  cACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTaaa                                                                                                                                   <  2:149308/147‑1 (MQ=255)
                  cACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTaaa                                                                                                                                   >  2:327755/1‑147 (MQ=255)
                  cACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTaaa                                                                                                                                   <  1:187164/147‑1 (MQ=255)
                  cACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACa                                                                                                                                 <  1:191201/149‑1 (MQ=255)
                   aCTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAg                                                                                                                                <  2:11344/149‑1 (MQ=255)
                   aCTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAg                                                                                                                                <  3:219071/149‑1 (MQ=255)
                          taCCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtg                                                                                                                          <  8:192786/148‑1 (MQ=255)
                           aCCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcc                                                                                                                        >  6:197082/1‑149 (MQ=255)
                             cTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGTGTTAAACATGTCTtgcctg                                                                                                                      >  8:115167/1‑149 (MQ=255)
                             cTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctg                                                                                                                      >  5:219729/1‑149 (MQ=255)
                              ttGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctgc                                                                                                                     >  4:277675/1‑149 (MQ=255)
                              ttGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgcctgc                                                                                                                     >  6:206578/1‑149 (MQ=255)
                                  aCCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATg                                                                                                                 <  5:197082/149‑1 (MQ=255)
                                  aCCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATg                                                                                                                 <  8:178726/149‑1 (MQ=255)
                                      ggCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgcagc                                                                                                             <  1:327755/149‑1 (MQ=255)
                                      ggCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGCCTTAATCAtttcttaa                                                                                                                                        <  1:46550/122‑6 (MQ=25)
                                      ggCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGCCTTAATCAtttcttaa                                                                                                                                        >  2:46550/1‑117 (MQ=25)
                                       gCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgcagca                                                                                                            >  6:136727/1‑149 (MQ=255)
                                       gCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgcagca                                                                                                            >  6:64484/1‑149 (MQ=255)
                                           aaCCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACAcc                                                                                                        >  2:314057/1‑149 (MQ=255)
                                            aCCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACAccc                                                                                                       <  8:289579/149‑1 (MQ=255)
                                              cGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAtt                                                                                                                                                  >  8:318253/1‑104 (MQ=255)
                                              cGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAtt                                                                                                                                                  <  7:318253/104‑1 (MQ=255)
                                              cGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCaa                                                                                                     >  4:116156/1‑149 (MQ=255)
                                                      gCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgc                                                                                              <  1:314057/148‑1 (MQ=255)
                                                           cTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCca                                                                                         >  3:118641/1‑148 (MQ=255)
                                                             ggAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCAt                                                                                                                  >  5:20352/1‑121 (MQ=255)
                                                             ggAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCAt                                                                                                                  <  6:20352/121‑1 (MQ=255)
                                                              gAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCcac                                                                                        <  3:116156/146‑1 (MQ=255)
                                                              gAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGc                                                                                     >  4:151941/1‑149 (MQ=255)
                                                               aacaacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCt                                                                                    >  8:344289/1‑149 (MQ=255)
                                                                  aacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTg                                                                                  >  6:246664/1‑148 (MQ=255)
                                                                    ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGTCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGta                                                                                >  4:120120/1‑148 (MQ=255)
                                                                    ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGta                                                                                >  8:224104/1‑148 (MQ=255)
                                                                     cTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCtt                                                                                   >  7:296793/1‑144 (MQ=255)
                                                                     cTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCtt                                                                                   <  8:296793/144‑1 (MQ=255)
                                                                          gCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTgg                                                                         <  6:219729/149‑1 (MQ=255)
                                                                            tGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCt                                                                       >  6:292461/1‑149 (MQ=255)
                                                                              ccGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCttt                                                                     <  5:292461/149‑1 (MQ=255)
                                                                                gcgcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTg                                                                   <  5:136727/149‑1 (MQ=255)
                                                                                gcgcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTg                                                                   <  5:246664/149‑1 (MQ=255)
                                                                                  gcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGc                                                                 >  1:129275/1‑149 (MQ=255)
                                                                                    tGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGTCCACAGCTTGTATACCTgg                                                                         >  1:62440/1‑139 (MQ=255)
                                                                                    tGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGTCCACAGCTTGTATACCTgg                                                                         <  2:62440/139‑1 (MQ=255)
                                                                                           ttttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAcc                                                        >  4:76177/1‑149 (MQ=255)
                                                                                           ttttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAcc                                                        >  1:151605/1‑149 (MQ=255)
                                                                                           ttttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAATCATTTTAGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTgg                                                                         >  6:165977/1‑132 (MQ=255)
                                                                                                gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTAcc                                                    <  5:206578/148‑1 (MQ=255)
                                                                                                 cAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTg                                                                   >  3:250764/1‑132 (MQ=255)
                                                                                                 cAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTg                                                                   <  4:250764/132‑1 (MQ=255)
                                                                                                  aaaGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATc                                                 <  6:1452/149‑1 (MQ=255)
                                                                                                     gTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAgg                                              <  4:118641/149‑1 (MQ=255)
                                                                                                       cTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAgg                                              >  1:123009/1‑147 (MQ=255)
                                                                                                         gCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCTCACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGc                                          >  5:275751/1‑149 (MQ=255)
                                                                                                          cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGaa                                            <  3:76177/146‑1 (MQ=255)
                                                                                                           ttAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCt                                        >  4:8859/1‑149 (MQ=255)
                                                                                                             cgaaaaTGAATAAGGCTTTAGTCTATGACCAAAGCCTTAATCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGc                                      >  5:14308/4‑149 (MQ=37)
                                                                                                               gAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCtt                                    >  4:299920/1‑149 (MQ=255)
                                                                                                               gAAATGATTAAGGCTTTGGTCATAAACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGAACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGCCAACCTACCATCAGGAAGCCTGCt                                     >  6:31955/1‑148 (MQ=255)
                                                                                                                 aaTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGTCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                                  <  3:120120/149‑1 (MQ=255)
                                                                                                                 aaTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                                  <  5:64484/149‑1 (MQ=255)
                                                                                                                 aaTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                                  <  3:8859/149‑1 (MQ=255)
                                                                                                                    gATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAcc                                                        >  4:76385/1‑124 (MQ=255)
                                                                                                                    gATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAcc                                                        <  3:76385/124‑1 (MQ=255)
                                                                                                                     aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCa                               <  2:129275/148‑1 (MQ=255)
                                                                                                                     aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCa                               >  3:218930/1‑148 (MQ=255)
                                                                                                                     aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCAt                              >  8:73538/1‑149 (MQ=255)
                                                                                                                     aTTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCAt                              >  1:31479/1‑149 (MQ=255)
                                                                                                                      ttAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATTCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc                             >  8:198183/1‑149 (MQ=255)
                                                                                                                       tAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCg                            >  8:207572/1‑149 (MQ=255)
                                                                                                                         aGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATTCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGt                           <  7:198183/148‑1 (MQ=255)
                                                                                                                          ggCTTTAGTCTATGACCAAAGCCTTAATCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCa                                                >  4:64234/1‑126 (MQ=38)
                                                                                                                          ggCTTTAGTCTATGACCAAAGCCTTAATCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCa                                                <  3:64234/126‑1 (MQ=38)
                                                                                                                           gCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAg                        <  2:151605/149‑1 (MQ=255)
                                                                                                                                          tAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGaa                                            >  6:167604/1‑114 (MQ=255)
                                                                                                                                          tAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGc          >  4:349971/1‑148 (MQ=255)
                                                                                                                                          tAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGc          >  7:238157/1‑148 (MQ=255)
                                                                                                                                            aaGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACTATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGCTAg       >  4:206804/1‑149 (MQ=255)
                                                                                                                                            aaGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCa                                                >  8:212316/1‑108 (MQ=255)
                                                                                                                                            aaGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCa                                                <  7:212316/108‑1 (MQ=255)
                                                                                                                                             aGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGCTAgg      <  3:151941/149‑1 (MQ=255)
                                                                                                                                                cTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCATCAGCGCCCTCGCCCAGCTAGgcag   <  4:218930/149‑1 (MQ=255)
                                                                                                                                                 ttATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGCTAGgcagc  >  5:323632/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGCCCAGCTAGGCAGC  >  NC_002947/5670325‑5670618

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: